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ABSTRACT 

A conceptual design of a direct Brayton-cycle marine power plant is presented. The 
design is a modification of the commercial MGR-GT, as proposed by James Staudt, sized to 
produce 40,000 shaft horsepower (SHP) and 5 MW of ship service electrical power. 

The requirements of a shipboard power plant are discussed and the design changes that 
must be made to the components of a commercial power plant in order to fit them into the 
demanding environment of a ship at sea are detailed. 

The final design consists of an 80-MWth passively-safe pebble bed reactor with an out- 
let temperature of 850*C. The reactor powers two separate closed Brayton cycle power con- 
version loops operating at a compressor discharge pressure of 8.2 MPa. Other features of the 
system are compact highly efficient heat exchangers, an advanced integrated electric 
propulsion system using solid state power converters and frequency changers, magnetic bear- 
ings, and high speed generators, both helium and water cooled. These features combined to 
produce system efficiencies exceeding 45%. 

Results show that the use of a direct cycle electric drive power plant for ship propulsion 
is an attractive alternative. The heat exchangers and rotating machinery can be made com- 
pact and light, but the direct use of the commercial reactor core is not viable because of 

. . J , . . 

weight and size. A small, fast spectrum core is one potential solution for the heat source if it 
can be produced within the constraints of passive safety. Also further work needs to be done 
to optimize the turbine-compressor-generator system. 

Thesis Supervisor: Dr. Lawrence M. Lidsky, 

Title; Professor of Nuclear Engineering. 
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Chapter 1 Introduction and Background 

1.1 Introduction. 

This paper investigates the feasibility of using a modular gas-cooled closed Brayton- 
cycle gas turbine reactor as a power plant for naval vessels by developing a self-consistent 
prototype design. Although this concept has been studied for over twenty years, recent 
advances in fuel design, materials, components, an ever increasing amount of gas reactor 
experience, plus recent political and social attitude changes have shown that this system war- 
rants renewed attention. 

1.1.1 Previous work. 

The idea of using a closed Brayton cycle with a gas cooled reactor as the heat source is 
not a new idea. The closed-cycle gas turbine was first introduced in 1936 [8] and the idea of 
coupling it with a gas-cooled reactor was introduced soon after. It was considered by the 
navy as early as the early fifties, when the navy began design work on the USS Nautilus 
(SSN-571). The Nautilus was eventually powered by a pressurized water reactor. In 1974 the 
Westinghouse Astronuclear Laboratory (WANL) presented a design concept for a 140,000 
SHP closed Brayton cycle plant. 

Figure 1-1 shows the WANL concept. The concept used a closed Brayton cycle with 
helium as the working fluid. The core is a graphite moderated, epithermal spectrum reactor, 
using TRISO fuel particles in extruded graphite fuel elements. The fuel is highly enriched 

The containment is shaped in an inverted ’T’ with two sections. The upper section 
contains the reactor core, control drums and the primary shield. The lower section contains 
two power conversion loops, each consisting of a turbine - compressor - heat exchanger 
package coupled to a superconducting generator. It also contains helium storage bottles, 
emergency cooling system, a fission product cleanup system, support stmcture, and second- 
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Figure 1-1. Isometric view of the WANL concept 140,000 HP Qosed-Brayton cycle 
power plant. [4] 

ary shielding. 

The priority of the WANL design was compacmess. The entire unit is only 32 feet 
high, 18 feet wide, and 34 feet long. To achieve this compact design expensive shielding 
materials, and superconducting generators were used to arrange the plant as tightly as possi- 
ble. With the compact shield and the highly power dense core, the plant cannot passively 
dissipate decay heat. This requires a quick reaction emergency cooling system to prevent 
core damage upon loss of cooling. Since water is not compatible with core materials, this 
emergency cooling system has to be gas based. Also, maintenance would be difficult 
because everything is packed so tightly together. Due to the high technical risk (along with 
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other factors) the WANL design was never funded, so development did not proceed. 

During the same time frame WANL also designed a smaller (25,000 shaft horsepower) 
plant concept based on the NERVA reactor core. This concept is shown in Figure 1-2. The 
system is basically the same as the larger module except the containment arranged longitudi- 
nally instead of vertically. The system is direct drive with a single power conversion loop 
consisting of a turbine - compressor - heat exchanger package coupled to a reduction gear 
and propulsion shaft. Again, funding and development did not proceed. 



(MERGENCY COOlINO 




Figure 1-2. Isometric view of the WANL concept 25,000 HP Qosed-Brayton cycle 
power plant [4] 

More recently, design studies of an MGR Brayton cycle power plant for commercial 
power generation have been conducted at MIT and in the United States and Germany. The 
Germans have acquired a wealth of experience in large closed-cycle gas turbines, having 
operated fossil fueled units since 1956. The largest, Oberhausen 2, operates at 50 MW(e) 
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and is a test facility for designing a large HGTR-GT. The German nuclear power industry 
has also operated large pebble bed reactors at the temperatures necessary to make the gas tur- 
bine economical. 

In the United States there is only one commercial gas cooled reactor in operation, the 
330-MWe Fort St. Vrain power station in Colorado. It is a nuclear steam system that has 
suffered many problems over its history. Many of the problems arising from water ingress 
into the helium coolant loop from the gas circulator bearings. 

The General Atomic Company has conducted an extensive design study of a large 
direct Brayton cycle plant. Their design focused on a plant with a power level ranging from 
800 to 1200 MWe with from two to four power conversion loops. All plant components 
were to be contained in a prestressed concrete reactor vessel. The program was terminated 
because of economic reason arising from the many technical problems which arose during 
the study. 

The MIT design is for a smaller, modular plant. (See sect. 1.2) It is the MIT design 
that the design for this reactor wUl be based. 

1.1.2 Marine Reactors. 

Nuclear power has long been recognized as an ideal power source for naval vessels. 
Although there are many advantages, the principle benefit is elimination of refueling or bat- 
tery charging. In general, surface ships carry enough fuel for five to six thousand miles at 
their most economical speed. Conventional submarines use a diesel engine to travel on the 
surface and batteries or some other air independent propulsion system when submerged. 
Non-nuclear submarines must carry both fuel and batteries, a double burden. The range of 
diesel submarines on the surface is about the same as a surface vessel, but current battery 
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technology limits the submerged endurance to a few hours. In order to achieve even this 
modest range, five to ten percent of total ship displacement is fuel and two to three percent of 
a submarine’s weight is battery. 

A nuclear power plant gives both surface ships and submarines virtually unlimited 
range, and it frees a large amount of space and weight that previously had to be devoted to 
fuel. Without the fuel load to consider, ships can be made smaller, or the same size ship 
could carry a larger payload. 

Nuclear reactors have been used in marine propulsion since the USS Nautilus was 
launched in 1954. Since that day almost all submarines have been nuclear powered along 
with many large surface combatants. Smaller ships, such as frigates and destroyers, have not 
been nuclear powered. The pressurized water reactor is too big, too heavy, and too expensive 
to justify installation in the smaller ships. 

In this country, the pressurized water reactor coupled to a Rankine power conversion 
cycle has been virtually the only nuclear power plant used to power ships. Although the 
PWR has been a proven, reliable power plant at sea for over thirty years, it suffers from the 
following problems which prevent nuclear power from reaching its fuU potential. 

• Low Thermal Efficiency. Because of the difficulty of significantly superheating a 
PWR, Thermal efficiency is limited to around 30%. This low efficiency results in a 
great deal of latent heat being rejected to the sea. Besides the obvious energy waste, 
the heat wake generated by the waste heat could increase the detectability of the ship. 
To be fair, all thermal power plants reject most of their generated heat, but higher effi- 
ciency means more usable power for the same heat loss. 

• Large Component Size. Due to the relatively low energy density of saturated steam, 
machinery, piping, turbines, and other steam system components are much larger than 
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in conventional superheated steam plant designs. This is a weight and volume penalty 
that reduces some of the benefits over fossil fueled ships gained by eliminating the 
fuel load. 

Lack of Inherent Safety. The U.S. nuclear power program has an undeniably excel- 
lent safety record. This is due to superb operator training, well engineered systems, 
and a management system that treats safety as the most important aspect of the 
operation. This is engineered safety, however, not inherent safety. The systems to 
ensure the engineered safety of the naval PWR can be five to ten percent of the reactor 
plant weight. An inherently safe system can possibly save some of this weight and 
space. All of the safety system weight and volume cannot be recovered since 
inherently safe nuclear systems generally are larger and heavier than similar systems 
in a standard reactor. 

Heavy Weight Concentration. The PWR is a high density system compared to non- 
nuclear ship systems. Heavy shielding, large pumps, thick pressure vessels and pip- 
ing, combine to produce the heaviest single ship system besides the hull. This 
generally requires that it be placed low and in the center of the ship. This increases 
shielding requirements, and makes the ship more difficult to arrange. 

High Maintenance Costs. Radiological controls, safety graded systems, and inten- 
sive training makes maintenance expensive and time consuming. 

An inherendy safe closed Brayton cycle power plant has the potential to correct all of 



the above deficiencies. 
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1.1.3 Marine Environment. 

There are many design differences between a reactor for ship propulsion and commer- 
cial power reactor used by utilities. Most of these differences are a result of their respective 
environments. The major environmental differences and their consequences are detailed 
below: 

MOTION: A marine reactor is constantly subjected to the motion of the vessel it is 
installed upon. It is not uncommon for a vessel to experience roll motions up to 50 degrees 
from vertical at periods of 10 to 30 seconds. Pitch and heave accelerations can also become 
significant depending on the location in the ship. On the other hand, except during earth- 
quakes, a commercial reactor can expect to have a stable platform for its entire operating life. 

This motion has several design consequences. Foundations and structural members 
must be heavier and stronger to prevent lateral movement as well as provide vertical support. 
Since ocean forces are cyclic more attention must be paid to structural fatigue in those sup- 
port members. Rotating machinery such as turbines, generators, etc. have to be oriented so 
that lateral bearing forces are reduced. This usually means that the large rotating machinery 
must be oriented fore and aft to reduce the effects of ship roll. Finally, components which 
require gravity for operation, such as control rods dropping by gravity in a scram, must take 
the motions of the ship into account in the design. 

SHOCK: Because a naval vessel has to go in harms way, shock is an important design 
consideration. The shock forces and motions produced by a weapons explosion or a collision 
are fundamentally different from earthquake loads. Shock forces have a higher frequency 
and magnitude than earthquake forces but are of a shorter duration. 

Shock has to be considered in every level of the design process. At the component 
level, individual components (electronics, turbines, foundations, etc.) have to be shock hard- 
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ened, while at the system level, the effects of shock induced motions and forces have to be 
considered in the arrangement of components, subsystems and the interfaces between those 
systems. It does not do any good to have shock hardened components if they bump into 
other components or tear out their connectors during a shock. The design consequences of 
shock are increased clearance between components, more and stronger support structures, 
and flexible connections. A full shock analysis is beyond the scope of this paper, however; 
qualitative shock considerations will be discussed where appropriate. 

Power Ramp. A commercial power reactor lives in a very stable power environment. 
It is taken up to full power and, baring casualties or emergencies, it stays at a constant power 
level all the time. When power changes are required, they are done gradually and slowly. 
This is not the operating climate of a naval power plant. Even discounting emergencies, nor- 
mal operations and maneuvering requires rapid and frequent power changes. 

This necessity to change power level rapidly and safely impacts the design in several 
ways. The reactor control system has to be more extensive than a civilian design in order to 
control power peaking and allow rapid power level changes. This means more control rods 
and/or a large negative temperature coefficient. On the non-nuclear side of the plant, tur- 
bines and other equipment must be able to respond rapidly to power changes, meaning again, 
more and faster control systems (higher capacity throttle valves, or a larger inventory control 
system.) An energy storage system (such as a large steam drum) can be used to make power 
increase transients smoother at the reactor. The storage system supplies power to the ship 
until the reactor can catch up. Likewise, for a rapid decrease in power, energy has to be 
dumped until the reactor can be powered down. 

Corrosion. Corrosion is a major concern in all power plant designs, commercial or 
marine. A marine reactor operates in the very corrosive environment of the sea. The power 
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plant itself is isolated from the sea and salt spray by the containment system so the corrosive 
effects of the surrounding environment on the design is minimal. Where the environment 
makes a difference is in the interface systems (heat exchangers and ventilation for example.) 
The corrosive environment affects the design in material selection for interface systems, and 
could require additional heat exchangers. 

1.2MGR-GT. 

The marine power plant design is based on the design of the Modular Gas Reactor - 
Gas Turbine (MGR-GT) as proposed by Lidsky and Staudt of the MIT Reactor Innovation 
Program. [8] The MGR-GT is a system intended for commercial electricity production. The 
focus of this paper is to take the commercial MGR-GT, size it to meet shipboard require- 
ments, and make any other changes necessary to adapt it to the harsh environment of a ship at 
sea. This section gives a brief description of the MGR-GT as a system. A more detailed 
description of individual parts will be in the chapters on that particular component. In the 
chapters that follow the components used in the MGR-GT will be used as the baseline for the 
marine version. 

The following is a list of the design objectives of the MGR-GT design. [8] 

• System placed in a below-grade sUo. 

• Flowpath permitting sweeping the reactor vessel with cool helium. 

• Minimum machinery-module vessel size, and no larger than the reactor vessel. 

• Flowpath maintaining low pressure difference across inner tube of concentric ducts. 

• Flowpath that minimizes system pressure losses. 

• Vertical turbomachine configuration, allowing easy access for machine removal from 
silo. 



Introduction and Background 



17 



• Net turbomachine thrust in a direction opposite machine weight to minimize thrust- 
bearing load. 

• All bearings in a low-temperature environment. 

• Generator submerged in high-pressure, low-temperature helium. 

• Plate-fin recuperator fins loaded in tension only. 

• Easy access to heat exchangers, especially the precooler, for inspection and mainte- 
nance. 

Figure 1-3 shows the layout of the MGR-GT in a below grade sUo. The plant is in two 
steel pressure vessels connected by concentric cross-flow ducts. All power system compo- 
nents except electrical power conversion and distribution equipment are contained in the 
silos. Shown in the drawing are the core, pressure vessels, ducting, generators, turbine, 
compressor, regenerators, precooler, and the inventory control vessels. 

The reactor heat source is a 200 MW pebble bed design using helium as a coolant. 
Table 1-1 lists the reactor characteristics: 

The MGR-GT is passively safe. It is designed so that the combination of a large nega- 
tive temperature coefficient and the TRISO ceramic fuel enables the core to survive a total 
loss of coolant accident with failure to scram without damage to the core or release of fission 
products. The temperature coefficient rapidly shuts down the reactor while the small size 
allows generated heat to dissipate into the surrounding environment. [8] 

The rest of the plant is a closed cycle gas turbine and associated equipment for control, 
start up, and electric power generation. The following gives a thumbnail sketch of the 
MGR-GT showing the major features 
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Table 1-1. Characteristics of the MGR-GT Reactor Heat Source. 



Core Type 


Pebble Bed cylindrical, graphite reflected on all 
sides 


Power Level 


200 MW(th) 


Active Core Length 


1 0 meters 


Core Diameter 


3 meters 


Pressure Ves.sel Length 


27 meters 


Pressure Vessel Diameter 


6 meters 


Fuel Type 


TRISO Coated Fuel Particles formed into 6 cm 
diameter spheres. 


Enrichment 


1 % 


Coolant 


Helium 


Reactor Exit Temperature 


850”C 


Moderator 


Graphite 


Refueling 


Continuous 



The compressor, turbine, and generator are on a single shaft supported by active mag- 
netic bearings. There are also auxiliary mechanical bearings used for stan-ups and 
emergency landings. 

A highly effective plate-fin regenerator is used to increase cycle efficiency. 

A shell and tube precooler using fresh water on the cold side. The water will either be 
cooled in a cooling tower or by another water source such as a river or pool. 

The generator operates at a high rotational speed to reduce size and weight and is 
cooled by high pressure helium. 

A load commutated inverter system is used for electrical frequency conversion. 
Inventory control is used for normal power level changes with turbine bypass to the 
precooler available for emergencies and shutdown cooling. 
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• System pressure (compressor outlet) approximately 8 MPa. 

Again, for a more detailed description of the individual systems see the chapter asso- 
ciated with the component. 
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Figure 1-3. 



The MGR-GT in a below-grade sUo [8] 
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1.3 Platforms. 

This is to be a modular design that should find use in almost any naval vessel with 
little or no modification to the basic characteristics of the ship. This will set the desired 
power level and will determine the optimum size for the system. Table 1-2 lists various 
naval platforms and their requirements. 



Table 1-2. Possible Platforms 


Platform 


SHP 


Shafts 


SHP/Shaft 


Power Level* 


q, = .3 


Tj, = .5 


Aircraft Carrier (CV) 


260000 


4 


65000 


680 


408 


Battleship (BB) 


212000 


4 


53000 


555 


333 


Cruiser (CG) 


80000 


2 


40000 


209 


125 


Destroyer (DDG) 


40000 


2 


40000 


209 


125 


Submarine (SS) 


40000 


1 


30000 


79 


48 


Frigate (FFG) 


40000 


1 


30000 


105 


63 


Auxiliary 


20000 


1 


20000 


52 


31 



Power level is the required plant power level (in MW) of a thermal power plant at the ther- 
mal efficiencies shown. Note: This total does not include electrical power generation or 
power needed for aircraft operations. 
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In Table 1-2 SHP is shaft horsepower. Shafts is the number of propulsion shafts, and 
power level is the thermal power level required to produce the required SHP if the ship were 
nuclear powered. Of the above platforms only the carriers, cruisers, and submarines have 
nuclear powered variants. The power level is calculated from the SHP data from the 
following formula: 

, SHP(7.457x 10^) 

Power level (MWth) = 

E Mechanical efficiency, accounts for the losses in the power transmission. It is 
also the ratio of power provided at the shaft to the power at the engine. T|„, = .95 
is used for the above analysis. 

T), = Thermal efficiency. Ratio of thermal power supplied by reactor to net work out. 
For PWR’s r\p = .30 

Based on the data in Table 1-2 a module size of 40000 SHP will be chosen. Using this 
size one module would be installed on frigates and submarines. Cmisers and destroyers will 
have two modules (one per shaft), whUe carriers and battleships would require 8 modules 
(two per shaft). The above result of one propulsion module per shaft for the smaller ships 
and two modules per shaft for the larger ships is the same as conventional designs. This is a 
convenient result since it means that if the reactor/engine room combination can be fit into 
the existing engine room volume of the smallest of the above ship types then the module 
design would fit into any ship. It would then simply be a matter of how many. In general, 
frigates and submarines have the smallest engine rooms for the installed shaft horsepower 
therefore the module design will be based on those ship types. 



Requirements and Design Philosophy. 



23 



Chapter 2 Requirements and Design Philosophy. 

This chapter covers the design boundaries and the priorities used to optimize the plant 
and make decisions in the face of conflicting requirements. The first section details the 
design philosophy and priorities while the second section details the bounds of the design. 

By bounds I mean the performance requirements and engineering details that must be met in 
order for the power system to be considered a successful design. 

2.1 Design Philosophy. 

All large system designs are study in compromise. It would be nice if all one had to do 
to design the "best" system would be to pick the best components and put them together. 

This is not the case however, and a good system design often requires a loss of capability in 
one area in order to improve the whole. Priorities must be set in order to make intelligent 
decisions about design choices. The following lists the design priorities in their order of 
importance. 



1. 


Safety. 


2. 


Size. 


3. 


Weight. 


4. 


System Efficiency. 


5. 


Technical Risk. 


6. 


Cost. 


7. 


Modularity. 



In the design process the above priorities allow trade-offs to be made in a consistent 



and systematic manner. 



Requirements and Design Philosophy. 



24 



2.2 Design Requirements. 

Based on the requirements of the ocean environment, principles of naval architecture, 
and existing ship designs, the marine variant of the MGR power plant will be subject to the 
following design constraints: 

Design Standards. U.S. Navy design standards will be adhered to for all components 
and systems. These specifications are spelled out in General Specifications for Ships (GEN- 
SPECS). GENSPECS includes Navy design standards for shock, roll, pitch, stability, and 
material quality among other things. 

Electric Drive. Because of the difficulty in producing effective rotating shaft seals for 
high pressure helium systems, the system must be sealed. To accomplish this, the only prod- 
uct of the power system will be electricity. The ship use electric propulsion and auxiliary 
systems. 

Risk. Where possible, "off the shelf components will be used to reduce technical risk. 
Previous design work on high-temperature gas turbine reactors was never funded, partly 
because the design required significant research and development work on either components 
or materials in order to be viable. Using known technology will increase the chance that pro- 
gram development and deployment could be funded at a later date. 

Size. The reactor and all associated support equipment must fit into the volume avail- 
able in existing ships. In general this consists of the reactor compartment and engine rooms 
of existing PWR powered ships. Based on the machinery box sizes of the platforms in 
section 1.3 the volume envelope for a submarine will be a cylinder 10 m in diameter by 20 m 
long. A frigate module will be used to set the size of the surface ship module. A reasonable 
average for the machinery box of a frigate size ship would be 10 m high by 10 m wide by 
20 m long. 



Requirements and Design Philosophy. 



25 



Safety. The goal is passive safety, including safety during a core flooding accident. If 
passive safety should be unobtainable for all accidents then the plant should be passively safe 
for the majority of accidents, and additional safety systems as needed are to be included in 
the design to ensure safety in the other accidents. 

Modularity. This is to be a modular design that will allow it to be included in almost 
any new construction naval vessel. Ideally the entire power generating plant (unfueled) 
would be able to be built off-site and transported to the shipbuilding site and installed within 
the vessel. The only difference between different ships would be the number of units 
installed. If the above is not feasible the module could be buUt into the ship using common 
components. 

Motion. Table 2-1 gives the GENSPECS equipment motion qualifications. At the 
ship angles listed power plant equipment must remain fully operational. The requirements 
are different for surface ships and submarines.’ 

Table 2-1. Ship angles at which equipment must remain operational. All angles in 
degrees. [9] 



Condition Submarines Surface Ships 

Submerged Surfaced 



Trim 


30 


7 


5 


List 


15 


15 


15 


Pitch 


10 


10 


10 


Roll 


60 


30 


45 



The requirements for aircraft carriers are less stringent than either submarines or other sur- 
face vessels, therefore plants which meet the above requirements will automatically be quali- 
fied for CVs. 
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In addition to the above, the plant must be able to survive a roll of 90 degrees, (ie. the 
equipment should be mounted so that it does not fall off its foundations if turned on its side 
momentarUy. 

Containment. The containment system must fully enclose all primary plant compo- 
nents. The main design goal is for the reactor compartment/engine room to be able to with- 
stand the pressure of a primary loop rupture without venting to atmosphere or other ship’s 
compartments. 

Refueling. The degree of difficulty of refueling will detennine the refueling cycle. It 
will be assumed that the core will not be continuously refueled for reasons listed elsewhere. 
If refueling is "difficult", i.e. requiring disassembly or cutting, the refueling cycle should be 
as long as possible, ten to twelve years with a minimum acceptable refueling cycle of four 
years* . This would allow refueling to coincide witli overhaul cycles. On the other hand if 
refueling is "easy", the refueling cycle can be reduced to one or two years. 

Power Requirements. The power plant must be able to provide all power require- 
ments of the ship. The power requirements of all U.S. naval vessels are as follows: 1) Pro- 
pulsion; 2) Ships service electrical (60 Hz 3-phase); 3) 400 Hz ships service electrical. 

Enrichment. There are no restrictions on enrichment except those dictated by other 
design constraints such as passive safety. 

Redundancy. A high level of reliabUity is required for all vital components. Where 
possible, redundant components will be used to ensure adequate reliability. If redundancy is 
not possible (such as the reactor itself) the component must be made as reliable as possible 



* The minimum refueling interval is set by the four year minimum allowable time period 
between overhauls as specified by GENSPECS. Current maintenance philosophy is to per- 
form major overhauls only when major equipment and system upgrades and conversions are 
needed. 
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and a backup system should be provided. Along with multiple redundant components the 
plant should be flexible enough so that the various components can be cross-connected to 
provide power in several configurations. 

Working Fluid. Helium will be the working fluid of choice. In large commercial 
high-temperarure reactors helium has proven to be tlie best working fluid. There are other 
possible working fluids, such as the other noble gasses, and their use will be discussed. 

Materials. Material selection will be the same as the MGR-GT, unless investigation 
shows that a panicular material is incompatible with the marine environment. 

Fuel. TRISO ceramic fuel will be used. The final form (prismatic, pebble bed, particle 
bed, etc.), and enrichment is discussed later in the paper. This design decision effectively 
rules out using a fast spectrum core. Although a fast core has many attractive features (com- 
pact, high power density, long life, no reactivity increase with water ingress) the focus of this 
work will be the MGR-GT core. (See Chapter 4.) 

Temperatures & Pressures. Temperatures and pressures at the controlling points in 
the cycle will be the same as in the MGR-GT. These are: 



Turbine Inlet Temperature 
Compressor Inlet Temperature 
Compressor Discharge Pressure 



850'C 
30‘C 
8 MPa 
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Chapter 3 Brayton Cycle Analysis. 

This chapter details the methods and assumptions used to analyze the closed Brayton 
cycle portion of the power plant. The object of this analysis is to estimate the cycle parame- 
ters including temperature, pressure, mass flow rate, and overall cycle efficiency. Symbols 
used are listed in the nomenclature section at the end of this chapter. 

3.1 Method of analysis. 

3.1.1 Cycle Selection. 

There are many possible machinery arrangements possible for a closed Brayton cycle 
power plant. These range from simple, single shaft systems, to intercooled, regenerated, 
multiple shaft systems with several turbine-compressor pairs and separate power turbines. In 
general, all of these systems can be analyzed to a good level of accuracy by only three differ- 
ent variations of the simple cycle; regenerated, intercooled, or reheated. Each variation adds 
to the efficiency of the basic cycle in different ways. 

The cycle chosen for this plant design is a regenerated cycle. This cycle combines 
good efficiency with the least complexity. This is shown graphically in Figure 3-1. The 
other cycles were rejected for various reasons. The reactor flow path for a cycle with reheat 
is complicated and difficult to achieve; therefore, it is not considered a viable alternative. 

The basic cycle is the simplest, requiring, the least equipment and cost. It is used extensively 
in aircraft and in open-cycle marine gas turbines, however; at the temperatures and pressures 
of this design the simple cycle has no significant efficiency advantage over current reactor 
designs. Intercooling alone does not increase effectiveness significantly over the basic cycle. 
In fact, at low pressure ratios, intercooling reduces cycle efficiency. The combination of 
intercooling and regeneration (ICR Cycle) does produce better efficiency than regeneration 



alone. 
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Figure 3-1. Effect of imercooling and regeneration on the theimal efficiency of closed 
Brayton cycles as a function of pressure ratio. [22] 

In a sensitivity analysis perfoimed by Staudt, he found that the ICR cycle was about 
3% more efficient than a non-intercooled cycle at the pressure ratio for maximum efficiency 
(which for a regenerative cycle is about 2) and the pressures and temperatures of interest. 
This 3% increase in efficiency is bought at the cost of at least one more turbomachine and 
heal exchanger along with more ducting for each power conversion loop. As in the 
MGR-GT this modest improvement is not considered worth the additional complexity. 
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3.1.2 Cycle analysis. 

Figures 3-2a and b show a recuperated Brayton cycle and defines the cycle points for 
this analysis. Each point in the cycle is defined uniquely by the pressure, temperature and 
mass flow rate, m, at that point’. The cycle analysis is performed by first defining known 
conditions at several points, along with the net power required and the performance efficien- 
cies of the various components. This becomes an iterative process since the performance of 
components is often a function of the mass flow rate, temperature, and pressure. 

Table 3-1 lists the set cycle parameters. They were chosen so that the cycle would 
resemble the MGR-GT as much as possible. The differences between the cycles are due to 
the lower power requirement of the marine version, plus the marine power plant has two 
power conversion loops per reactor while the MGR-GT has only one. 

Table 3-1. Initial Brayton cycle parameters and set points. 



Net Power Output per loop 


18.6 MW 


Compressor Discharge pressure (P2) 


8.2 MPa 


Compressor Inlet temperature (PI) 


30 'C 


Turbine Inlet Temperature (T4) 


850 ’C 


Compressor pressure ratio (r) 


2.05 


Recuperator effectiveness (£) 


0.95 


Turbine polytropic efficiency (T|p,) 


0.90 


Compressor polytropic efficiency (Tjp,,) 


0.90 


Generator efficiency (1)^) 


0.98 


Motor efficiency (Ti„) 


0.95 



Since helium can be treated as an ideal gas with constant specific heat, Cp, enthalpy (h) is 
directly proportional to absolute temperature (T) and is given by h = Cp T. 
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The reactor heat source will supply two identical CBC power generating loops. Each 
loop will provide half of total ship power requirements. Therefore each loop must provide 
20,000 HP for propulsion and 2500 MW of ships service electrical load. Net Power (W„ei) is 
the power at the turbine shaft for each loop. Converting to watts, and accounting for the effi- 
ciency of the power transmission train gives W„^,. 

The pressures and temperatures in Table 3-1 were based on the MGR-GT. The pres- 
sure ratio (r) and regenerator effectiveness were picked from Figure 3-3 which shows the 
pressure ratio for maximum efficiency based on relative pressure drop and recuperator 
effectiveness. The efficiencies, and T|^^, are somewhat lower than the efficiencies of well 
designed commercial turbomachines. The values chosen were conservative estimates of tur- 
bine and compressor efficiencies. [3] 

Once the above set points are chosen, the cycle is analyzed point by point until the 
temperature and pressure at each point is known. Each cycle point can be related to other 
points in the system by a series of well known thermodynamic principles and equations. The 
equations that follow are based on the assumption that helium behaves as an ideal gas and all 
losses in the system are included in the component efficiencies and pressure drops. 

Point 1-2. Adiabatic compression by the compressor. Points 1 and 2 are related by the 
work required by the compressor to produce the proper pressure ratio r. 
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Figure 3-2. A recuperated Brayton cycle. 

(a) Cycle diagram 

(b) H-S diagram 
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3-1. 



P 




3-2. 



T, 

W^, = C(T,-T,) 3.; 

C=mC^ 

Point 2-3. Heat addition in the Recuperator. Energy from the turbine exhaust is used 



to preheat the compressor discharge helium. The amount of heat transferred depends on the 
effectiveness of the heat exchanger used. The relative pressure drop between 2 and 3 (6P^r)is 
also determined by the regenerator characteristics. Note that this particular cycle is highly 
regenerated. There is more heat added to the working fluid in the regenerator than in the 
reactor. 



Point 3-4. Heat addition in the Reactor core. All energy input to the system occurs at 
tltis point. There is also a pressure drop which is given by equation 3-7 below. [8] This 



Reynolds number (Re) and other directly measurable core characteristics (such as height, and 
pebble diameter) This pressure drop and the drop across the regenerators are one of the main 



T5-T2 T5-T2 



3-4. 



P, = P,(l-5/>J 



3-5. 



equation uses a correlation for the core friction coefficient ('P) which is dependent on the 
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Figure 3-3. Brayton cycle efficiency as a function of pressure ratio, recuperator effec- 
tiveness and pressure loss. [8] 
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reasons to operate at high pressure. Not only does the heat transfer characteristics of helium 
improve at higher pressure but the relative pressure drop across components also drops as 
pressure is raised. 

a = C(T,-T3) 3.6. 



AP = ^ ( 1 3 3 .3 + 3 .663Re"’®') 

2d* 8' 



(l-<t>)VHG 

<1>^ pdl 



3-7. 



Re = 






f ‘ 1 






1 J 



3-8. 



Point 4-5. Adiabatic Expansion through the turbine. Initially the presence of two 
turbines will be ignored to simplify the calculations. Cycle point 4a will be estimated later 
after points 4 and 5 are found. The above produces a very good estimate as long as the 
pressure drop between the exit of one turbine and the entrance of the other is negligible. The 
power produced by the turbines is exactly the energy (per second) removed from the helium 
between points 4 and 5. The work produced also has to equal the sum of the compressor 
work and the net power. The last requirement provides the iteration condition for the system 
of equations. 

W, = C(T,-T,) = W^ + W^ 3 - 9 . 



T4 

T5 



(T-l) 

/p \“TV 






3-10. 



Mips . 
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Point 5-6. Heat transfer in the regenerator. Tlie temperature drop between points 5 
and 6 will be the same as the temperature rise between points 2 and 3. The relative pressure 
drop across the hot side i^Prh) will not be tlie same as the cold side since the regenerator may 
have a different surfaces on each side and the pressure is lower. (See sect. 5.1 and App. B.) 

P5 = P,(1-5/>J 3-11. 

Point 6-1. Heat rejection to the environment. At this stage the helium is brought back 
to the initial temperature and pressure in the precooler, rejecting its latent heat in the process. 
The relative pressure reduction across the precooler (5P^^)is generally very small and is a 
function of the heat exchanger chosen. 



a = C(Te-T,) 


3-12. 


P, = P,(l-h5P^) 


3-13. 



To perform the analysis the above equations were programmed into a spreadsheet. 

Then initial guesses for mass flow rate, and heat exchanger pressure drops were entered. 
Because the system of equations is overdetermined, (more equations than unknowns) the 
spreadsheet calculated two different values for Tj. The mass flow rate was then adjusted up 
or down until the calculated values for Tj matched. When the two values match the system is 
consistent. Because the new calculated mass flow rate may not be the same as the mass flow 
rate used in the heat exchanger analysis program COMPHX.BAS (App. B), the heat 
exchanger analysis program is run again using the new cycle values and mass flow rate. This 
process is repeated until a stable solution is achieved. This usually took no more than 4 or 5 
iterations. Once total turbine and compressor power are determined, the conditions at the 
power turbine inlet (point 4a) can be estimated by finding the helium conditions at the point 
where turbine work would equal compressor work plus ships service generator load. 



Brayton Cycle Analysis. 



37 



3.2 Results. 

Table 3-2 lists the results of the above calculations. It describes the expected helium 
conditions at the various points in the marine MGR-GT at full power with both power con- 
version loops on line. The effectiveness and pressure drops across the regenerator and pre- 
cooler were the values calculated using the heat exchanger analysis programs described in 
Appendix B. and the heat exchangers as described in section 5.1. 

Table 3-2. Conditions at the Marine MGR-GT cycle locations. 



Point 


Location 


Temp rC) 


Press. (MPa') 


1 


Compressor Inlet 


30.0 


4.00 


2 


Compressor Outlet 


143.9 


8.20 


3 


Reactor Inlet 


582.5 


8.13 


4 


Turbine Inlet 


850.0 


8.08 


4a 


Power Turbine Inlet 


717.0 


5.10 


5 


Turbine Outlet 


605.6 


4.09 


6 


Precooler Inlet 


167.0 


4.01 



Helium mass flow rate = 27.4 kg/s 



Turbine Power 
HP turbine 
Power turbine 
Total Power 



18.2 MW per side 

16.6 

34.8 



Compressor Power 



16.2 MW per side 



Cycle thermal efficiency 



48.8% 



Reactor Power Level 



38.1 MW per side 

76.2 MW total 
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3.3 Nomenclature. 

Symbol Meaning Units 

dh Hydraulic diameter m 

d^ Fuel pebble diameter m 

H Core Height m 

tit Mass flow rate kg/s 

Cp Specific heat J/kg*K 

G Mass velocity based on flow area kg/m^s 

G Average mass velocity based on frontal area kg/m^s 

P Pressure Pa 

T Temperature *K 

Q Heat flow W 

(r-reactor,c-cooler) 

W Power W 

(t-turbine,c-compressor) 

AP Pressure drop across pebble bed core Pa 

p Helium density kg/m^ 

ji. Dynamic viscosity Pa s 

Re Reynolds number 

4> Core-void volume fraction 

'F Core-friction coefficient 

e Heat exchanger effectiveness 

5P, Relative pressure drop across a heat exchanger where in place of x: 

rc = Regenerator-cold side 
rh = Regenerator-hot side 
pc = Precooler 
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Chapter 4 Reactor Design. 

This chapter details the reactor heat source design. Objectives are discussed first then 
fuel properties and choices followed by core structure and vessel design. Finally the final 
chosen design is summarize at the end of the chapter. 

4.1 Design Objectives and Considerations. 

The reactor heat source design follows closely the MGR-GT reactor design. Differ- 
ences between the two reactors are made so the marine design follows the constraints and 
requirements listed in Chapter 2. As will be seen later, this decision prevented a detailed 
study of other possible heat sources. In particular a gas-cooled fast reactor was not consid- 
ered at this time. Although a fast spectmm reactor has the potential to address many of the 
problems that will be discussed later, the focus of this investigation was to use the MGR-GT 
design as much as possible. The fast core wUl be the subject of future research. The 
MGR-GT reactor is a low-enrichment, graphite-moderated, thermal reactor. Modification of 
this design into a marine variant results in a large core." Attempts to make the core more 
compact, while still keeping within the design constraints from Chapter 2, proved to be very 
difficult. This is discussed in greater detail later in this chapter. 

There were three major core design considerations: 1) maintain full passive safety; 2) 
minimize core size; and 3) use a proven fuel design, to minimize technical risk. 

4.2 Fuel. 

For a thermal helium-cooled reactor design, at the temperatures of interest, there does 
not seem to be a better choice than TRISO based fuel. TRJSO fuel consists of micro-grains 
of fissile material surrounded by layers of refractory ceramic material. Figure 4-1 shows a 



Large in comparison to naval PWR reactor cores. 
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diagram of a TRISO fuel particle and how the panicles are incorporated into a graphite 
matrix to form a fuel pebble. Since the panicle-graphite matrix is an extmdable mixture it 
can be formed into almost any convenient shape, the most common being spheres, and hex- 
agonal prisms. The prismatic fuel is formed with channels for coolant flow and control rods. 
The TRISO fuel can also be fomied into a porous block with a very high effective surface 
area. This is the basis for a form of particle bed fuel used in the NERVA* reactor design. 

There is a large knowledge base associated with TRISO based fuel. In this country the 
Fort St. Vrain power plant uses prismatic fuel elements while in Germany the AVR., a small, 
15 MWe, pebble bed plant, in operation since 1967, has operated with helium outlet tempera- 
tures as high as 950*C. This knowledge base makes TRISO fuel a low risk choice. 



Pyr«ljrUc Carbon 
ailleoD CarbM# 
Pyroljrtio Carbon 
Carbon ^ffor 




17090 ru€l rarUo*M 



Oropbita •• bram 



Figure 4-1. Fuel Pebble sectioned to show TRISO fuel particles embedded in a graphite 
matrix. A cross section of a fuel particle is shown on the left. 



’The NERVA (Nuclear Engine for Rocket Vehicle Application) reactor is a high power den- 
sity reactor for space applications. 
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4.2.1 Fission Product Containment. 

With the TRISO coated fuel particles the two major parameters affecting fission prod- 
uct release is temperature and fuel quality. It has been determined that if fuel temperature 
remains below 1600*C there will not be significant fuel damage resulting in fission product 
release. This is shown in figure 4-2. Izenson [3] calculated that for a 200 MW pebble bed 
reactor with an outer radius of 3 meters, in an underground containment, peak fuel tempera- 
tures in the hottest region of the core did not exceed 1600'C during a depressurized loss of 
coolant accident (LOCA). This temperature is dependent on the heat transfer characteristics 
of the reactor structure and its ability to transfer heat to the ground through the containment. 
I confirmed this result with the marine MGR-GT design. (See section 4.4) 

The other significant source of fission product release is defective fuel panicles and 
entrained heavy metal atoms in the carbon layers outside the silicon diffusion layer. Given 
that the reactor does not operate at fuel temperatures above 1600*C, the defects are expected 
to be the major source of fission product release. [3] Any increase in fuel quality wUl 
decrease directly the circulating activity. Table 4-1 lists the circulating activity in current 
plants. 

Table 4-1. Activity Concentration in HTGRs with TRISO fuel. Ci/m^ [21] 



Isotope or Class 


AYR 


FSV 


Noble gases 


1.4x10'^ 


2.3x10' 


Tritium 


3.0x10-' 


5.0x10 


Sr’® 


— 


5.0x10' 


AgllOm 


2.8x10-'® 


— 


1.3. 


1.3x10-'® 


4.5x10 


Cs’^" 


8.0x10-" 


3.8x10 


Cs’” 


3.0x10" 


9.8x10 
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The values for the AVR are probably more applicable since it uses higher quality fuel 
and has not had the water ingress problems of Fort St. Vrain. It should be noted that even the 
higher levels found at Fort St. Vrain are orders of magnitude less than water cooled plants. 

Another possible avenue for fission product release is from particle damage from 
shock. Warships are subject to shock loading on occasion and it is imponant that the fuel 
remain intact when the ship is subject to shock loading (such as an explosion or collision.) 
Although the fuel compact may be damaged or broken by shock (see the next section) the 
real containment structure is the silicon carbide layer of the TRISO microsphere. This tiny 
layer is extremely tough and it is unlikely that any shock loading that the ship can survive 
would break microspheres in sufficient numbers to constitute a major release. Even though 
fission products are retained, if the fuel geometry were shattered the small panicles (dust and 
small fragments) might be a very severe contamination source. 
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Figure 4-2. Coated particle integrity against temperature. [19] 

4.2.2 Mechanical Properties. 

The graphite matrix fuel has very good mechanical properties for a reactor material. 

As listed in Table 4-2 it is dense, strong, chemically stable, and has a higher thermal conduc- 
tivity than stainless steel. Although it is not strong enough to be used as structural material, 
it has good shock resistance. Pebbles manufactured for AVR and other pebble bed reactors 
must withstand a drop of 30 feet onto a steel floor without damage. 
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Table 4-2. Typical Properties of Graphite -Matrix Fuel Compacts. [14] 



Graphite density 1.90 - 1.95 g/cm^ 

Crushing Strength* 8000 psi 

Thermal Expansion (to lOOO'C) 48.75xlO'^/*C 

Thermal Conductivity (radial) at 2000'C 0.3 W/(cm^)(’C/cm) 

Electrical resistivity 3.5xl0'^i2-cm 

Permeability (He (2) 1 atm) 2xlO'^cmVsec 

Pore Structure <1% of porosity due to pores > 1 p in diameter 

Radiation Stability <0.1% contraction after 6xl0” fissions/cm’ at 



1100- 1500*C 



4.3 Core Type. 

There are several types of core designs to consider. These are: 1) Pebble Bed; 2) Pris- 
matic; 3) Particle Bed. All of the above use the proven TRISO coated fuel particles and have 
advantages and disadvantages which will be discussed briefly below. 

4.3.1 Pebble Bed. 

The main advantage of the pebble bed core is refueling. The plant can be designed to 
be continuously refueled while on-line. With on line refueling there is a very low reactivity 
swing over the life of the plant. The fuel reaches a steady state level of bumup quickly after 
initial start-up. Although continuous, on-line refueling in a vessel at sea is not a good 
idea,**the capability exists for a partial core exchange at sea or alongside a tender to reduce 
reactivity swing over lifetime. Even when on-line or partial refueling is not being used, a 



'Values for crushing strength and thermal expansion are average values. Both properties are 
dependent on grain orientation, which is not well defined in a spherical compact. 

** The fuel handling and storage equipment would add too much, weight, volume, and com- 
plexity to the plant design to be considered worthwhile. 
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pebble bed core can be refueled in batch without disassembly of the reactor. The refueling 
evolution could be as simple as opening a series of valves and letting the spent core drain out 
of the bottom of the reactor. The new core could be loaded just as easily. In reality, this 
would have to be a tightly controlled evolution since the spent core would be extremely 
radioactive. With proper procedures and equipment I believe that the pebble bed core could 
be loaded and unloaded without major disassembly of either the reactor or the ship and, in 
the case of submarines, the large pressure hull cut required to remove the core. Without 
major disassembly the evolution should be much quicker and easier than a conventional 
refueling. Assuming the above is tme then the reactor need not be designed to go for long 
periods of time between refueling, thus less reactivity needs to be loaded into the core, sim- 
plifying the fuel design. Another advantage is that this type of core has been proven in oper- 
ation in the two German reactors AVR and THTR, thus reducing risk. 

The major disadvantage of this core type is a larger pressure drop across the core than 
either the prismatic or particle bed cores. However; as primary coolant pressure is increased 
the difference in pressure drop between the core types becomes less significant [8] 

Another disadvantage is that the pebbles are "loose" inside the core. Because of the 
motions of the ship loose pebbles can be tossed about inside the core, damaging themselves 
and reflector material. This motion can be corrected by filling the core to the top and ensur- 
ing the pebbles are fully packed. 

4.3.2 Prismatic. 

The prismatic core’s main advantage is its lower pressure drop. Even though raising 
system pressure decreases the difference between core types, the loss is still less than an 
equivalent pebble bed. The prismatic core also is a proven technology having been used in 
the Fort St. Vrain reactor. The fuel elements are securely locked in place, thus making it 
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more resistant to ship motion than the pebble bed. 

The major disadvantage of the prismatic core is in refueling. Because refueling 
requires at least partial disassembly of the reactor vessel, refueling would be a much more 
difficult evolution than in the pebble bed core. Because of this, the core would have to be 
designed to go at least eight, and preferably twelve years between refueling so that refueling 
would coincide with overhaul schedules. The excess reactivity that would have to be carried 
around to accomplish this would significantly affect the fuel design. The usual method 
employed to extend core life whUe keeping the core compact is to increase the enrichment 
and add burnable poisons. While there is no question that this could be done, it is unclear if 
it could be accomplished while still keeping the reactor safe from a core flooding accident. 
Further study is needed to answer this very important question. 

4.3.3 Particle Bed. 

The particle bed core is a very attractive alternative, promising high power densities 
with a low pressure drop. It shares the major disadvantage of the prismatic core in that it 
requires partial disassembly to refuel. Its other disadvantage is that it is new technology and 
not yet a proven concept. A small, high-power-density, fast-spectrum particle bed core, has 
great potential. 

4.4 Reactor Safety 

Passive safety is the major design requirement for this reactor, and it mms out to be the 
limiting factor on enrichment, and thus core size. Although the plant could be designed for 
engineered safety, as PWR’s are operated today, the safety equipment will add weight, vol- 
ume, and complexity to the final system, all of which are to be avoided if possible in a ship 
design. Passive safety also makes sense from a public safety consideration and it would be 
politically more attractive, (thus increasing the probability of it getting funded). On the down 
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side passively safe cores are larger thus requiring larger pressure vessels and more shielding. 
The question is which weighs less or is smaller, a large safe core or a smaller core with 
installed safety equipment. With the above in mind safety issues must be addressed in the 
following areas: 1) Fission product containment during normal operation and accidents; 2) 
Reactivity insertion due to water ingress; 3) Primary coolant retention during a depressuriza- 
tion accident; 4) Radiation shielding of personnel. Fission product retention during normal 
operation discussed above. Accidents, water ingress, and shielding are discussed in the 
sections that follow. 

4.4.1 Los.s of Coolant. 

For most reactor plants a catastrophic loss of coolant accident (LOCA) is the basis 
event. In the commercial MGR-GT design this accident limits the average power density to 
6 or 7 MW/m^ and the radius is limited by the thermodynamic propenies of the materials in 
the core and surroundings. 

To investigate this design point in the marine variant I developed the program 
HEAT.BAS to perform a one-dimensional transient heat conduction analysis. Tlie objective 
was to verify that the fmal core configuration would not exceed 1600’C during a LOCA. 
HEAT.BAS was written in QuickBASIC and runs on IBM style personal computers. It is 
fully described in Appendix A. along with a sample input fUe and source code. 

4.4.1. 1 Problem Assumptions. 

The following summarizes the assumptions used in modeling the reactor and initial 



conditions for the problem. 



memm 
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1. Surface of the shielding water tank remains at a constant 50°C. The shielding water 
tank is a large water filled tank which has both its own cooling system and is in direct 
thermal contact with the sea through the hull of the ship. Given the long time frame 
associated with the postulated accident, 50” is probably a conservative estimate. 

2. No credit is taken for axial heat conduction, heat loss through the ends of the core, or 
any natural circulation of helium within the core. Again this assumption will produce 
a conservative estimate. 

3. The power density, and initial temperature distribution will be the same as the 
MGR-GT since the core radius is approximately the same and the exit temperature is 
the same. The values used are given in Table 4-3. They represent the steady state 
conditions at full power at the axial position in die core where the highest temperature 
is reached after the casualty. This is approximately where z/H = .2; z is distance from 
core entrance and H is core height. [7] 

4. The volume fraction of pebbles is = 0.61. 

core 

5. The pebble diameter is 6 cm. 

6. The pebble bed emissivity is 0.8 and the emissivity of outer surface of the core barrel, 
the outer and inner surface of the pressure vessel, and the inner surface of the shield- 
ing water tank is 0.6. 

7. The accident is assumed to start at time 0 with depressurization of the reactor core. 
Any helium cooling due to expansion is ignored and the reactor power falls off 
according to the decay heat fraction (ie. the reactor trips at time 0 with the core at 
atmospheric pressure.) 
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Table 4-3. 



Initial power density and temperature distribution used in the transient acci- 
dent analysis. 



Radius (m) 


Temperature 


0.05 


582 


0.1 


580 


0.2 


579 


0.4 


572 


0.5 


565 


0 . 6 


560 


0.7 


552 


0.8 


544 


0.9 


542 


1.01 


539 


1.06 


538 


1.12 


537 


1.24 


536 


1.30 


531 


1.35 


530 


1 . 45 


425 


1.55 


320 


1.57 


310 


1 .77 


200 


1.80 


185 


1.85 


170 


1.87 


50 



Power Density (W/cc) 

8.21 
8.15 
8.10 
7 . 60 
7.30 
7.25 
7.24 
6.80 
6.67 
7.33 

7.00 
6.87 
6.90 

3.00 
0 

0 

0 

0 

0 

0 

0 

0 



9. 



All heat generation is in the active core, radiation heating in the reflector or vessels is 
not included. 

Dimensions and materials used are shown in Figure 4-3. 
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Material Legend 



Figure 4-3. One-dimensional model for the Marine MGR-GT reactor showing dimen- 
sions and materials used. 

4.4.1. 2 Results. 

HEAT.BAS was run using the above assumptions, dimensions, and materials. Figure 
4-4 shows the centerline temperature plotted against time since event initiation. After the ini- 
tial sharp temperature increase, the heat generation decays until heat rejection into the sur- 
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rounding water equals the heat addition. For the baseline case (linear heat generation rate = 
39 MW/m) peak temperature reaches 1400*C at t=39 hours. The temperature then gradually 
falls off, until at 100 hours after the event, centerline temperature drops to just under 1300°C. 




Tim© (hours) 



Figure 4-4. Core centerline temperature as a function of time after a depressurized loss 
of flow accident (baseline power level 39 MW/m). 

Notice that even with the conservative model, centerline temperatures never get close 
to the 1600’C required for fuel damage. There is a safety margin of 200 degrees at the base- 
line power level. 
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Figure 4-5 shows the radial temperature distribution for various times after the casu- 
alty. The top trace corresponds to the maximum centerline temperature 39 hours after event 
initiation. Of note here is that even though the center is at 1400*C the pressure vessel remain 
well below any temperature where damage could occur. Note that this is not the point of 
maximum vessel temperature. That would be near the top of the pressure vessel as the heat 
in the core rises. 

A simple sensitivity analysis was also performed to determine what power level the 
reactor could operate at without exceeding the 1600*C limit. This is show in Figure 4-6. The 
three cases correspond to linear power level of 39, 49, and 59 MW/m (100%, 115%, and 
125% baseline power level). The initial temperature distribution is the same for all tliree 
cases’ and the power density was increased by the same percentage at each node. Each case 
was run until centerline temperature reached a maximum. The results were then plotted and 
compared to baseline. 

As the power level increased, the reactor reaches maximum centerline temperature 
quicker and at a higher value (1560*C at 37.7 hours @ 49 MW/m, 1703*C at 36.2 hours @ 

59 MW/m). The results indicate that it would be safe to operate the reactor at 49 MW/m. 
However, since operating at 49 MW/m does not allow a reduction in reactor size because of 
criticality, the 200 degree safety margin afforded by operating at the lower power level will 
be used as a hedge against uncertainty. 



* It was assumed that the temperature distribution would not be significantly different at dif- 
ferent power levels since the reactor operates at a constant exit temperature with the mass 
flow rate setting power level. 
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4-5. 



Radial temperature distribution histoiy following a depressurized loss of 
flow accident. Temperature distribution peaks at 39 hours after event. 
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Tim© (hours) 



Figure 4-6. Effect of increasing linear power level on the core centerline temperature as 
a function of time after a depressurized loss of flow accident. The traces end 
at the time of peak temperature. The power level indicated is the average 
power produced by one meter of core length. 

4.4,2 Shock Loading. 

The shock resistance of graphite moderated cores is an open question. Nuclear graphite 
when new has good strength and shock resistance (see sect. 4.2.2.) Whether or not these 
properties persist as they age in a neutron flux remains to be seen. 

Of the three core types mentioned above the pebble bed is the most shock resistant. 
Assuming a shock load sufficient to break core graphite but not other core structures (control 

I 

rods, piping, supports, etc.) the pebble bed’s propenies will not change much. As mentioned 
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above, damage to the pebbles does not cause loss of fuel containment. The broken pebbles 
will simply act like smaller pebbles. The above does not apply, of course to the graphite 
reflector, which could be significantly damaged by shock. For that reason, in the marine 
variant, the amount of graphite reflector is kept to a minimum, serving more as a thermal 
buffer for the core barrel than a reflector. Instead, the shielding water tank also serves as a 
reflector. The axial reflectors are unfueled pebbles so they are also immune to shock. With 
this configuration even if the radial graphite reflector is damaged the pieces should not move 
much and will tend to act as another unfueled pebble. 

In a prismatic core the fuel is cooled by helium flowing through channels molded into 
the prismatic elements. Shock damage to this type of configuration could change the flow 
characteristics tremendously. Again, containment is not lost, but fragments could block flow 
passages in other elements causing local hot spots and possible fuel failure. 

The conclusion is that the marine variant MGR-GT core should be very shock tolerant. 
4,4.3 Water Ingress. 

Water ingress is a serious casualty in a graphite moderated reactor for several reasons. 
The first of which is that the superheated water vapor could react with the hot graphite caus- 
ing degradation of the reactor structure and potential release of entrained fission products. 

The second and more serious problem is the large insertion of reactivity caused by the water. 
Izenson comments in his paper that the reactivity effects of water ingress were mostly a 
function of the heavy metal content of the fuel. By limiting the heavy metal content of the 
fuel to about ten grams per pebble and about eight percent enrichment, water ingress of any 
amount did not cause a reactivity increase greater than the delayed neutron fraction. [3] This 
is shown in Figure 4-7 as the increase in reactivity (AA'^y) as a function of the amount of water 
entering the core. Reactivity insertion was calculated for both the base core and a core con- 
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taining a small (0.1% wt) Gadolinium poison. 




Mass Water (kg) 

Figure 4-7. Reactivity effect of water ingress on the base and poisoned cores of the 

MGR-GT. The upper trace is the base core and the lower trace is the poi- 
soned core. Enrichment = 7.5% wt U‘^^. The poisoned core has 0.1% wt 
Gadolinium. [3] 

Although Gadolinium lowers the reactivity insenion the effect is slight and any advan- 
tage gained in the water ingress is offset by an increase in reactor power peaking. [3] 

The Westinghouse Astronuclear Laboratory (WAJ4L) conducted a feasibility study in 
1979 of a large highly enriched gas-turbine reactor for a low specific weight power plant. 
They concluded that in that particular core the reactivity insertion due to core flooding was 
sufficient to cause the core to go super prompt critical. This is a serious consequence for a 
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reactor that is surrounded by water, and it tends to support Izenson’s results that high enrich- 
ment increases susceptibility to damage from water ingress. The WANL reactor design oper- 
ated in the epithermal region so their solution was to use a spectral shift type poison such as 
Gadolinium to reduce the reactivity. [4] Another possibility is to go to a fast core and again 
use thermal poisons to shut down the reaction when water enters. Much future work needs to 
be done in this area. 

From the above, the marine core will be limited to the same 7.5% enrichment of the 
MGR-GT. 

4.4.4 Primary Coolant Retention. 

An accident which causes the sudden depressurization of the primary system is poten- 
tially serious due to spreading of contamination, and the resulting pressurization of the reac- 
tor compartment/engine room. This is most significant for submarines because venting to the 
atmosphere is not an option. Assuming a catastrophic breach in the pressure boundary, 
turbine deblading or pipe rupture for example, with the entire helium inventory, including the 
inventory control vessels, venting into the compartment. The helium would expand into 
about 3 to 4 times their initial volume. With the primary system at 80 atmospheres, the com- 
partment will equalize at between 16 to 20 atmospheres (530 to 660 feet of water). Modem 
submarine pressure hulls should easily withstand this pressure*. Thermal shock to the 
pressure boundary should not be a great problem since in general the reactor compartment is 
surrounded by fluid which will keep structural temperatures down. 



Although submarine pressure hulls are buUt to withstand external pressure, ring-stiffened 
cylinders are stronger with internal pressure than external pressure since buckling is not a 
problem with a cylinder in tension. 
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4.5 Core Size. 

Core size is a complicated function of several factors; these are, criticality, thermal 
hydraulics, core life, control margin, and safety. The interaction of these factors will set or 
limit reactor characteristics. Given a pebble bed graphite moderated core with TTUSO fuel 
particles, this section details how the reactor was sized and the assumptions used. 

4.5.1 Criticality. 

Reactor criticality calculations were performed using the Nodal Graphite Code (NGC) 
developed by Ediz Tanker. NGC is a two-dimensional nodal code for solving group diffu- 
sion equations. [15] The code uses a quadratic nodal scheme which pemiits mesh spacings 
twice as large as those required for finite difference codes. The program is written in 
FORTRAN and runs on desk-top personal computers. NGC requires the user to divide the 
core into a convenient mesh and provide the group cross sections for each node. Once the 
above is done the program will calculate k^ff. 

The cross sections used to estimate k^ff for this design were the group cross sections 
used in tlie MHTGR by GA Technologies. The MHTGR is a five region annular, cylindrical 
core, graphite reflected on all sides, with an inner reflector. The cross sections were aver- 
aged over the MHTGR spectrum to produce the group constants. Although the spectrum for 
the MHTGR is not the same as the Marine MGR-GT it provides a starting point and it will 
give general limits for core and reflector size. 

For my NGC calculations I used two-group cross sections and two material regions, 
active core, and graphite reflector. For the core region I used the MHTGR core cross sec- 
tions. For the radial and axial reflector I used the cross sections for MHTGR region five. 
The following summarizes the results using this configuration. 
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Bare core. Bare core calculations were performed analytically and with NGC to get a 
starting point and to confirm NGC results. Core radius had to be at least 1 meter in order to 
achieve criticality (k^ff =1). Criticality was not very sensitive to core length and as long as 
the length was greater than about 2 meters the core could be effectively treated as an infinite 
cylinder. 

Reflected core. Once bare core calculations were completed reflectors of various sizes 
were added to try to bring up to 1.25. A k,ff of 1.25 was used as a rough margin to give 
the core some extra reactivity for Xenon override, and to allow control margin and bum-up. 
Table 4-4 lists the results. 

Table 4-4. Initial Core Size Estimates. 



Active Core Radius 


1.25 m 


Active Core Length 


3.00 m 


Radial Reflector Thickness 


50 cm 


Axial Reflector Thickness 


75 cm 


Approximate enrichment 


7% 


k.ff 


1.22 



The above results were more reassuring than useful. What they showed is that the 
critical dimension is radius. Length had very little effect as long as it was greater than the 
core diameter. Also the addition of a reflector did not decrease total reactor radius. The 
smaUest total radius (core & reflector) was the bare core. There was no combination of 
reflector and core that I could find which resulted in a smaller total radius than the bare core. 
The effect of the reflector was to either maintain the same k^jf with a smaller fueled region, or 
increase k^ff while keeping the diameter of the fueled region constant. They were also 
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reassuring in that even though the cross sections used were for an annular prismatic core, the 
dimensions were similar to the dimensions of the MGR-GT. Since length only seemed to 
effect overall power level and not k,ff, the above results led to the decision to use the MGR 
core dimensions as the basis of the marine design. 

4.5.2 Core length. 

Reactor length was set primarily by the core power level. From Chapter 3, the required 
reactor power level to produce 40,000 SHP is about 80 MW(th). From 4.4 the core can oper- 
ate safely at 39 MW/m of core length at the chosen diameter. Simple arithmetic results in an 
active core length of about 2 meters. However, since this is less than the core diameter (2.7 
meters) setting the reactor length at 2 meters would result in excessive neutron leakage and a 
lower k,ff. For this reason I set the reactor length at 3 meters. The increased length wUl 
allow the reactor to operate at a lower power level than it is capable of. 

The benefits of the lower power level are: 

1. A longer core life since tliere is less bum-up. 

2. Increased safety margin. 

3. More control drum worth due to the increased length and better aspect ratio. A 
long thin core is better for periphery control than a short fat one because the rods 
are longer and the radial leakage is higher in the thin cores. A short fat core 
would be better if internal control rods were used. 

4. Heat transfer is enhanced since tlie heated length is greater. 

The major drawbacks are: 

1 . The larger core is harder to arrange within the ship envelope. 

2. More shielding is required to cover the increased length 
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3. Core pressure drop is directly proportional to core length, so a long thin core has 
a higher relative pressure drop than a short fat one at the same pressure and mass 
flow rate. 

Another alternative is to increase system power level to use the extra length. At 40 
MW/m a 3 meter core is capable of being safely operated at 120 MW. This would corre- 
spond to a system power level of 60,000 SHP. Since the power conversion equipment (tur- 
bines, generators, heat exchangers, etc.) would have to be enlarged to match this is not a 
viable option for a small ship design. It may be a good choice for the carrier design since it 
would drop the total number of plants needed from eight to five or six. 

4.6 Radiation Shielding. 

Because of the low activation cross section of helium, it is not expected that the pri- 
mary system outside of the reactor itself will require shielding from manned areas. Any fis- 
sion product plate-out or radioactive dust which may accumulate on piping or in components 
will cause radiological precautions to be needed during maintenance but the component 
pressure vessel and insulation along with the reactor companment bulkheads should be suffi- 
cient for the primary system. 

Primary reactor shielding, intended mainly to stop neutrons, will come from the shield- 
ing water tank which surrounds the reactor. The shielding water tank also acts as a heat sink 
as explained above. Additional shielding is provided by fuel oil* or other tanks separating 
the reactor compartment from manned spaces. Where intervening tanks are not possible, 
such as the shielded tunnel in a submarine or above the reactor compartment in a surface ves- 
sel the shielding will be extensive. 



All nuclear powered ships have emergency power generators (diesel or gas turbine) for use 
when reactor power is unavailable. 
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Because of the large size of the reactor core the shielding will be a major component of 
the plant weight. To estimate the shield weight reference 17 will be used. In reference 17 
several shield designs were presented for a gas-cooled, intemiediate spectmm reactor, at a 
power level of 100 MW thermal. A lead-water shield is used to reduce dose rate to 0.5 
mR/hour at the edge of the shield. Water is used as the thermal neutron shield and lead is 
used for the gamma shield. The shield consists of 25 cm of lead at the pressure vessel sur- 
rounded by 180 cm of water. Using the core dimensions listed below in Table 4-5 the shield 
dimensions can be estimated. The thermal neutron shield will be a 1.8 meter thick water tank 

6.7 meters high. The tank inner radius is 1.85 meters and the outer radius is 3.65 meters. 
Surrounding the core tank is a 25 cm thick lead y-shield. The lead is positioned several centi- 
meters away from the inner wall of the shielding tank to allow fluid circulation. This results 
in approximately 230 metric tons of water and 260 metric tons of lead, a total radial shield 
weight of about 490 metric tons. This figure can be taken as a practical upper limit of shield 
weight since in a actual shield design the reactor would not be shielded equally all around. 
Heavy shielding would go around manned spaces, such as the overhead in a surface ship and 
the shielded passageway in a submarine, while other areas would get little or no shielding 
(the bottom for example). 

4.7 Reactor Design Summary 

This section summarizes the design of the reactor heat source. Figure 4-8 shows a cross 
section through the marine MGR-GT reactor core. Figure 4-9 shows two side profiles. The 
left hand side is a section through a control drum and the right hand side is through the co- 
axial helium ducts. Co-axial ducts are used to reduce insulation requirements on ducting and 
to lower the temperature of pressure boundaries. Notice that with this flow arrangement the 
reactor vessel is continuously swept with ’cool’ compressor discharge helium. 
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Figure 4-8. Section through Marine MGR-GT core showing pressure vessels and control 
drums. The reactor is surrounded by the shielding water tank and is sepa- 
rated from the tank by a 2 cm air space. The hatching in the central core 
denotes pebble, not hexagonal, fuel elements. 



Table 4-5 summarizes the marine reactor characteristics and Table 4-6. gives the 
weight summary. Weights were estimated using the geometry shown in Figures 4-8 and 4-9 
along with known density for the various materials. This weight estimate does not include 
foundations. 
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Figure 4-9. 



Side profile of the Marine MGR-GT showing gas flow path. Pressure vessel 
and Control drum are cooled by compressor discharge flow. 
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Table 4-5. Marine MGR-GT Reactor Summary. 


Pressure Vessel 




Height 


6.7m 


Outer Diameter 


3.7 m 


Thickness 


8 cm 


Active Core 




Type 


Pebble Bed 


Height 


3 m 


Diameter 


2.7 m 


Fuel 


~ 92645 Fuel Pebbles 
6 cm diameter 
loaded with 7 g Uranium 


Enrichment 


~ 7.5% 


j Control 


1 6 Reflector/Absorber 
Dmms. 



Table 4-6. Reactor Weight Summary. 


! Component 


Weight (metric tons) 


1 Core pebbles 


17.3 


Graphite reflector 


4.5 


Core Barrel 


10.5 


Control drums 


9.0 


Pressure Vessel 


44.7 


Upper Reflector Pebbles 


8.8 


Lower Reflector Pebbles 


8.8 


Total 


103.6 
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Chapter 5 Mechanical Design. 

5.1 Heat Exchanger Design. 

System heat exchangers form a crucial part of the overall design. Because of the rela- 
tively poor heat transfer characteristics of gasses, compared to water, heat exchangers of high 
effectiveness, low pressure drop, and small size are difficult to create. Small size is 
obviously necessary to reduce both weight and volume. High effectiveness and low pressure 
drop are critical factors in cycle thermal efficiency. This was shown in Chapter 3, Figure 
3-3. In order for this design to be successful reasonable sized heat exchangers with good 
effectiveness and pressure drops need to be found. 

This section details the design of the main system heat exchangers, the regenerators and 
the precoolers. 

5.1.1 Regenerator. 

The regenerator is the only gas-gas heat exchanger in the entire system. Shell and tube 
designs have traditionally been used in commercial applications because they are rugged and 
relatively easy to inspect and repair. Their main disadvantage is size. Because of the need 
for compactness, traditional shell and tube heat exchangers will not work. 

In recent years the technology necessary to design and buUd compact plate-fin heat 
exchangers for use in high-temperature, high pressure difference applications has greatly 
improved. Their greatest use has been in the large gas-turbine powered, gas pipeline com- 
pressors. These regenerators have accumulated thousands of hours of operating experience 
with not one failure. [8] Using the same technology on the marine MGR-GT should be a 
relatively simple task, especially since the marine recuperator wiU be operating in a clean 
helium environment. The clean environment means that the fouling caused by combustion 
gasses cannot occur, thus the heat exchanger passages can be made as small as necessary to 
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meet design objectives. 

The regenerators were designed using the effectiveness-NTU method as detailed in 
Kays and London’s Compact Heat Exchangers. [11] The mechanics of this method were 
incorporated into a heat exchanger analysis program COMPHX.BAS. COMPHX.BAS and 
the effectiveness-NTU method are detailed fully in Appendix B. 

The first step in using the program was to define general cycle parameters (Chapter 3) 
such as required effectiveness, maximum pressure drop, and entrance and exit temperature 
and pressure. Then several heat exchanger surfaces were selected from the surface geome- 
tries listed in Reference 1 1 These surfaces are listed for reference in Appendix C. 

The main selection criteria was the compactness ratio, a, which is defined as the ratio 

of surface area to heat exchanger core volume, the higher the a, the more compact the heat 
exchanger. For this study, every heat exchanger surface with a greater than 800 ftVft^ (2600 
mVm’) was analyzed using COMPHX and the Brayton cycle analysis spreadsheet from 
Chapter 3. Of the five surfaces used the best performance came from surfaces 46.45T and 
1/9-24.12 in Kays & London*. The characteristics of these two surfaces are shown in Figures 
D-1 and D-2. Surface 46.45T has plain fins which extend the full length of the heat 
exchanger, while in surface 1/9-24.12 the fins are short offset strips which serves to interrupt 
the flow. 

For the initial analysis at least, the same surface was used on both hot and cold sides of 
the heat exchanger. After the results of the initial analysis were done it was noticed that most 
of the total pressure drop through the heat exchange was on the hot side. The fluid velocity 
on the hot side is higher than the cold side because of the lower pressure (therefore density.) 



These surfaces are numbered 1-18 and 3-17 in Appendix C. 
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To try to mitigate the hot side pressure loss, a surface with larger flow passages was used on 
the hot side (1/8-20.06) with surface 1/9-24.12 on the cold side. While this configuration 
successfully lowered the pressure loss, it resulted in a much larger volume and weight, ren- 
dering it unacceptable. Table 5-1 summarizes the results of the analysis for the three sur- 
faces, 46.45T, 1/9-24.12, and the combination 1/9-24.12 & 1/8-20.06. Stainless steel was the 
material used with a separating plate thickness of . 1 cm. 

Table 5-1. Regenerator performance analysis and sizing results. 



0 


1/9-24.12 


46.4.57 


Combination 


Size Summary 


Core length (m) 


1.5 


1.5 


1.75 


Frontal area (m^) 


2.5 


2.5 


3.5 


Height* (m) 


2.4 


2.4 


2.6 


Core volume (m^) 


3.75 


3.75 


6.13 


Weight (Mtons) 


10.4 


14.2 


17.2 


Performance Summary 


Effectiveness {%) 


94.9 


93.2 


94.9 


NTU 


18.6 


13.7 


18.9 


Pressure drop (%) 


Hot-side 


1.81 


1.06 


.35 


Cold-side 


.71 


.42 


.93 


Total 


2.52 


1.48 


1.28 


System Performance 


Massflow (kg/s) 


27.4 


26.6 


26.47 


Thermal efficiency 


.488 


.487 


.497 



The performance of the three regenerators is very close, especially between 46.45T and 
1/9-24.12. The plain fin surface is less effective and weighs more, but has a lower pressure 
drop than the strip fin surface. As a result they have an almost identical effect on system 



Height includes cross-flow headers. 
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performance. The best system thermal efficiency is achieved by the combination heat 
exchanger. However to achieve a 95% effectiveness it had to be 7 tons heavier and twice the 
volume of 1/9-24.12. In a tight design, the marginal efficiency increase does not justify the 
weight and size. Though either of the other two would be acceptable, the lower weight of the 
strip-fin exchanger gives it the advantage. 

5.1.2 Precooler. 

Unlike the regenerator the precooler has helium gas on only one side. With water on 
the other side, corrosion and fouling must be taken into account. For this reason a simple 
cross-flow shell and tube heat exchanger will be used. The helium will be on the shell side 
with water flowing through the tubes. Shell and tube heat exchangers work well with a high 
pressure difference between the two fluids and can easily be made robust and shock resistant. 
With this arrangement the tubes can be easily finned to augment the heat exchange surfaces 
(although it will be seen later that this does not produce the most compact design) and the 
water plenum and tubes can be accessed relatively easily for maintenance. 

The use of tubes also facilitates locating and fixing leaks. A plate-fin exchanger would 
be more compact but a leak between the two sides would be difficult to locate and even more 
difficult to repair. In the regenerator a leak wUl reduce system performance by a small 
amount but it will not damage the system. A leak in the precooler, however; will either result 
in potentially contaminated helium escaping to the environment or worse, water ingress into 
the primary system. A leaky tube is easy to locate and repair (plug or replace.) 

Once the choice of heat exchanger type is sizing is done using the same effectiveness- 
NTU method employed in sizing the regenerators. The difference between the two being two 
different fluids and cross-flow instead of counter-flow. The major concerns in the design 
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are: 1) minimum size, 2) minimal pressure drop, 3) reasonable water mass flow. For the pur- 
poses of this analysis the parameters listed in Table 5-2 considered as constants in the prob- 
lem. 

Table 5-2. Input parameters for use in the Precooler design. 



Parameter 


Value 


Helium Conditions 




Massflow 


27.4 kg/s 


Temperature: Inlet 


167 “C 


Outlet 


30 *C 


Inlet Pressure 


4.01 MPa 


Pressure Drop 


.25 % 


Water Conditions 




Inlet Temperature 


20 *C 


Massflow 


250 kg/s (-4000 gpm) 


Heat exchaneer Material 


304 Stainless steel 



With the input conditions defined the next step is to select several candidate surface 
geometries from Reference 1 1 Only heat exchangers with circular tubes were considered for 
reasons given above. Four configurations were selected based again on the largest values of 
a. Two of them were bare circular tubes with different diameters and tube pitch, one con- 
tained circular tubes with circular fins and the last used circular tube with continuous fins 
(like a car radiator). Figures D-3 through D-6 show the surface geometry and performance 
characteristics. The characteristics are included as data statements in the PRECOOL.BAS 
source code. 

The program PRECOOL.BAS was then run to size the precooler using the four sur- 
faces. The program is explained in detail in Appendix B, however what it basically does is 



Mechanical Design. 



71 



pick an initial guess at the dimensions of the heat exchanger, then analyze it using the helium 
and water parameters from Table 5-2. If the results do not match requirements the dimen- 
sions are adjusted until they do. Table 5-3 lists the results. 

Table 5-3. Precooler performance and sizing results. 





CF-8.72 


S 1.50-1.00 


S 1.50-1.25 


8.0 3/8T 


Tube length (m) 


1.93 


1.92 


1.20 


1.87 


Width (m) 


.82 


.73 


1.19 


.72 


Depth (m) 


2.59 


.99 


.756 


3.38 


Frontal area (m^) 


1.58 


1.40 


1.43 


1.43 


Volume (m^) 


4.09 


1.39 


1.08 


4.84 


Number of tubes 


4215 


5289 


11901 


4603 


Weight (kg)“ 


3564 


534 


346 


3528 


Effectiveness 


.932 


.932 


.932 


.932 



Based on the results in Table 5-3, a precooler using surface S 1.50-1.25 is clearly the 
best choice. It is somewhat surprising that the bare tubes did better than the finned surfaces. 
The explanation for this is that the fmned surfaces rely on the fins for most of the heat trans- 
fer. The tubes are relatively widely spaced in comparison to the bare tube surfaces. The rel- 
atively poor thermal conductivity of stainless steel reduces the effectiveness of the fins. 
Therefore the primary heat exchange surface is the tubes themselves. The bare tube surfaces 
use smaller tubes at a fine pitch so they have more effective heat exchange area per unit vol- 
ume than the firmed surfaces. If a more conductive material were used (such as aluminum) 
the situation would be different. 



’ Weight is weight of tubes and fins. It does not include headers, support plates, or the outer 
envelope. The weight of the other items will be roughly proportional to volume. 
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5.2 Bearing Design. 

One of the major problems in closed cycle systems is keeping the circulating fluid 
clean. Contaminants in the system are constantly circulating unlike a open cycle where they 
quickly leave the circuit. Keeping the system clean is especially imponant in this reactor 
design. If lubricants from the bearings of the rotating machinery were to enter the circuit the 
results could range from contamination and fouling of flow passages in heat exchangers and 
the core to chemical reaction with core graphite and reactivity excursions. Much of the poor 
availability of the Fort St. Vrain nuclear power plant can be traced to water ingress from the 
water-lubricated bearings on its gas circulators. [8][18] 

5.2.1 Magnetic Bearings Properties 

Although fluid lubricated bearings are generally very reliable and the technology well 
established prevention of leakage requires complicated seals and extensive supporting equip- 
ment. Gas bearings have been used with great success in small high speed systems, however, 
the technology base for the large gas bearings that this reactor design requires does not exist. 
Because of the severe potential problems associated with mechanical bearing lubrication sys- 
tems (either water or oil). The only bearing system considered for this design was active 
magnetic bearings. 

Active magnetic bearings work by levitating the rotating equipment with a magnetic 
field instead of a fluid film. The technology is not new, and several companies (mostly 
French) currently market bearings and control systems large enough to handle the needs of 
the marine MGR-GT. Table 5-4. lists the major advantages of magnetic bearings in this 
application. 

Bearing Configuration. 
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Table 5-4. Advantages of Active Magnetic Bearings. [18] 



• Elimination of system contamination by lubricant ingress 

• Elimination of fire and explosion hazards 

• Elimination of costly lubrication equipment and complex sealing and gas buffering sys- 
tems 

• Lower power consumption and overall system volume than comparable fluid lubricated 
bearings 

• Reduced bearing frictional parasitic losses 

• Unlimited bearing service life (no contact or wear surfaces) 

• High rotational speed capability 

• Vibration free operation (ability to pass critical speeds) 

• Alignment and balancing simplification 

• Continuous monitoring of rotor status 

• Responsive to rapid system transients 

• Reduced maintenance 

• Simplicity of operation and control (computer system) 

• Potential for very high reliability. 



Magnetic bearings consist of two types: radial bearings, which maintain the rotor posi- 
tion in the center of the bearing; and axial bearings, which maintain the longitudinal position 
of the rotor. Radial bearings resemble electric motors in constmction. (Fig. 5-1) The outer 
stator is electrically divided into four electromagnetic quadrants. Each quadrant contains 
north and south poles which attract the rotor when current is supplied to its poles. In an axial 
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bearing, a flat solid ferromagnetic disk, fitted radially to the shaft, is used for the rotor. Elec- 
tromagnets are positioned on both sides of the disk. By varying the current flowing to each 
magnet the disk is attracted to either side thus creating a double acting thrust bearing. (Fig. 
5-2)[20] 




Figure 5-1. Radial Magnetic Bearing. [20] 




Figure 5-2. Axial Magnetic Bearing. [20] 
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Since the magnetic bearing is a non-contact bearing, auxiliary mechanical bearings 
have to be included as a backup in case of power failure, bearing overload, or other bearing 
failure. They would only be in use during shutdown, or other temporary or emergency situa- 
tions. The bearings have to be dry-lubricated to prevent the very same lubricant ingress 
problems that the magnetic bearings were being used for in the first place. Figure 5-3 shows 
how the auxiliary bearings would be employed in a magnetic bearing machine. 




Figure 5-3. Basic arrangement of auxiliary bearings in machine with magnetic bearings. 
[ 8 ] 

Bearing Control. 

Since attractive magnetic bearings are unstable an active control system is used to keep 
the bearing in position (Fig 5-4.) Integral to the bearing position sensors are used to continu- 
ously determine the exaa location of the rotor. This position information is then fed back 
into a closed-loop control system with a proportional, integral, derivative controller. The 
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control system generates an error signal which goes to power amplifiers. The resulting cur- 
rent change corrects any error in the rotor position. The damping and stiffness of the system 
is a function of parameters programmed into the controller. Unlike mechanical bearings, 
these parameters can be changed as conditions require, making the bearing very responsive, 
esp>ecially at shaft critical speeds. 




Figure 5-4, Diagram of an active magnetic bearing control system. [8] 

Vibration Reduction and Automatic Balancing. 

Another attractive feature for magnetic bearings is their ability to automatically balance 
a rotor through inertial axis control, and active vibration reduction. Figure 5-5 shows a 
representation of how automatic balancing works. The rotor effectively is allowed to rotate 
around its inertial axis instead of its geometric axis. As long as the imbalance is not too 
great, the result is very quiet and smooth operation with no forces transmitted into the sup- 
ports. Also by placing vibration sensors on the bearing the control system can create oppos- 
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ing forces to cancel the vibration. In tests on a auxiliary turbine for a forced draft blower 
vibration and sound level reduction of between 10 and 25 dB were achieved.[20] This is 



very attractive for a submarine application where noise reduction is very important. 



Rotor Alone 




Unbalance: 0.07 oz. 



Conventional Bearings 



Active Magnetic Bearing 
Inertial Axis Control 





Transmitted Imbalance 
Force Per Bearing = 400 lbs. 



Transmitted loibalance 
Force Per Bearing = 0 lbs. 



d = Rotor Diameter = 2 in. 

Rotor Weight: 22 lbs. 

Rotor Speed: 60,000 RPM 

G: Geometric Axis 

I: Inertial Axis 

e = Difference Between Axes: 0.4 mils 



Figure 5-5. Automatic balancing: a diagrammatic representation of how inertial axis 
control works. [20] 



Reliability. 

Although reliability data on installed magnetic bearing systems are still being accumu- 
lated, it being a new technology, the operating data as of 1985 on industrial applications, 
showed that in over 38,000 hours of operation there were 26 failures (12 on one machine). In 
all cases the power amplifier was at fault and in no case was there any damage other than the 
failed electrorucs. [8] Therefore the reliability of the system seems to be a function of the 
reliability of the electronics. This is fortunate, from a maintenance point of view, since it is 
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much easier to carry extra electronic modules and replace them at sea than to repair bearings. 
As experience is gained, the problems with the early systems should be corrected and mag- 
netic bearings should become very reliable. 

5.3 Turbomachine design. 

The turbine and compressor design will foUow the design method proposed by Staudt 
for the MGR-GT. The method used to design the turbines for the marine variant is to take 
the baseline machine from Reference 8 and apply some scaling relationships. For this pur- 
poses of this paper both turbines in the split shaft arrangement are assumed to operate at the 
same speed. This is not necessarily true but as will be seen later both shafts have to drive 
generators, and it is much more difficult to design an acceptable high speed generator than it 
is to design a high speed turbine. 

5.3.1 Baseline Turbomachine. [8] 

The MGR-GT turbomachine combination is used as the baseline for the scaling rela- 
tionships of the next section. Table 5-5 lists the characteristics. 

Table 5-5. Characteristics of the MGR-GT turbomachinery. 





Turbine 


Compressor 


Stages 


6 


15 


RPM 


10000 


10000 


Tip Diameter (cm) 


86.2 


73.3 


Bladed length (cm) 


63.0 


130 


Poly. Eff. (Expected) 


.931 


.937 


Poly. Effi (minimum) 


.91 


.91 
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5.3.2 Scaling Relatiunships. 

In Reference 8 Staudt derives a series of scaling relationships for turbomachines. The 
relationships allow one to study the effects of changing plant parameters (pressure, mass flow 
rate, etc.) on turbomachine size. Once a ’good’ design is found, good estimates can be made 
of other machines with different system parameters. The basis of the equations that follow 
are that blade velocity triangles and blade stress remains constant. 

The first set of scaling ratios relate the blade tip diameter (D,) and ipm (N) to changes 
in system pressure, power level held constant: 



The second set of scaling ratios relates the blade tip diameter and rpm to mass flow rate 
m (power level) and number of stages (n), pressure held constant: 




5-1. 




5-2. 




5-3. 




5-4. 







5-5. 
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As can be seen from the above that increasing the system pressure reduces machine 
size and increases the speed. At constant pressure, lowering the mass flow and thereby the 
power level produces the same effect, smaller size and faster. 

The final scaling equation relates rpm, and number of stages (n) and mass flow 
between any two machines with the same inlet and exhaust temperatures and pressure ratios. 



5.3.3 Turbomachine Design Results. 

Using the above equations will produce turbomachines similar in performance to the 
baseline machine. This is in absence to any consideration of what is attached to the machine. 
Recalling from Chapter 3 the expected cycle conditions for the marine version, the pressure 
and temperature were the same as the MGR-GT but the mass flow was only 27.4 kg/s verses 
150 kg/s in the MGR-GT. The results using two different scaling methods (same number of 
stages, and set RPM) are listed below. 

Table 5-6. Results of using scaling relationships to size marine MGR-GT turbomachinery. 




5-6. 



= constant 



Turbine 



Compressor 



Stages 



6 



20 



15 



72 



RPM 



23307 



7200 



23307 



7200 



Tip Diameter (cm) 



36.0 



96.1 



31.3 



81 
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From the above it can be seen that both designs produced by scaling relationships were 
unacceptable. The first design is nice and compact, but at 23,000 rpm, designing an accept- 
able generator becomes very difficult. The second design produces a marginally acceptable 
turbine, but the compressor at 72 stages would be too long and unwieldy to be used. 

Tlte conclusion is that the design parameters that formed the basis of the MGR-GT 
design would not be acceptable in the smaller version. A solution would be to reduce system 
pressure, this is also unacceptable. A lower pressure would greatly increase the core pressure 
drop and heat exchanger size. Discussions with several mechanical engineers have led me to 
the conclusion that an acceptable turbine/generator package can be designed, but it would 
take a dedicated design that is beyond the scope of this work. For the purposes of the rest of 
the analysis I will assume the turbine-compressor package is the same size as the MGR-GT 
and rotating 7200 rpm. It is probably not a good assumption that the turbomachinery can be 
the same size at 7200 rpm as it was at 10,000. However, several suitable 7200 rpm genera- 
tors were found in the literature, and none which had a higher rotational speed. This does not 
imply that a suitable generator does not exist, just that it was not found. 

To estimate the size of the two turbines in the split shaft design it was assumed the two 
turbines were one turbine with the shaft split at the proper location. The high pressure (HP) 
turbine must power the compressor and the ships service generator. This is a total power 
requirement of 18.8 MW. The power turbine must power only the propulsion generator, for a 
total requirement of 16.0 MW. Using the above assumptions and the characteristics of the 
MGR-GT turbomachinery, the marine variant wUl have the following characteristics. 
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Table 5-7. Characteristics of the marine MGR-GT turbomachinery. 



Turbine 





Compressor 


Power 


Compressor 


Stages 


4 


3 


15 


RPM 


7200 


7200 


7200 


Pressure* (MPa) 


8.1 


5.1 


8.2 


Temperature (*C) 


850 


717 


30 


Tip Diameter (cm) 


68.4 


86.2 


73.3 


Bladed length (cm) 


42.0 


31.5 


130 


Poly. Eff. (minimum) 


0.9 


0.9 


0.9 



’ Pressure is at turbine inlet or compressor outlet. 
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Chapter 6 Electrical Design. 

This chapter details a possible electrical system design. A full electrical analysis of the 
system including optimization and component design was not performed. The intent of this 
chapter is to discuss technology and the broad requirements that the electrical system would 
have to meet. Possible equipment configurations will be discussed along with advantages 
and disadvantages found. 

6.1 Design Considerations. 

In designing the power system the following design consideration were adopted: 

• Generators must rotate at the same rpm as the prime mover. (No reduction gears 
allowed) 

• Ships service power must be ’clean’, ie. none or few harmonics, constant voltage and 
frequency independent of ship maneuvering. 

• Minimum size and weight. 

• No rotating mechanical seals subject to full system helium pressure are allowed, (to 
reduce helium leakage) 

• Components should require no more than ’medium-risk’ in development. Thus super- 
conducting motors and generators will not be considered. 

• Both the propulsion bus and ship service bus will operate at the same frequency. This 
will allow propulsion generator to power ship service loads if necessary. 

6.2 Integrated Electric Propulsion. 

There are several methods to transfer power from the turbines to the ship. They include 
direct drive, direct drive through reduction gears, direct energy conversion, and integrated 
electric propulsion. Direct drive can be ruled out immediately because the slow speed of the 
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rotating machinery (~180 rpm) would require huge turbomachinery. Direct drive through 
reduction gears would allow both the turbines and propeUors to operate at their most efficient 
RPM. This method was also rejected because of the weight of the reduction gear and support 
equipment plus the difficulty of producing effective high pressure helium mechanical seals*. 
Direct energy conversion of the scale required has not been demonstrated, and in keeping 
with the philosophy to minimize risk, rejects this method. This leaves integrated electric 
drive. A proposed single line diagram for a single shaft AIEP system is shown in Figure 6-1. 




Figure 6-1. Integrated electric propulsion plant system, single line diagram, [adapted 
from Ref. 22 ] 

*In commercial power plants much of the helium leakage can be traced to the mechanical 
seals. Also seal lubrication could represent a potential contamination source for the primary 
system. 
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In an integrated electric drive ship, both ships service power and propulsion are derived 
from the same prime mover, in this case, the closed cycle turbines. This method of propul- 
sion has been a source of great debate in recent years and it is not the intent of this paper to 
enter that debate. Advanced Integrated Electric Propulsion (AIEP) simply seems to be the 
best method for this particular choice of prime mover. Although electric propulsion has been 
used since before World War II, recent advances in power electronics and electric machine 
design has sparked renewed interest in this propulsion form. Some of the more significant 
advantages of AIEP are listed in Table 6-1. 

Table 6-1. Electric propulsion system benefits. [22] 



• Arrangement flexibility for propulsion plant components 

• Variable reduction ratio between the prime mover and thruster 

• Increased control flexibility 

• Reversibility for unidirectional prime movers without gears or CRP propellers 

• Ease of electrical cross current capability 

• Distributed prime mover operating hours 

• Survivability Improvements 

• Facilitates ships service power form propulsion prime mover 

• Ease of automation 

• Ability to create a totally enclosed power system to reduce helium loss 



There are three main components that must be considered. The generators, which con- 
vert mechanical rotation to electric power; the power conversion equipment, solid state fre- 
quency converters to provide maneuvering and reversing capability to the third component. 



Electrical Design. 



86 



the main propulsion motor. The motor that will be discussed is for either single shaft ships or 
twin shaft, twin reactor plant systems. It will not be appropriate for aircraft carriers or any 
other ship which requires more than one power module per shaft. 

6.3 Generator. 

The design of the system generators depends on several factors. Tliey must be cooled, 
be compatible with the prime mover, and they must generate power that the system can use. 
The more efficient the cooling system the smaller the design. With the advent of reliable, 
efficient, solid state, frequency converters and power conditioners the generator can, for the 
most part, generate any convenient frequency and the power converters can do the rest. As 
discussed in section 5.3 the optimum turbine design requires a very high rpm, therefore the 
generators will be designed with the highest rpm possible. 

The optimum use of equipment and space would be to have one generator coupled to a 
single turbine on each power loop. These two generators would power a single ships power 
bus, which would supply all ships loads, propulsion or ships service. In practice this has not 
been possible because of the harmonics and transients generated when the ship maneuvers. 
Modem electronic systems require cleaner power than is possible with this arrangement. 
Transients propagate too quickly for power conditioners to catch when the systems are linked 
electrically. With the arrangement shown in Figure 6-1 the only linkage between the two 
systems is thermodynamically through the reactor. This creates a natural buffering system 
which smooths out propulsion transients. The system therefore requires two different gener- 
ator designs; a 2.5 MW ships service generator, and a 17 MW propulsion generator. 

6.3.1 Generator Cooling. 

The cooling system for the generators has to effectively cool the generators without 
allowing either primary system helium to leak to the environment or outside agents, such as 
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water or oil, to leak in. This can be accomplished by one of two methods. The first is to cool 
the generator with high pressure helium. The second is to use water cooling in a pressurized 
casing. 

Helium Cooled. Large commercial generators have cooled with hydrogen for since 
the 50’s. The hydrogen cooling technology is well understood and could be applied directly 
to a helium cooled design. Helium has only slightly worse heat transfer characteristics than 
hydrogen and pumping power and windage losses would be similar. With the generator cas- 
ing pressurized to near Brayton cycle system pressure, the rotating seal between the generator 
and compressor or turbine could be very simple and cooled by gas leakage. This was the 
method proposed for the MGR-GT [8] Some disadvantages of this system are that it requires 
a separate helium system with cooler and circulator within the generator casing, and it results 
in a larger, heavier machine than the all water cooled machine. 

Water Cooling. Recently (70’s) techniques have been developed using water to 
remove the heat generated in both the stator and rotor. Water cooled machines have been 
used by central station utility generators for several years to increase power density. Figure 
6-2 shows the improvement in power density realized by various cooling techniques verses 
the year of introduction. The increase in power density afforded by water cooling is quite 
dramatic. Another plus is that NAVSEA has determined that the risk of water cooled motors 
and generators is acceptable, thus clearing the way for their use aboard ship. [22] 
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Figure 6-2. Historical generator power density trends. [22] 

Water cooling in no way prevents the casing from being sealed and pressurized with 
helium to prevent leakage around the mechanical seal between the generator and the prime 
mover. The risk of water ingress into the primary system is minimal for two reasons. The 
first is the water will be at a lower pressure than the helium so any leakage will be helium 
leaking into the water not the other way around. The second is that the generator is separated 
from the primary system by a gas cooled seal around the shaft. The only way water could 
leak into the primary system is for both the generator cavity and the primary system to loose 
pressure and the generator cavity fill to the level of the seals. The helium cooled design also 
has the potential for water ingress since the helium cooler is integral to the generator casing. 

Because of the significant weight and volume savings water cooled stators and rotors 



will be employed in this design. 
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6.3.2 Propulsion (ienerator. 

Each propulsion generator has to provide 16 MW. It would have to rotate at 23,000 
rpm to optimize the turbine. To design an acceptable generator at such a liigh speed is prob- 
ably possible but no previous design could be found to support this. The ’best’ rpm would 
have to be a trade off between the turbines and the generators. Several designs with similar 
power levels at speeds up to 7200 rpm were found and it was decided to use one of these. 
This decision puts a severe constraint on the turbine design, as mentioned earlier, but I feel 
designing a slower speed high-pressure helium turbine would pose less technical difficulties 
than a high speed motor design. A good compromise will probably be in the range of 1 1 to 
15 thousand rpm. The generator design chosen comes from a design program written by Jim 
Davis of MIT [1] His program performs an optimization analysis of water-cooled, electric 
motors and generators for use in ship propulsion. He analyzed 25,775 HP machines with 
speeds ranging from 180 rpm to 7200 rpm. It was one of the 7200 rpm design that were cho- 
sen. Table 6-2 lists the generator characteristics. 
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6.3.3 Ship Service Generator. 

Because of the small size of the ship service generator water cooling adds unnecessary 
complexity to the design. The ship service generators will therefore be helium cooled. This 
makes the ship service generator almost the same size as the propulsion generator. Design 
characteristics were taken from Reference 24 which describes the Toshiba super motor drive 
system. The difference between motors and generators are generally a matter of application, 
not design. The motor listed is only intended to get a idea of the size and weight not to be a 
detailed application. Its characteristics are also listed in Table 6-2. 

Table 6-2. Design characteristics of the Marine MGR-GT generators. 



Machine type 

Number of Pole pairs 
Shaft RPM 
Power (HP/MW) 
Synchronous frequency (Hz) 
Rotor radius (m) 

Active length (m) 

Overall length (m) 

Overall diameter (m) 

Weight (kg) 

Full load efficiency 



Ship Service 


Propulsion 


Synchronous 
Helium Cooled 


Synchronous 
Water cooled 


4 


4 


7200 


7200 


3400/2.5 


25,775/19.2 


480 


480 




.204 




1.08 


2.0 


1.94 


1.0 


.57 


4000 


3136 


.97 


.98 



I'PNlW 
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6.4 Propulsion Motor. 

The propulsion motor is the largest mechanical component of the entire system. The 
requirement for high power with low speed generally requires a large diameter motor. Again 
with water cooling technology the motor can be made smaller. A design for a 60 Hz 40,000 
HP motor is described in reference 23 The characteristics of this motor is listed in Table 6-3. 
Table 6-3. Design characteristics of the Marine MGR-GT Propulsion Motor. 



Machine type 
Shaft RPM 
Power (HP) 

Synchronous frequency (Hz) 
Overall length (m) 

Overall diameter (m) 

Weight (kg) 

Full load efficiency 



Synchronous-Water cooled 



0-180 

40,000 

60 

4.9 

3.81 

60100 

.96 



6.5 Power Conversion Equipment. 

In recent years there has been a revolution in power electronics. Rectifiers, controlled 
converters, inverters, and cycloconverters based on liquid cooled thyristor stacks, have all 
been developed to the point where there use in a shipboard environment will involve little 
risk. U.S. manufacturers have supplied many power conditioners to applications in industry 
at similar power levels for many years therefore little problem is seen in performing the 
power conversion necessary to control the motor or provide ship service power. In reference 
23 a water cooled frequency changer/power converter for a surface ship application. The 
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characteristics of the power converter is given in Table 6-4. Each power convener module 
consists of six semi-conductor devices arranged in a stack. Each thyristor is located between 
two liquid cooled heat sinks. The cooling system uses deionized water to control corrosion. 
Each of these units can handle 15 MW so only one unit is required per module. 

Table 6-4. Power Converter Characteristics. [23] 



Type: Converter-Inverter with DC link, line commutated. 

Rating: 15 MW, 6300V ac supply 

Weight: 22,700 lbs. (10300 kg) 

Auxiliaries: Deionized cooling water; uninterruptable power supply 

L - 260" (6.6 m) 

W-54" (1.4 m) 

H - 84" (2.1 m) 



Size: 
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Chapter 7 Control and Control Systems. 

This chapter will briefly discuss the methods of control available for this power plant. 
Reactor control will be discussed first then power plant control and speed control. 

7.1 Reactor Control. 

The primary method of reactor control will be through the negative temperature coeffi- 
cient. This will ensure that reactor power will follow load over the entire power range. 
Reactivity control is needed at start up, to set reactor outlet temperature and to compensate 
for bum-up and Xenon build up. Once temperature is set and the reactor is at steady state it 
will automatically and quickly follow load changes while maintaining basically constant tem- 
perature. 

Because of the height of the reactor core and pressure vessel control rods could not be 
used. For that reason, reflector/absorber drums are installed around the edge of the reactor 
core. The drums extend the length of the core and are 40 cm in diameter. They are con- 
structed of nuclear graphite in a light steel container. The steel suppons the graphite and pro- 
vides a wear surface for the drum supports. The current design uses 16 drums spaced evenly 
around the edge of the core. Figure 4-8 and 4-9 show the proposed lay out. This design is 
for illustration only. Reactor control calculations were not performed because it was not 
deemed necessary at this point in the design. The control system shown illustrates some fea- 
tures which will probably be necessary once more a more detailed analysis is performed. 
Table 7-1 lists some of these design considerations. 

Table 7-1. Reactor plant control system design considerations. 

• Since the core is batch refueled a large reactivity control margin is needed to compen- 
sate for bum-up and allow control at end of life. 
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• The drums are relatively large so the reactivity difference between positions is signifi- 
cant. 

• The drive mechanism is in a pressurized housing with no mechanical seals, this is to 
prevent helium leakage. 

• The control drum are cooled by the compressor discharge flow which sweeps the reac- 
tor vessel. 

• Control drives should be designed so that power is required to keep the drums in the 
’out’ position. In that way if control power is lost the drums will automatically rotate 
into the core, (passive scram on loss of control power) 

• Power peaking will be high in an edge controlled reactor. If later studies show unac- 
ceptable peaking with peripheral control drums, in core control rods will be necessary. 
Possibly inserted from the side. Material selection for in core rods will be difficult 
given the high temperatures. 

7.2 Power Plant Control. 

Since the reactor heat source provides an effectively constant turbine inlet temperature, 
the only two effective control methods available are turbine bypass or inventory control. 

Both methods work by adjusting the mass flow rate through the system. Bypass control is 
effective at all power levels, however since the heat source sees basically a constant load, 
bypass control is very inefficient during off design point operations. Inventory control, on 
the other hand, reduces mass flow by lowering system pressure. The advantage is plant effi- 
ciency remains relatively constant over the full range of inventory control. The major disad- 
vantage is that it requires a large storage volume to contain the removed inventory. For any 
practical system both bypass and inventory control systems must be used. 
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The range of control of each system depends on the volume available for control inven- 
tory vessels, the speed of response desired, Control valve size, and the operating profile. The 
operating profile is important because it determines how much time the system spends at off 
design points. According to GENSPECS a typical naval vessel is expected to have the fol- 
lowing speed-time profile while at sea. 

Table 7-2. Speed-time profile for naval vessels. 

Speed Time at or below 

speed (%) 

1/4 9% 

1/2 23% 

3/4 38% 

Full 25% 

Flank 5% 

Since the ship spends most of its time at 3/4 power and 68% of the time above the 50% 
power level, 50% is a reasonable lower limit for inventory control with bypass control below 
50%. This wUl require inventory control vessels capable of storing half of the full system 
inventory. A passive transfer system (Bleed off at compressor discharge, injection at 
compressor inlet) would reach its storage capacity when storage vessel pressure is equal to 
compressor discharge pressure. Since most of the primary volume is at or near compressor 
discharge pressure this method of inventory control would require a storage volume approxi- 
mately equal to the entire primary system volume to achieve 50% control. This is obviously 
unacceptable in a volume limited ship design. 

Required storage volume can be reduced if helium is pumped out of the system using 
transfer compressors. This allows the inventory control vessels to be at a higher pressure, 
thus reducing volume. The power necessary to mn the transfer compressors is the control- 
ling factor in this design. For a power decrease, the system is operating at a higher power 
level than is required. This power can be used to run the transfer compressors. For a power 
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increase, the storage vessels will be at a higher pressure than the primary system so no pump- 
ing power is needed until the primary system and the storage vessels equalize. A storage sys- 
tem pressure of twice the system pressure at 50% power level, reduces the required storage 
volume by half. Although the detail of such a process have not been worked out the trends 
are at least in the right direction. The capacity and response time for such a system will have 
to be the subject of a future study. 
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Chapter 8 Design Summary. 

The previous chapters detailed some of the design decisions that went into this marine 
version of a commercial power plant. This chapter summarizes those designs and attempts to 
put the whole system in perspective. The major components and overall plant perfonnance 
and capabilities will be discussed. Since the major intent of this study was to see what had to 
be done to put the plant to sea, a formal comparison with current naval nuclear power plants 
will not be done. Final conclusions and recommendations for future work will be made in 
the next chapter. 

8.1 Reactor Compartment Arrangement. 

The reactor compartment was constrained from Chapter 2 to be a cylinder 10 m in 
diameter and 20 m long. The reactor compartment had to completely contain the primary 
system and yet meet aU the other design constraints. Figures 8-1 thru 8-3 show the proposed 
reactor compartment arrangement. Figure 8-1 is a section through the reactor vessel showing 
the reactor, shielding tank and a shielded side passageway. Figure 8-2 is a top view showing 
a possible arrangement of the turbine-generator modules, the reactor and shielding. Figure 
8-3 is a side profile of one of the machinery modules. Figure 8-3 also shows one of the pos- 
sible flow arrangements that allow all pressure retaining surfaces to be swept with ’cool’ 
helium. 

The original intent was to make the system modular so that it could be transported to 
the buUding ways and installed on basically any ship. Because of the reactor plant size and 
weight of shielding this goal cannot be achieved. Also, since the shielding water tank needs 
to be in thermal contact with the sea it must be custom fit to each individual ship. What can 
be done though, is to put the turbomachinery together as a unit in their own containment ves- 
sel. This module can be buUt, transported, and installed as a unit ready to be piped into the 
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reactor. 




Figure 8-1. 



Section through Reactor Compartment. (Submarine) 
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Figure 8-2. Top view of marine MGR-GT Reactor compartment. 
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Figxire 8-3. Side profile of marine MGR-GT RC showing helium flow path. 
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8.2 Component Summary 

This section summarizes the major component designs; 

Reactor. The reactor is an 80 MWth helium-cooled pebble bed reactor. The major 
characteristics are: 

• Passive safety. The system is capable of withstanding a loss of coolant accident with- 
out system damage, by rejection of decay heat into the shielding water tank and ulti- 
mately into the sea. 

• Coolant pressure at 8.2 MPa with a core exit temperature of 850°C. 

• Coolant flow path which sweeps the pressure vessel with relatively low temperature 
helium. This cools the pressure vessel and allows current ASME pressure vessel 
codes to be used. 

• Potential for quick batch refueling without major disassembly. 

• High negative temperature coefficient at all power levels promotes stable operation 

• Reflector/absorber drums used for reactivity control. 

• Size: Vessel is 6.7 m high and 3.7 m in diameter. 

Shielding. Reactor shielding consists of a 1.8 m water tank which acts as a neutron 
shield. The reactor vessel is surrounded by a 25 cm thick lead gamma shield. The shield is 
designed to limit exposure to .5 mR/hr at the shield surface. 

Electric Plant. The electric plant is an integrated ship propulsion system that provides 
both propulsive power and ships service power from the same prime mover. Power is gener- 
ated at high frequency (480 Hz) and is converted to 60 Hz using solid state frequency conver- 
ters and power conditioners. The main propulsion motor (40,000 HP) is driven by two 
propulsion generators through a water-cooled, thyristor based, frequency converter. Both the 
main motor and the propulsion generators use water cooled stators and rotors to reduce size 
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and weight. The ships service generators are cooled by high pressure helium. The electrical 
system provides flexibility in that ships propulsion generators can provide ships service 
power if required. Table 8-1 summarizes the electrical power plant components. 

Closed Drayton Cycle. The system is powered by a highly recuperated Closed Bray- 
ton Cycle. The equipment configuration chosen is a split shaft system with the compressor, 
ship service generator, and the high pressure turbine on one shaft, and the low pressure, or 
power turbine, and propulsion generator on the other shaft. All components are enclosed in a 
common casing and are supported by active magnetic bearings. The split shaft arrangement 
decouples ship service power generation from the propulsion bus, and it facilitates reactor 
start up by allowing the ship service generator to act as a motor to power the compressor, 
thus circulating helium until the reactor gets to power. The power for start-up would be pro- 
vided by the emergency generators or shore power. The system heat exchangers feature low 
specific pressure drop, <3%, and high effectiveness. The regenerator is a compact plate-fin 
arrangement while the precooler is a shell and tube heat exchanger to facilitate repairs and 
inspections. The Brayton cycle components are also summarized in Table 8-1. 
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Table 8-1. Marine MGR-GT Plant Parameters and Equipment Summary. 

T urbomachinerv 



Stages 


Turbine 

HP 

4 


Power 

3 


Compressor 

15 


RPM 


7200 


7200 


7200 


Pressure* (MPa) 


8.1 


5.1 


8.2 


Temperature (*C) 


850 


717 


30 


Tip Diameter (cm) 


68.4 


86.2 


73.3 


Bladed length (cm) 


42.0 


31.5 


130 


Poly. Eff. 


.9 


.9 


.9 


Surface (ref 11) 


Heat Exchangers 

Precooler 

S-1.50-1.25 


Regenerator 

1/9-24.12 


Volume (m^) 


1.08 




(both sides) 
3.75 


Length (m) 


.756 




1.5 


Frontal Area (m^) 


1.43 




2.5 


Weight (kg) 


1500 




10400 


Relative Pres. Drop 


.0025 




.0252 


Effectiveness 


.93 




.949 



Electrical Components 





Generator 


Propulsion 




SS 


Prop. 


Motor 


Cooling 


He 


Water 


Water 


Number of Pole Pairs 


4 


4 




Shaft RPM 


7200 


7200 


0-180 


Frequency 


480 


480 


60 


Power (HP/MW) 


3400/2.5 


25,775/19.2 


40,000 


Length (m) 


2.0 


1.94 


4.9 


Diameter (m) 


1.0 


.57 


3.8 


Weight (kg) 


4000 


3136 


60100 


Full load efficiency 


.97 


.98 


.96 



Pressure is at turbine inlet or compressor outlet. 
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8.3 Weight Summary. 

This section summarizes the weight estimates for the reactor compartment components 
studied above. These weights are for the component only. They do not include foundations, 
supports, piping, insulation, valves, or the weight of fluid and support equipment. 

Table 8-2. Reactor Companment Component Weight Summary. 



Component 


Number Installed 


Total Weight 
(tons) 


Remarks 


Reactor 


1 


105 




Shielding 


1 


490 


Lead-water shield com- 
pletely surrounding reactor. 


Regenerator 


2 


20.8 


10.4 each 


Precooler 


2 


3 


1.5 each 


Turbomachinery 


2 


26 




Subtotal 




644.8 





As can be seen above the estimated weight of the reactor plant components alone is 
over 600 tons. Most of this is reactor shielding. This weight must be reduced in order to 
have a viable design. 
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Chapter 9 Conclusions and Closing Remarks. 

There have been many studies of using nuclear gas turbines on board ship. They have 
all concluded that the system warranted further study, but the technology just was not to the 
point where it was feasible at the time. This particular design can say the same thing to a 
point. The technology to support a closed Brayton cycle and electric propulsion either exists 
or will exist in the near future. The ability to create low-pressure drop, highly effective com- 
pact heat exchangers and the continued refinement and development of magnetic bearings 
removes the major disadvantages of the closed Brayton cycle that many previous studies 
have focused on. Advances in power electronics brings effective electric propulsion technol- 
ogy within reach. Everything is in place except the heat source. 

As this design stands the reactor heat source is too big and too heavy to be practical on 
smaller ships (submarines, frigates, and destroyers.) While it is true that some weight and 
volume is saved outside of the reactor, any gain is more than offset by the weight and size of 
the core and its required shielding. In Reference 17 Hsu states that an acceptable maximum 
weight for a reactor subsystem (core and shielding) is 300 tons. Assuming this is tme the 
marine MGR-GT is over 200 tons too heavy. Unless the reactor can be made smaller this 
system will not be suitable for shipboard use. 

I feel the major driver in the reactor size is passive safety (including water ingress). 
While there is no question that passive safety in a civilian power plant is worth almost any 
design compromise to achieve, the same thing may not be true of naval power plants. 
Although passive safety should certainly be a goal, it cannot be sought at the expense of all 
else. The bottom line is that any good ship design is the result of compromise. The system 
must be safe, but it must also work. 
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9.1 Areas for Future Study. 

Most of the future study needs to be on size reduction for the reactor heat source and 
concurrently the radiation shield. The following summarizes the areas that effect core size 
and some possible solutions. 

Reactor Size. Reactor criticality studies were performed using a two group approxi- 
mation and the NGC program. The cross sections used in the calculations came from the 
MHTGR. The MHTGR is a five region cylindrical reactor with an annular core and using 
prismatic fuel. Although the flux averaging used to calculate the cross sections cannot be the 
same for this design they were considered close enough for preliminary' studies and to get 
rough sizes. Before more detailed studies could be conducted more accurate cross sections 
applicable to pebble bed reactors need to be developed. 

Enrichment. Cross sections used in NGC corresponds to an enrichment of 7%. Once 
more accurate cross sections are developed, the effects of various enrichments can be stu- 
died. The goal here is to reduce reactor size and to ensure enough excess reactivity is loaded 
at beginning of life to last until the next refueling and to provide Xenon override. 

Water Ingress. The water ingress problem is the one potential Achilles’ Heel in the 
system. Much further smdy is required in this area, for if the core cannot tolerate a core 
flooding casualty, it cannot be considered a feasible design. 

Fast Core. Of all the possible solutions to the problem of core size I feel that going to 
a fast spectrum core could be the best solution. Fast cores are small and very power dense 
and generally have a high bumup and a long life. Water ingress is not a problem as far as 
reactivity is concerned, if the spectmm is softened by water ingress the reactor will simply 
shut down instead of the reactivity insertion caused in the Pebble bed core. The main ques- 
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tion that needs to be answered is whether such a core can be made passively safe. Without 
the thermal mass of the graphite temperature changes will probably be very rapid. The 
question is can heat be radiated away rapidly enough to prevent fuel damage. 

Plant Control. More work needs to be done on how inventory control on a system 
with more than one turbine-compressor group would work. Questions to be answered are: 
what are the response times, can the two systems be operated independently or must they 
always portion the load equally, what is the optimum volume and pressure of the storage ves- 
sels, and the mass capacity of the compressors to get the desired response. 

High-Speed Generafor/Lower Speed Turbine. It was discussed in Chapter 5 that 
applying the turbine scaling equations produced either turbomachines that were too large or 
rotated too fast to design a reasonable generator without using reduction gears. More study 
needs to be done on optimizing the generator-turbine combination. There should be some 
combination of rpm, number of stages, and system pressure that produces an acceptable gen- 
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Appendix A Heat Transfer Program. 

The depressurized loss of coolant accident was analyzed using the program 
HEAT.BAS. HEAT.BAS is a general purpose code that solves the one-dimensional transient 
heat conduction problem for any axi-symmetric cylindrical body. It is written in Quick- 
BASIC Version 4.5 and runs on personal computers. The program is based on the method 
outlined by Sanchez in his report Pu.ssive Heat Removal: Sensitivity Study for Modular 
Pebble Bed Reactors. [7] 

A.l Program Theory and User Guide. 

The program uses the modified explicit method to solve the transient problem. In this 
method the object to be analyzed is partitioned into radial nodes. A heat balance between 
nodes is performed and the change in temperature for each node at the desired time step is 
calculated. This process is repeated for each node at every time step. After each time step 
the nodal temperatures are saved and become the initial conditions for the next iteration. 

The program assumes that the outermost node is at constant temperature and the center- 
line temperature is the same as the temperature of the first node. This corresponds to a con- 
stant temperature exterior boundary condition and a zero heat flux central boundary 
condition. This sets the boundary conditions for the problem. 

The Equation A-1 is the heat balance equation used in the normal explicit method. It 
basically states that the temperature at any node j at one time step in the future is the temper- 
ature it is now plus the change in temperature due to the net heat added to or removed from 
the node. There are three heat addition terms. The first two are the heat flow from the two 
adjacent nodes and the third is the heat generation term. 

T/r + 60 = T/0+-^^{K;_,[T._,(r)-T/0] -hK;:„[T.„(0-T/r)] +Q.) 



A-1. 
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is the effective heat transfer coefficient from node i to node j, Cy = (pC^V), and 

contains the material properties of node j. Qj is the heat generation rate and includes the 
decay heat fraction f(t), t in seconds. 

Q; = 9c/(0Vy A-2. 

fit) = ,0622[r^-(r + 3.184 x 

Vy = 7r(r;-r;.,) A-3. 



Two types of heat transfer between nodes are allowed, conduction, and radiation. The 
program handles this by using two different methods to calculate the effective linear thermal 
conductance between nodes (K/). The program detennines the type of heat transfer from the 
input file. 

K; for conductive heat transfer between nodes is estimated using the following formula 



[ 10 ]: 



K^ = 



27ik 

In- 



A-4. 



which is the formula for heat conduction through a hollow cylinder, k is the thermal 
conductivity of node], fj and ij are the radii of nodes i and j. 

For radiation heat transfer between nodes (such as across a helium gap) the equation for 
radiative heat transfer between concentric cylinders is used. (Equation A-5) [10] 




iKvp T?-Tj 



A-5. 
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where r is nodal radius, a is the Stefan-Boltzmann constant, and 0 denotes the eniissivity of 
the surface. Since radiation is a surface effect, the effective volume of a radiation node is the 
not the volume between the nodal radii (Equation A-3) but it is based on the formula 
(equation A-8). The desired time step is substituted for t„„ and the equation is solved for 
volume. This method ensures a stable solution and models the heat transfer more accurately. 

The modification to the explicit method uses the following equation to achieve a stable 
solution using larger time steps then would normally be allowed by the nomial explicit 
method. [7] 

Tjit + 5r) = T//) + (K;: _ ,(T. ^.(/) - T/0) 

A-6. 

+ k; , ,(T. „(/) - T//)) + Q .) + Z.(T/r) - T/r - 5r))) 

where Tj(t) is the temperature at node j at the current time, Ty(r + 5t) is the temperature at 

node j one time step in the future, Tj{t — 5t) is the temperature at node j one time step in the 
past. 

Zj is a factor which is supposed to provide a stable solution at any time step. However, 
during trial runs in HEAT.BAS the range of stability was increased by only a factor of about 
two. For example, if a two second time step was stable using the unmodified method then a 
four second time step was usually stable with the modified equation. Zj is a smoothing factor 
which performs a weighted average of the temperature change over more than one time step. 
This smoothing factor does not provide absolute stability but it does delay the onset of insta- 
bility and reduces the oscillation amplitude after instability develops. 

Zj is defined as follows: 
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A-7. 



0 




A-8. 



is the time step which will provide a stable solution. It is a function of the material 
properties of the node and the nodal volume Vj. 

A.2 Materials. 

HEAT.BAS uses the following codes to designate the material for each node. The pro- 
gram uses the material code to calculate the temperature dependent thermal property for each 
node. Where correlations or mathematical models are available they are used. If no equation 
is available, tabular data is used and the properties are determined by linear interpolation for 
the desired temperature. 

Table A-1 List of material codes and material used in HEAT.BAS. 



The method of calculation and equations used for each of the materials listed in Table 



Material Material 
Code 



0 Core pebbles 

1 Graphite 



2 304 stainless steel) 

3 2 1/4 Cr Mo steel 



4 Helium 

5 Reflector to core interface 

6 Test Material (Constant Material properties) 



A-1 are as follows. 
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Core Pebbles. The core pebbles’ thermal conductivity is calculated using the modified 
Zehner-Schuluender model. This model is valid for high temperatures and includes radiation 
between the pebbles and conduction through the pebbles. [7] 



Kff- 



(i-VT^)<))+ 
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1 ^ 1 




1 B. J 


' ' 0-kJ 
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Cp = (l-(l))C*+<})-C^' 


A- 10. 




A-11. 



where 

T absolute temperature in *K. 

()> porosity of pebble bed. 

E emissivity of pebbles. 

CT Stefan-Boltzmann constant. 

Cp specific heat in J/kg*K. (Superscript indicates either graphite or helium), 
effective thermal conductivity in W/m-*K. 

"kf thermal conductivity of pebble (graphite) in W/m-'K. 

p density in kg/m^. (subscript indicates either graphite or helium). 

Graphite, 304 Series Stainless Steel and 2 1/4 Cr Mo Steel. The heat transfer properties of 
graphite and steel are provided in tabular form. The program uses a subroutine to look up the 
specific heat and the thermal conductivity of the material as a function of temperature. Lin- 
ear interpolation is used to obtain values between the tabulated temperatures. Densities were 



Thermal Conductivity (W/m-K) 
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considered constant at 1700 kg/m^ for graphite, 7800 kg/m^ for 304 stainless steel, and 7675 
kg/m^ for CrMo steel. Figures A-1 and A-2 show the specific heat and thermal conductivity 
of the various materials as a function of temperature. [7] 




O Graphite + CrMo Steel O 304 SS 



Figure A-1 



Thermal conductivity of Graphite, 304 Stainless Steel and 2 1/4 Cr Mo Steel 
as a function of temperature. [7] 



Specific Heat (kJ/Kg-K) 
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□ Graphite + CrMo Steel O 304 SS 



Figure A-2 Specific Heat of Graphite, 304 Stainless Steel and 2 1/4 Cr Mo Steel as a 
function of temperature. [7] 

Helium. The thermal properties of helium at atmospheric pressure were calculated using the 
following correlations and equations. [7] 

K(P,T) = .002682(1 + .001 123P)T‘''‘'-'^*’' A-12. 

C^ = 5195.J/kg°K A-13. 



P(P,T) = 



48.14(PA') 




A-14. 






Heat Transfer Program. 



117 



where: 



K(P,T) 



P(P,T) 



= pressure in bars. 

= absolute temperature in 'K. 

= specific heat in J/kg’K. 

= thermal Conductivity in W/m-'K. 
= density in kg/m^. 
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A.3 HEAT.BAS Input File. 

The program works with an input file that must provide the following information for 



each node. 

1 . The outer radius of the node in meters. 

2. The initial temperature in *C. 

3. The material code. (See Table A-1) For nodes with conductive boundary condi- 
tions the code is for the material to the inside of the nodal radius. If the node is a 
radiating surface (such as the outer boundary of the pressure vessel) the code is 
for the material of the surface not the helium. 

4. The emissivity on the interior surface of the node. (0 is used for a conductive 
boundary condition.) 

5. The emissivity on the outer surface of the node. 

6. The initial heat generation rate (power level) of the node in W/m^ 

Table A-2 shows a sample input file. The file is an ASCII file with the entries either in 
columns (as shown in Table A-2) or separated by commas. The nodal radii do not have to be 



in order but no two nodes are allowed to have the same radius. 



Table A-2. 


Sample 


input file for HEAT.BAS 






0.05 


582 


0 


0 


0 


8210000 


0.1 


580 


0 


0 


0 


8150000 


0.2 


579 


0 


0 


0 


8100000 


0.4 


572 


0 


0 


0 


8010000 


0.5 


565 


0 


0 


0 


7600000 


0.6 


560 


0 


0 


0 


7300000 


0.7 


552 


0 


0 


0 


7250000 


0.8 


544 


0 


0 


0 


7240000 


0.9 


542 


0 


0 


0 


6800000 


1.01 


539.3 


0 


0 


0 


6670000 


1.06 


538.5 


0 


0 


0 


7330000 


1.12 


537 


0 


0 


0 


7000000 


1.24 


536 


0 


0 


0 


6870000 


1.3 


531 


0 


0 


0 


6900000 


1.35 


530 


5 


0 


0 


3000000 


1.45 


425 


1 


0 


0 


0 


1.55 


320 


1 


0 


0 


0 


1.57 


310 


2 


0 


0.6 


0 


1 . 77 


200 


3 


0.6 


0 


0 


1.8 


185 


3 


0 


0 


0 


1.85 


170 


3 


0 


0.6 


0 


1.87 


50 


3 


0.6 


0 


0 
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A.4 Source Code Listing for HEAT.BAS 



f ★ 



• ★ 



one dimensional heat conduction for a cylindrical core 



The program requires the initial temperature and power distribution 
** and solves for the centerline temperature after shutdown with no 
' * normal heat removal 






DECLARE SUB DATAPREP () 

DECLARE SUB HEATTRAN (K! , rl I , rl. Tempi!, Temp!, el!, e2 ! , qlj!, jKl ! ) 
DECLARE SUB HePROP (t!, p! , rho ! , Cp ! , Kf) 

DECLARE SUB MPROP (t I , Cp ! , K! , prop ( ) AS ANY) 

DECLARE FUNCTION INTERP! (rl, a!, bl) 

DECLARE FUNCTION Keff! (t ! , por!, emis!, d!, K! ) 

DECLARE SUB SORTNODE () 

DECLARE SUB MATPROP (n%, t!, rho!, Cp ! , K! ) 

DECLARE FUNCTION DECAYPOWERl (tl) 



^♦♦^♦initialization section 



DEFINT I-N 

CONST pi - 3.141592654#, S * 5.6697E-08, TRUE - -1, FALSE 
hepres - 1; rhog - 1700: por - .39: d - .06: emis - .8 



TYPE infilerec 



radius 

inittemp 

mat type 

bcin 

bcout 

initpower 



AS SINGLE 
AS SINGLE 
AS INTEGER 
AS SINGLE 
AS SINGLE 
AS SINGLE 



END TYPE 



TYPE matrec 



Temp 


AS 


SINGLE 


K 


AS 


SINGLE 


Cp 


AS 


SINGLE 



END TYPE 



0 



♦♦♦♦♦ dimension arrays 
DIM SHARED graph (1 TO 16) AS matrec 

DIM SHARED ss304 (1 TO 10) AS matrec, CrMo (1 TO 8) AS matrec 
DIM Temp(lOO), Told(lOO), Tnew(lOO), heat (100) , r(lOO) 

DIM mtype(lOO) AS INTEGER 

DIM bcin (100), bcout (100), v(lOO) 

DIM inrec AS infilerec 
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/ ★★★★★Load material property arrays **** 

FOR i - 1 TO UBOUND (graph) : ' Graphite 

READ graph (i) . Temp, graph(i).K, graph(i).Cp 
NEXT i 

FOR i -= 1 TO UBOUND (ss304) : ' Stainless steel 

READ ss304 (i) . Temp, ss304 (i).K, ss304 (i).Cp 
NEXT i 

FOR i - 1 TO UBOUND (CrMo) : ' Cr Mo Steel 

READ CrMo (i) . Temp, CrMo(i).K, CrMo(i).Cp 
NEXT i 

SCREEN 0, 0, 0: CLS 

PRINT "ONE DIMENSIONAL TIME DEPENDENT HEAT TRANSFER PROGRAM” 

PRINT ”FOR A CYLINDRICAL CORE.” 

PRINT 

t ★★★★★ load problem ***** 

name$ - "input. 000” 
out$ - "heat. out” 

PRINT "WHAT IS THE NAME OF THE INPUT FILE <”; name$; ”> ”; 

INPUT X$ : PRINT 

IF LEN(x$) <> 0 THEN name$ - x$ 

OPEN name$ FOR INPUT AS #1 

f ★★★★★ find number of nodes ***** 

FOR n - 0 TO 100 

IF EOF(l) THEN EXIT FOR 
LINE INPUT #1, r$ 

IF r$ - "" THEN EXIT FOR 
NEXT n 
CLOSE 1 

OPEN name$ FOR INPUT AS #1 

f ★*★★★ load node properties ***** 

FOR i - 1 TO n 

INPUT #1, r(i), Temp(i), mtype(i), bcin(i), bcout(i), heat(i) 

Temp(i) - Temp(i) + 273.15 
NEXT i 

Temp(O) - Temp(l) 

PRINT "What is the name of the output file<”; out$; ”> 

INPUT X$ 

IF x$ «* ”” THEN x$ ■* out$ 

OPEN X$ FOR OUTPUT AS #2 

INPUT "What is the time step"; dt 

CALL SORTNODE: ' ★*★★★ sort nodes by radius 

CALL DATAPREP : ' ***** find effective nodal volumes 

f *************i(***********it************************************************* 

' ****** input file loaded and sorted -- start calculations 

t - 0: ' ★★★★★set initial time 



***** start of time loop 



***** 
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DO: ' ★★★★★ start with the exterior node and work toward the center 

t » t + dt : ' ★★★★★ elapsed time in seconds 

time » t / 3600: ' convert to hours 

f ★★★★★ check for a radiation outer boundry and get 

' heat transfer coefficients for outer node 

IF bcin(n) - 0 THEN q - .5 * (Temp(n - 1) + Temp <n) ) ELSE q * Temp (n) 

CALL MATPROP (mtype (n) , q, rho, Cp, K1 1 ) 

HEATTRAN K1 1 , r(n - 1), r (n) , Temp(n - 1), Ternp(n), bcout (n - 1), bcin(n), q2j, jK2I 
q2j - -q2j 

' ★★★★ loop thorugh nodes 

FOR j - n - 1 TO 1 STEP -1 
il - j - 1 
IF bcin(j) - 0 THEN 
/ ***** Conduction only ***** 

CALL MATPROP (mtype ( j) , .5 * (Temp(il) + Temp(j)), rho, Cp, K1 I ) 

HEATTRAN Kll, r(il), r(j), Temp(il), Temp(j), bcout(il), bcin(j), qlj, jKl I 
ELSE 

f Radiation node ***** 

CALL MATPROP (mtype ( j) , Temp ( j ) , rho, Cp, Kl I ) 

HEATTRAN Kll, r<il), r(j), Temp(il), Temp(j), bcout(il), bcin<j), qlj, jKl 1 
END IF 

q * heat(j) * DECAYPOWER (t ) * v(j): ' ★★★** calculate decay heat 

f ★★★★★ calculate new temperature for node j 

IF jKl I + jK2l - 0 THEN tmax - dt + dt ELSE tmax - rho ^ Cp * v(j) / (jKl I + jK2 1 ) 
IF dt / tmax > 1 THEN Zj - .5 * (dt / tmax - 1) ELSE Zj - 0 

Tnew(j) - Temp(j) + (dt * (qlj + q2j + q) / rho / Cp /v(j)) / (1 + zj) + Zj * 

(Temp( j) - Told( j) ) 

jK2 I - jKl I : ' ***** setup for next node 

q2j - -qlj 
NEXT j 

* ***** of node loop *** 

r***************ffir********************************************************** 

r ****** New temperature distribution has been found ***** 

' print the distribution once an hour to an output file 

' and check for temperature reduction to end program 

Tnew (n) “ Temp (n) 

Told(n) * Temp(n) 

FOR i « 1 TO n 

Told(i) - Temp(i) 

Temp(i) « Tnew(i) 

NEXT i 

' print every hour 

IF FIX (time) <> FIX (time - dt / 3601) THEN 
PRINT #2, time; 

FOR i - 1 TO n: PRINT #2, USING "#####.#”; Tnew(i) - 273.15; : NEXT i 

PRINT #2, 

PRINT time, Tnew(l) - 273.15, Tnew(l) - Thold, telaps - TIMER: Thold - Tnew(l) 



Heal Transfer Program. 



122 



telaps - TIMER 
END IF 

Told(O) - Temp(O) 

Temp(O) - Tnew(l): ' approximate central boundry condition 

LOOP UNTIL Told(O) >- Temp(O) AND time > 1 

^ ^nd of time Xoop*****^********^********* *****^***^******* 

f Central temperature has started to drop end program 

PRINT 

PRINT "Max Centerline temperature is”; Told(O); ” at t time; "hours” 

PRINT #2, time; 

FOR i - 1 TO n: PRINT #2, USING "#####.#"; Tnew(i) - 273.15; : NEXT i 



► Qf niain program* 

***** Data atatamant® for graphite ** 



/ 

DATA 


T 

473.15, 


120 


K 


Cp 

, 1150.00 


DATA 


573.15, 


110 


.2543 


, 1344.455 


DATA 


673.15, 


99. 


89896 


, 1507.8 


DATA 


773.15, 


90. 


76186 


, 1620.885 


DATA 


873.15, 


84 . 


06132 


, 1721.405 


DATA 


973.15, 


77. 


96992 


, 1771.665 


DATA 


1073.15 




72.48766, 1834.49 


DATA 


1173.15 


f 


67. 61454, 1872.185 


DATA 


1273.15 


f 


62.74142, 1935.01 


DATA 


1373.15 


r 


58. 47744, 1972.705 


DATA 


1473.15 




54.21346, 2010.4 


DATA 


1673.15 


f 


49.34034, 2060.66 


DATA 


1873.15 


f 


46.90378, 2098.355 


DATA 


2073.15 




43.85808, 2123.485 


DATA 


2273.15 




40.81238, 2148.615 


DATA 


2473.15 




38.37582, 2173.745 


$ 


★ ★★★★ 


data 


atataxnanta for 304 ai 


! 

DATA 


T 

273.15, 


K 

14.8 , 


cp 

470.00 


DATA 


373.15, 


16. 


199, 


504.4492 


DATA 


473.15, 


17. 


7383, 


531.1502 


DATA 


573.15, 


19. 


2043, 


552.0922 


DATA 


673.15, 


20. 


6703, 


567.2752 


DATA 


773.15, 


22. 


063, 


578.7933 


DATA 


873.15, 


23. 


3824, 


591.3585 


DATA 


973.15, 


24. 


6285, 


602.8766 


DATA 


1073.15 


, 25. 


9479, 


617.536 


DATA 


1173.15 


, 27. 


194, 


634.2896 


t 


***** 


data 


atatamanta for 2 1/4 


$ 

DATA 


T 

273.15, 


36.2 


K 


cp 

440. 


DATA 


373.15, 


37.0661, 


475.384 
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DATA 


473.15, 


37.2755, 


516.217 


DATA 


573.15, 


37.0661, 


558.6205 


DATA 


673.15, 


35.2862, 


585.319 


DATA 


773.15, 


33.8204, 


660.703 


DATA 


873.15, 


32.1452, 


751.792 


DATA 

END 


973.15, 


30.0512, 


915.124 



StTB DATAPREP 

' Data prep subprogram calculates effective nodal volume 

SHARED Tempo, Told(), heat () , r(), mtype() AS INTEGER, n 
SHARED bclnO, bcoutO, v(), dt 

FOR i - 1 TO n - 1: ' ★★★★★ start loop through nodes 

il - i - 1: ipl - i + 1 

Told(i) - Temp(i): ' Initialize temp history array 

IF bcin<i) <> 0 THEN 

/ ♦’^’^♦♦find effective volume of a radiation boundry surf ace 

CALL MATPROP (mtype (1) , Temp(l), rho, Cp, K1 I ) 

HEATTRAN Kll, r(il), r(i), Temp(il), Temp(i), bcout (11) , bcln(l), qlj, jKl I 
CALL MATPROP (mtype (ipl) , Temp (ipl), rhol, Cpl, K2 I ) 

HEATTRAN K2I, r(i), r(ipl), Temp(l), Temp(ipl), bcout (i) , bcin(ipl), q2j, jK2 I 
v(i) - (jKl I + jK21) / rho / Cp dt 1 . 25 

v(i + 1) - -v(i) : ' ★★★★★ subtract volume from next node 

ELSE 

' Conduction node 

v(i) - v(i) + pi (r(i) 2 - r(il) 2) 

END IF 

PRINT i, r (i) , v(i) 

NEXT i 

Told(n) - Temp(n) 
t ★★★★★ print file header 

FOR i - 0 TO n 

PRINT #2, i; : NEXT i: PRINT #2, 

FOR i - 0 TO n 

PRINT #2, r(i); : NEXT i: PRINT #2, 

PRINT #2, 0!; 

FOR i - 1 TO n 

PRINT #2, USING *’#####.# Temp(i) - 273.15; : NEXT i: PRINT #2, 

END SUB 

FUNCTION DECAYPOWER (t) STATIC 

t Function calculates the decay power fraction as a function 

' of time in seconds 

IF t < 1 THEN 

DECAYPOWER - EXP ( -2 . 8 0 94 8 9 ♦ t) 

ELSE 

DECAYPOWER - .0622 ♦ (t (-.2) - (t + 3.184E+07) (-.2)) 

END IF 

END FUNCTION 
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SUB HEATTRAN (K1 , rl, r, T^mpZ, T^xapO , •!, m2, qj, JKl ! ) STATIC 

★★★★★ calculate the heat transfered from adjecent node 
input variables 

kl Thermal conductivity of node 

rl radius of inner node 

r radius of node 

Tempi temperature of inner node 
Tempo temperature of node 

el boundry condition into node (emissivity) 

e2 boundry condition out of inner node (emissivity) 

output variables 

qj heat transfered from inner node 

jKl 1 effective linear heat transfer coefficient 



IF Tempi - Tempo - 0 THEN 

qj - 0: jKl 1 “0: ' ★★★★★ no temperature difference between nodes 

ELSE 

IF e2 - 0 THEN 

' ***** Conduction only 

qj - 2 * pi * Kl / LOG(r / rl) * (Tempi - TempO) 

ELSE 

f radiation boundry with conduction through helium 

'radiation first 

qj - 2 * pi * rl * S / (1 / el + (1 - e2) / e2 / r * rl) * (Tempi ^ 
CALL HePROP ( . 5 * (Tempi + TempO), hepres, x, y, Khl) 
qj - qj + 2 * pi * Khl / LOG (r / rl) ♦ (Tempi - TempO) 

END IF 

jKl! “ qj / (Tempi - TempO) : 

END IF 
END SUB 

SUB HePROP (t, p, rho, Cp, Kl ) STATIC 

' ★★★★★ Helium properties as a function of temperature and pressure 

' input variables 
' t Temperature in 'K 

' p pressure in bars 

' output variables 
' rho density in kg/m^3 

' K! thermal conductivity 

' Cp Specific heat 

rho - 48.14 * (p / t) / (1 + .446 * p / t ^ 1.2) 

K! - .002682 * (1 + .001123 * p) * t ^ (.71 * (1 - .0002 * p) ) 

Cp - 5195 
END SUB 

FUNCTION Keffl (t, por, emis, d, K! ) STATIC 

f This function computes the effective thermal conductivity 

' through a pebble bed using the modified Zehner- 

' Schulender equation 

' input variables 



- Tempo 4 ) 
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f 


t 


Temperature *K 


t 


por 


Porosity of pebble bed 


0 


emis 


emisssivity of pebbles 


0 


d 


diameter of spheres 


0 


K! 


Thermal conductivity of graphite at temp t 



Bz - 1.25 ((1 - por) / por) (10 / 9) 

b-4*S*d*t^3 
gammaf — K! / b / pi 
a * 2 / emis - 1 
c - SQR(1 - por) 

K«ff! - b* ((1-c) / a + c/ a’^ (Bz + 1) / Bz / (1 + 1/ a/ gammaf)) 

END FUNCTION 

SOB MATPROP (n, t, rho, Cp, K! ) STATIC 

/ Subroutine to find the various material properties given temp 

' the following code is used 
' N— 0 Core pebbles 

' 1 graphite reflector 

' 2 core barrel (currently 302 stainless) 

' 3 Pressure vessel (currently 2 1/4 Cr Mo steel 

' 4 helium 

' 5 Reflector to core 

' 6 Test Material (Constant Material properties) 

SHARED por, emis, d, hepres, rhog 
SELECT CASE n 

CASE 0: ' core pebbles 

' to find )ceff for the pebble bed the porosity , emissivity, and 
' pebble diameter are defined in the main module 
CALL MPROP(t, Cpg, Kg 1 , graph ()) : ' graphite properties 

K! - Keffl (t, por, emis, d. Kg!) 

CALL HePROP(t, hepres, rhoh, Cph, Kh!): ' Helium properties 

rho « rhog (1 - por) + rhoh * por 
Cp » (1 - por) * Cpg + por * Cph 
CASE 1: ' graphite 

CALLMPROP(t, Cp, K!, graphO) 

CASE 2: ' 302 stainless 

CALLMPROP(t, Cp, K! , ss304 ()) 

rho - 7800 

CASE 3: ' 2 1/4 Cr Mo steel 

CALL MPROP(t, Cp, Kl, CrMo()) 
rho - 7675 
CASE 4: ' heliiam 

CALL HePROP (t, hepres, rho, Cp, Kl ) 

CASE 5: ' reflector to core interface node 

CALL MPROP (t, Cpg, Kg! , graph ()): ' graphite properties 

pi - por / 2 

K! * Keff! (t, pi, emis, d. Kg!) 

CALL HePROP (t, hepres, rhoh, Cph, Kh!): ' 



Helium properties 
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rho “ rhoh * pi + rhog * (1 - pi) 

Cp * (1 - pi) * Cpg + pi * Cph 
CASE 6: ' t#at material (constant material properties) 

K! - 50 
Cp - 2000 
rho * rhog 
CASE ELSE 
PRINT n 
STOP 

END SELECT 

END SUB 

SUB MPROP (t, Cp, Kl, propO AS matreo) STATIC 

! ««««« Finds the material properties as a function of temp 

n - UBOUND (prop) 
tl - prop(l).Temp 
IF t < tl THEN 

Cp « prop(l) .Cp 
Kl “ prop (1) . K 
EXIT SUB 

ELSEIF t >- tl AND t <- prop (n). Temp THEN 
FOR i - 2 TO n 
il - i - 1 

IF prop(i) .Temp >- t THEN 

ratio - (t - tl) / (prop(i) .Temp - tl) 

Cp - ratio * (prop(i).Cp - prop(il).Cp) + prop(il).Cp 
Kl - ratio * (prop(i).K - prop(il).K) + prop(il).K 
EXIT SUB: * <- normal exit 

ELSE 

tl — prop(i).Temp 
END IF 
NEXT i 
END IF 
BEEP 

PRINT ’’ERROR, Temperature out of range of tabulated values” 

PRINT "Temperature t; ” Allowable range is”; prop (1 ). Temp; "to"; prop(i).Temp 

END 

END SUB 

SUB SORTNODE 

f this subroutine sorts the nodes in order of radius 

' it also checks for duplicate radii and gives a warning 

r 

SHARED Tempo, heat(), r(), mtype() AS INTEGER, n 
SHARED bcinO, bcout() 

FOR i - 0 TO n 
FOR j « i + 1 TO n 

IF r(j) - r(i) THEN 
BEEP 
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PRINT "Input file contains two nodes with the seme radius”; r(i) 
PRINT "program cannot continue." 

END 
END IF 

IF r(i) > r(j) THEN 
SWAP r (i) , r < j) 

SWAP Temp ( i ) , Temp ( j ) 

SWAP heatU), heat<j) 

SWAP mtype(i), mtype(j) 

SWAP bcout(i), bcout(j) 

SWAP bcin(i), bcin(j) 

END IF 
NEXT j, i 
END SUB 
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Appendix B Heat Exchanger Analysis. 

This section outlines the methods used to design the main system heat exchangers, the 
regenerators and the precoolers. The regenerators were analyzed using the program 
COMPHX.BAS and the precoolers were designed using the program PRECOOL.BAS. Both 
programs were written in QuickBASIC version 4.5 and run on IBM XT type personal com- 
puters. 

B.I Regenerator Design Program COMPHX.BAS 



COMPHX.BAS perfonns a complete analysis of a compact, plate-fin, counter flow 
heat exchanger with cross flow headers. It performs this analysis using the effectiveness- 
NTU method [11]. Where effectiveness is defined by either equation B-1 or B-2 and NTU 
(Number of heat Transfer Units) is a non-dimensional parameter defined in equation B-3. 



£ = 



~C , 



T,. -Tu 



B-1. 



T,. -r. 



r. -r. 



B-2. 



Tu -T. 



Where C is a flow stream capacity rate defined as C = riiCp, m is the mass flow rate. 



Cp is the specific heat, A is the minimum flow area, and U,„ is an average heat transfer 
coefficient. 



NTU = 



C . 

mtn 



B-3. 



COMPHX.BAS was derived from a program written by Jon Ness as described in 
Reference 12. It was originally written in FORTRAN IV and was used to predict the 
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performance of recuperators installed on open cycle marine gas turbines. I have rewritten the 
program in QuickBASIC to mn on personal computers and have added the abUity to model 
non-uniform hot side gas flow conditions, either helium or a fuel-air mixture as the working 
fluid, and input in either metric or English units. 

B.1.1 Mel hod of Analysis. 

There are two different types of heat exchanger analysis. The first (or design problem) 
takes as input parameters the desired effectiveness, mass flow rate, and pressure drop and 
produces as output the size required. The second method (or analysis problem) is the inverse 
of the first method. In this method the size, mass flow rate, and inlet conditions are given 
and the heat exchanger performance is calculated. The choice of method depends on the spe- 
cific problem, the first method tends to be better in the initial design phase, while the second 
is better at predicting the performance of a given design. COMPHX.BAS uses the second 
method. 

The analysis procedure is summarized below: 

1. Define cycle conditions. Fluid types (helium or air-fuel), mass flow rates (tii), cold 

side inlet pressure (F^..^) and temperature hot side inlet temperature and out- 
let pressure along with flow velocities define the problem initial conditions. 

2. Select heat exchanger properties. The user selects the counter flow length L, total 
frontal area A^^, material properties (metal density p and thermal conductivity k), plate 
thickness a, and core matrix fin geometr>'. The fin geometry is specified separately for 
the hot and cold sides since different fin geometries could be used on each side. All 
necessary surface characteristics (such as plate spacing b, hydraulic radius r^, fin thick- 
ness 5, ratio of heat transfer area to volume between plates (compactness) p, and ratio 
of fin area to total area, are from reference 1 1 and are listed in Appendix C This 
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3 . 



4. 



5. 



completes the problem definition. 

Set up iterations. Calculate the heat transfer and free flow areas on both sides, and 
make initial guesses for effectiveness and pressure drop. These initial guess will be 
iterated through the following steps until the processes converges on the fmal values. 
Determine temperature dependent fluid properties. Average hot side and cold side tem- 
peratures are used to calculate fluid properties. The assumed effectiveness is used to 
calculate temperatures based on equation B-1 and fluid properties are calculate using 
correlations in references 1 1 and 13. 

Calculate Reynolds number on both sides. Reynolds number is defined as: 



N = 



4/-.G 



B-4. 



where G is the flow stream mass velocity G = rii/A^ and A,, is the free flow area on one 
side. 

6. Determine Stanton number Nst, Colburn factor j = N^Np^, and friction factor f from 

the tabulated data from reference 1 1 . The program has tabulated data for Colburn and 
friction factor as a function of Reynolds number for each of the surfaces in its data 
base. The desired values at the required Reynolds number are interpolated between the 
tabulated values. 

7. Calculate the heat transfer coefficient h using the equation 

/i=NSTG-Cp B-5. 

8. Calculate the fm effectiveness Tj^d surface effectiveness Tjo using; 
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tanhw/ 




B-8. 






9. Calculate the overall coefficient of heat transfer U, based on the cold-side area. 



1 1 



1 



a 
■+- 



B-9. 



Uav T{oA (AJA,)r\ohhH k 



10. Calculate the new NTXJ and effectiveness. For a counter flow heat exchanger effec- 
tiveness is related to NTU by the following [1 1][12]: 



e = 



1 — e 



NTU(I-C 1 



B-10. 



1 — C;.e 



-NTU(1-C,) 



In the case where C, = 1, such as a closed cycle gas turbine with no bypass flow, the 
above reduces to the following. 

NTU 



e = 



1 -NTU 



This new value for effectiveness is compared with the initial guess. If the effectiveness 
has changed the new value is used as the initial guess and the process starts again at 
step 4 above. 
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11. Once effectiveness has converged the pressure drop is calculated. The pressure drop is 
composed of three major components: the header pressure drop, core pressure drop, 
and the bend pressure drop. These pressure losses are dependent on the geometry of 
the heat exchanger so for this program the configuration shown in Figure B-1 is used. 
Although the final heat exchanger may not look exactly like B-1 it should be close 
enough for this stage in the design process. 




Figure B-1 Regenerator Arrangement Used for Pressure Drop and Weight Estimates. 

[ 12 ] 



Again the calculated pressure drop is compared to the initial estimate. If it has changed 
the new pressure drop is used as the initial estimate and the process goes back to step 4. 
There is therefore two major loops; an outer pressure drop loop and an inner effective- 
ness loop. When both loops are satisfied the fluid conditions are set and all tempera- 
tures and pressures are known. 
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Of the three components of pressure drop the core loss is the most complicated. The 
following formula is used to estimate the core drop; 















--1 



V, 






A, V, 



ij 



B-12. 



The terms in brackets represent the four elements of the core drop; (1) entrance effects, 
(2) flow acceleration, (3) core friction, and (4) exit effects. In the above the K factors are 
entrance and exit loss coefficients, v, and Vj is the fluid specific volume at the entrance and 
exit, o is the ratio of free flow area to frontal area, and f is the friction factor from the 
tabulated data. 

Pressure loss in the crossflow headers is handled by increasing the core loss by assum- 
ing the headers effectively make the core longer. The headers are triangular wedges so the 
average length is easily obtained. 

Finally, the bend loss occurs where the fluid has to change direction as it moves from 
the headers into or out of the core. It is estimated by assuming it is a miter elbow with loss 
coefficient Kt,. Kb is curve fit from the literature and is used in the same manner as the 
entrance and exit coefficients above. [12] 

Afrer the above is completed the program estimates the weight of the heat exchanger. 
This estimate is done in five parts: (1) weight of cold side fins, (2) weight of hot side fins, (3) 
total weight of separating plates, (4) enclosure weight, and (5) header weights. 

The weight estimate for the fins and plate is based on the geometry of the fins (from the 
surface data base), the density of the metal, and core volume and frontal area. From the 
above data an average core density is determined. Multiply this density by the core volume 
to obtain core weight. Header density is assumed to be the same as the core so header weight 
is density times header volume. 
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Finally, the enclosure weight is based on the surface area. The area density of the 
enclosure is assumed to be 15 Ib/ft^ and it includes metal, insulation, and supports. [12] This 
number obviously will not be good for every application, but it was based on marine gas tur- 
bine installations so it should be close. 
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B.1.2 Input file for COMPHX.BAS 

The input file is an ASCII file which contains all the input data necessary to run the 
program. It can be created with a pure ASCII text editor or the easist way is to run the pro- 
gram and input the data manually. COMPHX.BAS will save the input data as a data fUe. 

The file contains the following information: 

Line Data 

1 Comment 

2 Type of units (1 = English, 2 = Metric) 

3 Cold side surface type, surface number, fluid code (1 = air-fuel, 2 = helium) 

4 Hot side surface type, surface number, fluid code 

5 Metal density, thermal conductivity, and plate thickness. 

6 Frontal area, counter-flow length. 

7 Cold side mass flow rate, inlet temperature, inlet pressure, and header velocity. 

8 Hot side mass flow rate, exit pressure. 

9 Number of horizontal nodes, number of vertical nodes. 

10-? Hot side velocity and inlet temperature for each node. 

Last Cold side and hot side fuel-air ratio. 
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Table B-1 Sample Input File for COMPHX.BAS 

STRIP FIN SURFACE 1/9-24.12 
2 

3,17,2 

3,17,2 

7832.8,16, .1 
2.5, 1 .5 

27.4,417.1,8.2,30 

27.4.4.01 

1.1 

1,878.7 

0,0 
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B.1.3 Sample Output from COMPHX.BAS 

COMPACT HEAT EXCHANGER ANALYSIS 

RUN # 1 INPUT FILE 3_17.IN 

1x1 NODES 

STRIP FIN SURFACE 1/9-24.12 



CORE HEAT TRANSFER SURFACE 




COLD-SIDE 


HOT-SIDE 


FIN TYPE 




STRIP 


STRIP 


FIN AND SURFACE NUMBER 




3-17 


3-17 


Fins/cm 




9.496044 


9.496044 


FIN LENGTH (cm) 




4.374007E-02 


4.374007E-02 


PLATE SPACING (cm) 




.1905004 


.1905004 


HYDRAULIC RADIUS (cm) 




3.025146E-02 


3.025146E-02 


FIN THICKNESS (cm) 




1.016002E-02 


1 .016002E-02 


COMPACTNESS (mVcu m ) 




2831.365 


2831.365 


FIN/TOTAL AREA (mVm^) 




.857 


.857 


FLUID TYPE 




HELIUM 


HELIUM 


LENGTH = 1.5 m 




VOLUME = 


3.749999 cu m 


FRONTAL AREA - 2 . 5 m^ 




HEIGHT « 


2.377325 m 


METAL DENSITY - 7832.8 


kg/cu 


m WEIGHT - 


10448.79 kgs 


PLATE THICFCNESS - . 1 cm 




K * 


16 W/m-*K 


COLD HEADER DESIGN DETAILS 












INLET 


EXIT 


DIAMETER (m ) 




.3504919 


.5268331 


VELOCITY (m /sec) 




30 


27.41736 


HEAT EXCHANGER CONDITIONS 












COLD-SIDE 


HOT-SIDE 




INLET EXIT 


INLET EXIT 


PRESSURE (MPa ) 


8.2 


8.141543 


4.082453 4.01 


TEMPERATURE (*K) 


417 


.1 855.1212 


878.7 440.6789 


MASSFLOW ()cg/sec) 




27.4 


27.4 


PRESSURE DROP (%) 




.7128874 


1.806785 


PRESSURE DROP (in H20) 




234.7874 


296.2564 


HOT SIDE WALL TEMPERATURE 


TK) 




867.5189 429.4602 


REYNOLDS NUMBER 






1134.089 1822.023 


Tw-Tg (*K) 






11.18117 11.21867 



EFFECTIVENESS = 94.89191 % 

TOTAL PRESSURE DROP - 2.519673 % 



NTU 



18.57682 
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B.1.4 Source Code Llsfiny for COMPMX.BAS 

' WRITTEN BY RICHARD D. LANTZ . 3-14-89 

' THIS PROGRAM PERFORMS COMPACT HEAT EXCHANGER PERFORMANCE AND 

' SIZING CALCULATION USING METHODS AND EMPIRICAL DATA FROM KAYS 

' AND LONDON'S "COMPACT HEAT EXCHANGERS" . THIS PROGRAM IS BASED 

' ON A PROGRAM WRITTEN BY JON NESS AS PUBLISHED IN DTNSRDC 

' REPORT PAS 82 - 41. 

' THE PROGRAM WAS WRITTEN IN QUICK BASIC AND IT ALLOWS 

' INTERACTIVE MODIFICATION OF INPUT PARAMETERS, THE USE OF AN INPUT 

' FILE, ALLOWS SAVING OF INPUT PARAMETERS AS AN INPUT FILE, ALLOWS 

' INPUT AND OUTPUT TO BE IN EITHER METRIC OR ENGLISH UNITS AND 

' ALLOWS A NON-UNIFORM GAS INPUT VELOCITY AND TEMPERATURE 

' DISTRIBUTION. THE WORKING FLUID IS EITHER FUEL-AIR OR HELIUM. 



k'kitit'k'k-kititit-k'kitititititiftiitit'lt'k-k-k'ltit'kit-kit’kitititiflt'kiflt-k'k-kitit-k-kit-kliitififitltit'tfk'k'ititific 



►INITIALIZATION SECTION* 



DECLARE FUNCTION TANH! (X!) 

DECLARE FUNCTION YN% () 

DECLARE FUNCTION MINI (XI, Yl) 

DECLARE SUB SURF (T%, NS%, S$ , AA! , BBl, SF 1 , PS 1 , RH I , DELI, BET I, FRI, WF I , WP 1 ) 
DECLARE SUB STAT (TYPE$, nn%, Al()) 

DECLARE SUB TRANSP (Tl, FARl , CP 1 , TK!, MUl, MWI, GTYPE%) 

DECLARE SUB INTERP (Al(), RE I , NSTI, FI) 

DECLARE SUB BENDLOS (XI, Yl) 

CONST PI - 3.14159, R - 15451, GC - 32.2 



CONST S 


V tl » 


TRUE - -: 


1, FALSE 


- 0, 


YES - 


-1, NO - 


0 














COMMON 


SHARED 


/TRANSPl/ 


Al, 


A2, 


A3, 


A4, 


A5, 


A6, A7, 


A8, 


Bl, 


B2, 


B3, 


B4, 


B5, 


B€ 


COMMON 


SHARED 


/TRANSP2/ 


Cl, 


C2, 


C3, 


C4, 


C5, 


C6, C7, 


Dl, 


D2, 


D3, 


D4, 


D5, 


D6, 


D7 


COMMON 


SHARED 


/TRANSP 3/ 


El, 


E2, 


E3, 


E4, 


E5, 


E6, E7, 


FI, 


F2, 


F3, 


F4, 


F5, 


F6, 


F7 



B6, B7, B8 



OPTION BASE 1 
DEFINT I-K 

DIM TA(2), TG(2), TW(2), TGRW(2), XNRA(2), XNRG(2) 

DIM SURFACE$(4), STATS$(4) 

DIM STATA(18, 4), STATG(18, 4) 

READ Al, A2, A3, A4, A5, A6, A7, A8 

DATA 1. 0115540E-25, -1 . 4526770E-21, 7 . 6215767E-18, -1 . 5128259E-14 
DATA -6.7178376E-12, 6 . 5519486E-08, -5 . 153 687 9E-05 , 2 . 5020051E-01 
READ Bl, B2, B3, B4, B5, B6, B7, B8 

DATA 7.2678710E-25, -1 . 3335668E-20, 1 . 0212913E-16, -4 . 2051104E-13 
DATA 9. 9686793E-10, -1.3771901E-06, 1.2258630E-03, 7.3816638E-02 
READ Cl, C2, C3, C4, C5, C6, C7 

DATA -6.2176401E-22, 7.1827364E-18 

DATA -3.1410386E-14, 6.7214720E-11, -7 . 5336781E-8 , 6 . 1 97 9074E-5 , -4. 81E-3 
READ Dl, D2, D3, D4, D5, D6, D7 

DATA 1 . 0404582E-19, -7 . 5213293E-16 




n r 
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DATA 2.1637607E-12, -3 . 15 930 96E- 9 , 2 . 4649233E-6, -9 . 0800204E-4, 1 . 1073E-1 
READ El, E2, E3, E4, E5, E6, E7 

DATA 2 . 4724974E-21, -1 . 675 6272E-17 

DATA 4 .1505396E-14, -3. 9906519E-11, - 9 . 13 47177E- 9, 8 . 8 743855E-5 , 2 . 98E-3 
READ FI, F2, F3, F4, F5 , F6, F7 

DATA -2 . 02551 69E-19, 1 . 4196996E-15 

DATA -3 . 9713025E-12,5 . 6582466E-9, - 4 . 3 41 4 613E- 6, 1 . 813500 9E-3 , -3 .3929E-1 

FOR I -= 1 TO 4: READ SURFACE$ (I) , STATS$(I): NEXT I 

DATA "PLAIN" , "PIAIN . STA" , "LOUVERED" , "LOUVERED . STA” 

DATA "STRIP", "STRIP. STA", "WAVY", "WAVY. STA” 



DIM C (2, 8) , 


- C$(2, 8) 








FOR I - 1 TO UBOUND(C, 2) 








FOR J - 


1 TO 2 








READ 


C(J, I), C$(J, 


I) 






NEXT J, I 












DATA 


1, 


"in". 




.3937, 


"cm" 


DATA 


1, 


"ft". 




3.28084, 


"m " 


DATA 


1, 






10.76391, 




DATA 


1, 


”cu ft”. 




35.31467, 


”CU 1 


DATA 


1, 


” • R” , 




1.8, 


"*K" 


DATA 


1, 


”psia" , 




145.0377, 


"MPa 


DATA 


1, 


"lb”. 




2.2046, 


"kg” 


DATA 


1, 


"Btu/hr-f t- 


*R", 


.578176051, 


"W/m- 



DIM FTYPE$ (2) 

READ FTYPE$(1), FTYPE$(2) 

DATA "AIR-FUEL", "HELIUM” 

READ RHO, Kl, ERRLIM, FTYPA%, FTYPG%, UT% 

DATA 489., 12., ,001, 1, 1, 1 

READ E, DELPA, DELPG, IRUN 
DATA .5, 0.01, 0.03, 1 

READ TYPA%, NSA%, TYPG%, NSG%, NV%, NH% 

DATA 3,14, 3,15, 1, 1 
READ RLENI, AFRA, FARG, FARA 

DATA 0.9275, 0.643, .0145, 0.0 

READ WA, PINA, TINA 

DATA 90.0, 116,4, 1040.5 
READ WG, PEXG, TING 

DATA 101.45, 14.70, 1646.4 
READ VINA, Tp 

DATA 90.0, .012 

DIM VING(NV%, NH%), MDOT(NV%, NH%) , TIN(NV%, NH%) 

SCREEN 0, 0, 0: WIDTH 80: CLS 

PRINT " WELCOME TO THE PERFORMANCE PREDICTION " 

PRINT " PROGRAM FOR COMPACT PLATE-FIN HEAT EXCHANGERS" 
PRINT " HEAT EXCHANGER TYPE TO BE SPECIFIED AS FOLLOWS- " 
PRINT " WHERE TYPE: " 

GOSUB FINTYPE: 
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PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 



" THE PROGRAM WILL. ACCEPT AN INPUT FILE OR THE HOT AND" 
" COLD SIDE PARAMETERS CAN BE ENTERED MANUALLY." 

” OTHER INPUT IS SELF-EXPLANATORY.” 

" THE PROGRAM IS BASED ON A COMPUTER MODEL FOUND IN" 

” DTNSRDC PAS 82-41.” 

" FOR A LISTING OF HEAT EXCHANGER SURFACE PARAMETERS" 

" RUN THE PROGRAM SURFLIST.” 



t ★★★mpQT SECTION 



TOP : 



PRINT 

PRINT " DO YOU WANT TO USE AN INPUT FILE? 

IF YN% * YES THEN 

PRINT "FILE NAME <"; DOC$; ”>”; 

INPUT F$ 

IF F$ <> THEN DOC$ - F$ 

ON ERROR GOTO FILEERRl 
OPEN DOC$ FOR INPUT AS #1 

PRINT "READING RUN # IRUN; " FROM FILE 
LINE INPUT #1, COMMENT? 

INPUT #1, UT%: ' 

INPUT #1, TYPA%, NSA%, FTYPA% : ' 

INPUT #1, TYPG%, NSG%, FTYPG%: ' 

INPUT #1, RHO, K!, Tp; ' 

INPUT #1, AFRA, RLENI: ' 

INPUT #1, WA, TINA, PINA, VINA: ' 

INPUT #1, WG, PEXG: ' 

INPUT #1, NH%, NV%: ' 

XMAX - 0! 



DOC$ 

type of units 

cold side surface 

hot aide surface 

heat exchanger properties 

frontal area and core length 

cold side conditions 

hot side conditions 

number of nodes 



' ** load hot side velocity, temp, and mass flow arrays 

REDIM VING(NV%, NH%) , TIN(NV%, NH%) , MDOT(NV%, NH%) 

FOR J - 1 TO NV% 

FOR I - 1 TO NH% 

INPUT #1, VING(J, I), TIN(J, I) 

VING(J, I) - VING(J, I) * C(UT%, 2): ' CONVERT TO ENGLISH UNITS 
TIN(J, I) « TIN(J, I) * C(UT%, 5) 

MDOT(J, I) • VING(J, I) / TIN(J, I) 

XMAX - XMAX + MDOT(J, I) 

NEXT I , J 
TING « 0! 

FOR J - 1 TO NV% 

FOR I - 1 TO NH% 

MDOT(J, I) - MDOT(J, I) * WG / XMAX ^ C(DT%, 7) 

TING - TING + TIN(J, I) * MDOT ( J, I) / C (UT%, 7) 



NEXT I, J 
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TING - TING / WG 

INPUT #1, FARA, FARG: ' fuel air ratio on both sides 

CLOSE 1 
ELSE 

DOC$ - ’’NONE" 

END IF 

********** IKPDT FILE HAS BEEN READ - ALLOW MODIFICATION 



DO 

PRINT 

PRINT "THE INPUT FILE HAS THE FOLLOWING LABLE" 
PRINT COMMENT$ 

PRINT "IS THIS CORRECT? "; 

I - YN% 

IF I - NO THEN 

PRINT "ENTER THE NEW LABLE" 

LINE INPUT COMMENT$ 

END IF 

LOOP WHILE I - NO 



DO 

IF UT% « 1 THEN A$ - "ENGLISH” ELSE A$ - "METRIC" 

PRINT 

PRINT A$; " UNITS ARE BEING USED. IS THIS CORRECT? "/ 

I - YN% 

IF I - NO THEN UT% - 2 \ UT% 

LOOP WHILE I - NO 

CALL SURF(TYPA%, NSA%, "COLD", AXA, BXA, SFA, BA, RHA, DELA, BETA, FRA, WFA, WPA) 
CALL SURF(TYPG%, NSG%, "HOT", AXG, BXG, SFG, BG, RHG, DELG, BETG, FRG, WFG, WPG) 

DO : CLS 

PRINT "THE CURRENT HEAT EXCHANGER MATERIAL PROPERTIES ARE AS FOLLOWS;" 

PRINT 

PRINT "DENSITY ", , RHO; C$ (UT%, 7); "/"; C$ (UT%, 4) 

PRINT "THERMAL CONDUCTIVITY ", K! ; C$ (UT%, 8) 

PRINT "PLATE THICKNESS", Tp; C$ (UT%, 1) 

PRINT : PRINT "IS THE ABOVE CORRECT? " ; : I - YN% 

IF I - NO THEN 
PRINT 

PRINT "DENSITY C$ (UT%, 7); "/"; C$ (UT%, 4); ")"; 

INPUT RHO 

PRINT "THERMAL CONDUCTIVITY ("; C$(UT%, 8); ")"; : INPUT Kl 

PRINT "PLATE THICKNESS ("; C$(UT%, 1); ")"; : INPUT Tp • 

END IF 

LOOP WHILE I - NO 

RHO - RHO * C(UT%, 7) / C(UT%, 4) 



Tp 


C(UT%, 


1) 


Kl 


* C(UT%, 


8) 
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DO: CLS 

PRINT ”THE HOT SIDE FLUID IS FTYPE$ (FTYPG%) 

PRINT "THE COOL SIDE FLUID IS FTYPE$ (FTYPA%) 

PRINT "IS THE ABOVE CORRECT? "; : I - YN% 

IF I “ NO THEN 
PRINT 

PRINT "THE CHOICES OF FLUID TYPE ARE:” 

U% - UBOUND (FTYPE$) 

FOR J - 1 TO U%: PRINT " J, FTYPE$ ( J) : NEXT J 

PRINT ; PRINT 
DO 

PRINT "THE HOT SIDE FLUID IS (1 U%; ”) : INPUT J 

LOOP WHILE J < 1 OR J > U% 

FTYPG% -= J 
DO 

PRINT "THE COLD SIDE FLUID IS U%; ”) ”/ : INPUT J 

LOOP WHILE J < 1 OR J > U% 

FTYPA% « J 
PRINT 
END IF 

LOOP WHILE I « NO 



DO: CLS 
PRINT 

PRINT ” THE CURRENT CORE SIZE IS AS FOLLOWS:” 

PRINT "FRONTAL AREA - "; AFRA; ” ("; C$ (UT%, 3); ")" 

PRINT "COUNTERFLOW LENGTH - RLENI; " (”; C$ (UT%, 2); ")”: PRINT 

PRINT "IS THE ABOVE CORRECT? "; : I - YN% 

IF I » NO THEN 

PRINT : PRINT "WHAT IS THE FRONTAL AREA "; C$ (UT%, 3) ; 

INPUT AFRA 

PRINT "WHAT IS THE COUNTERFLOW LENGTH "; C$ (UT%, 2); 

INPUT RLENI 
END IF 

LOOP WHILE I •= NO 

AFRA - AFRA * C(UT%, 3) : ' CONVERT TO ENGLISH UNITS 

RLENI - RLENI * C(UT%, 2) 



DO: CLS 

PRINT "THE CURRENT COLD SIDE CONDITIONS ARE AS FOLLOWS:”: PRINT 
PRINT "MASS FLOW RATE - 
PRINT "VELOCITY - 

PRINT "INLET TEMP - "; 

PRINT "INLET PRESSURE - 

PRINT "IS THE ABOVE CORRECT? : I -= YN% 

IF I « NO THEN 



WA; C$(UT%, 7) 


; ”/s" 


VINA; 


C$ (UT%, 


2); "/s" 


TINA; 


C$ (UT%, 


5) 


PINA; 


C$ (UT%, 


6) : PRINT 



PRINT 



PRINT "ENTER THE COLD SIDE CONDITIONS.’ 



PRINT 


"MASS FLOW RATE ("/ 


C$ (UT%, 


7) ; 


”/s) 


PRINT 


"VELOCITY ("; 


C$ (UT%, 


2) ; 


”/s) 


PRINT 


"INLET TEMPERATURE {"; 


C$ (UT%, 


5) ; 


H ^ fff • 



INPUT WA 
INPUT VINA 



INPUT TINA 





i 
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PRINT "INLET PRESSURE ("; C$(UT%, 6); : INPUT PINA 

END IF 

LOOP WHILE I - NO 
WA - WA ^ C(UT%, 7) 

TINA - TINA C(UT%, 5) 

VINA * VINA * C(UT%, 2) 

PINA •» PINA * C(UT%, 6) 



DO: CLS 

PRINT "THE CURRENT HOT SIDE CONDITIONS ARE AS FOLLOWS:": PRINT 



PRINT "TOTAL MASS FLOW RATE - 
PRINT "EXIT PRESSURE 
PRINT " IS THE ABOVE CORRECT? 

IF I - NO THEN 

PRINT : PRINT "ENTER THE HOT SIDE CONDITIONS 
PRINT "MASS FLOW RATE ("; C$ (UT%, 7); "/s)"; 

PRINT "EXIT PRESSURE ("; C$(UT%, 6); ")"; : 

END IF 

LOOP WHILE I - NO 

WG - WG * C(UT%, 7): ' ***CONVERT TO ENGLISH UNITS 

PEXG - PEXG * C(UT%, 6) 

XMAX - 01 



WG; C$(UT%, 7); "/s" 

PEXG; C$<UT%, 6); PRINT 

: I - YN% 



" : PRINT 
: INPUT WG 
INPUT PEXG 



CLS 

PRINT "THERE ARE NH%; " HORIZONTAL AND NV%; " VERTICAL NODES." 

PRINT : PRINT "IS THIS COF^ECT? 

IF YN% - YES THEN 
CLS 

PRINT "THE CURRENT HOT SIDE VELOCITY AND TEMPERATURE PROFILE IS AS FOLLOWS:" 
PRINT "NODE", "VELOCITY", "TEMPERATURE" 

PRINT , C$(UT%, 2); "/sec", C$ <UT%, 5) 

FOR J - 1 TO NV% 

FOR I - 1 TO NH% 

PRINT J; I; ")", VING ( J, I) / C(UT%, 2), TIN(J, I) / C(UT%, 5) 

NEXT I, J 
PRINT 

PRINT " WANT TO CHANGE ANY HOT SIDE VELOCITIES OR TEMPERATURES? 

IF YN% - YES THEN GOSUB GASVELPR02 
ELSE 

GOSUB GASVELPRO: 

END IF 

PRINT 

IF FTYPG% - 1 OR FTYPA% - 1 THEN 
DO: CLS 

PRINT "THE CURRENT HOT SIDE FUEL-AIR RATIO IS FARG 
PRINT "THE CURRENT COLD SIDE FUEL-AIR RATIO IS "; FARA 
PRINT : PRINT "IS THE ABOVE CORRECT? "; : I - YN% 

IF I - NO THEN 
PRINT 

INPUT "HOT SIDE FUEL-AIR RATIO IS? " ; FARG 
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INPUT ”COLD SIDE FUEL-AIR RATIO IS? FARA 
END IF 

LOOP WHILE I - NO 
ELSE 

FARG -= 01 
FARA « 01 
END IF 



***** CALCULATION SECTION***** 

***•«***«««**★* heat exchanger core dimensions ******************** 

CALL STAT ( ( STATS $ (TYPA%) ) , NSA%, STATA() ) 

CALL STAT ( (STATS$ (TYPG%) ) , NSG%, STATG()) 

CALL TRANSP (TINA, 01, CPA, KAl, MUA, MA, FTYPA%) 

RU - R / MA 

RHOINA - PINA / TINA / RU * 1441 

DINA - SQR(41 * WA / RHOINA / VINA / PI) 

♦♦*★**★★★★****★ HEAT TRANSFER AND FREE FLOW AREAS ***************** 

**** INITIAL LOOP VALUES**** 

VOL - AFRA * RLENI 



ALHA - 


BETA 


* BA 


/ 


(BA + 


BG + 


2! 


* Tp) 


ALHG - 


BETG 


* BG 


/ 


(BA + 


BG + 


21 


♦ Tp) 



AA - ALHA * VOL 
AG - ALHG * VOL 
ANODE - AFRA / NV% / NH% 

SIGA - ALHA * RHA 

SIGG - ALHG * RHG 

ACA - SIGA * AFRA 

ACG - SIGG * AFRA 

GA - WA / ACA 

ACANODE » ACA / NV% / NH% 

ACGNODE - ACG / NV% / NH% 

WANODE - WA / NV% / NH% 

FLA - VOL / AFRA 
XNCFL - SQR(AFRA) 

CMIN « 1! 

CA - 11 
CG » 11 

DELPT - DELPA + DELPG 
EOLD - 01 
DPT - 01 

START OF PRESSURE LOOP **** 
DO 

TBULKEA » 01 

TBULKEG - 0! 

DPT - DELPT 

PEXA “ PINA * (11 - DELPA) 
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PING - PEXG * (II + DELPG) 

FOR J - 1 TO NV% 

FOR I - 1 TO NH% 

★★★★★ START OF E-NTU LOOP 
DO 

FIND EXIT TEMPERATURES FOR EACH NODE 
TEXA - E * (TIN(J, I) - TINA) * CMIN / CA + TINA 

TEXG - E * (TINA - TIN(J, I)) * CMIN / CG + TIN(J, I) 

TAVG - (TEXG + TIN(J, I)) * .5 

TAVA - (TEXA + TINA) ^ .5 

***** AVERAGE CORE FLUID PROPERTIES 
CALL TRANSP(TAVA, FARA, CPA, KAI, MUA, MA, FTYPA%) 

CALL TRANSP(TAVG, FARG, CPG, KG I , MUG, MG, FTYPG%) 

**** AVERAGE CORE REYNOLDS NUMBERS 
GG - MDOT(J, I) / ACGNODE 
NRA - 4! * RHA * GA / MUA 

NRG - 4 ! * RHG * GG / MUG 

*****STANTON NUMBER, COLBURN FACTOR, AND FRICTION FACTOR 
CALL INTERP (STATAO , NRA, COLBFA, FA) 

IF FA - 01 GOTO ERRl: 

NPRA - CPA * MUA / KAl : ' prandel nuinber 

NSTA * COLBFA / NPRA .666: ' Stanton number 

HA - NSTA ♦ GA * CPA * 3600!: ' heat trans coeff 

CALL INTERP (STATG 0 , NRG, COLBFG, FG) 

IF FG - 0! GOTO ERRl: 

NPRG - CPG * MUG / KG I 
NSTG - COLBFG / NPRG .666 
HG - NSTG * GG * CPG * 3600! 



effectiveness 



MAL - SQR( (2! * HA) 


/ (K! 


* DELA 


/ 


12!) ) 


* (BA 


/ 


24!) 


MGL - SQR((2! * HG) 


/ (K! 


* DELG 


/ 


12!)) 


* (BG 


/ 


24!) 


ETAFA - TANH! (MAL) 


/ MAL 














ETAFG - TANH! (MGL) 


/ MGL 














SURFACE 


EFFECTIVENESS 












ETAOA - 1! - FRA * 


(1! - 


ETAFA) 












ETAOG - 1! - FRG * 


(1! - 


ETAFG) 













OVERALL COEFFICIENT OF HEAT TRANSFER 
RA - 1! / (ETAOA * HA) + 1! / ((AG / AA) * ETAOG * HG) +1! / (K! 

UA - 1! / RA 

EOLD - E 

★★★★ NTU and heat EXCHANGER EFFECTIVENESS 
CA - WANODE * CPA * 3600! 

CG - MDOT(J, I) * CPG * 3600! 

CMIN - MIN! (CA, CG) 

CR - MIN! (CMIN / CA, CMIN / CG) 

NTU - AA * UA / CMIN / NV% / NH% 



(Tp / 12!)) 
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IF CR - 1 THEN 

E -= NTU / (1 + NTU) 

ELSE 



X - 


EXP ( 


-NTU ^ 


(1! - CR)) 








E - 


(11 


- X) / 


(1! - CR * 


X) 






END IF 














TEXA - 


E 


(TIN(J, 


I) - TINA) 


♦ CMIN 


/ 


CA + TINA 


TEXG - 


E ^ 


(TINA - 


TIN(J, D) 


* CMIN 


/ 


CG + TIN(J 



check for convergence 
LOOP WHILE (ABS (EOLD / E - 1 ! ) > ERRLIM) 

★★★★bottom of e-ntu loop 

TBULKEG - TBULKEG + MOOT ( J, I) * TEXG 
TBULKEA - TBULPCEA + WANODE * TEXA 
NEXT I , J 

LOOP COMPLETE - GET BULK PROPERTIES FOR REST OF PROGRAM 
AND CALCULATE EFFECTIVENESS £ NTU 
TEXG - TBULKEG / WG 

TEXA - TBULKEA / WA 

TAVA - (TEXA + TINA) * .5 

TAVG - (TEXG + TING) * .5 



CALL 


TRANSP (TAVA, FARA, CPA, 


KA! , 


MUA, 


MA, 


FTYPA%) 


CALL 


TRANSP (TAVG, FARG, CPG, 


KG! , 


MUG, 


MG, 


FTYPG%) 


GG - 


WG / ACG 












NRA - 4! * RHA * GA / MUA 
NRG - 4! ♦ RHG * GG / MUG 












CALL 


INTERP (STATAO , NRA, 


COLBFA, 


FA) 






CALL 


INTERP ( STATG ( ) , NRG , 


COLBFG, 


FG) 







IF COLBFA - 01 OR COLBFG - 01 GOTO ERRl 
CA - WA * CPA * 3600! 

CG - WG * CPG 3600! 

CMIN - MIN! (CA, CG) 

CR - MIN! (CMIN / CA, CMIN / CG) 

' ★★★★★★ EFFECTIVENESS AND NTU FOR WHOLE HX 

E - (CA / CMIN) * (TEXA - TINA) / (TING - TINA) 

IF CR - 1 THEN 

NTU - E / (1 - E) 

ELSE 

NTU - LOG((l! - E) / (1! - CR ^ E) ) / (CR - 1!) 

END IF 

/ ★★★★★ pressure drop FOR WHOLE HX CALCULATED HERE’^*^ 

/★★★★★★★★★★★★★★★ inlet and exit loss COEFFICIENTS ★★★★★★★★★★★★★★★★★★★★ 

/ 

' CCA/G JET CONTRACTION -AREA RATIO 

CCA - .61000000001# - .14442945071# ^ SIGA + 1.0080347435# * SIGA ^ 2 
CCA = CCA - 1.7317560083# * SIGA ^ 3 + 1.1559407939# ^ SIGA ^ 4 
CCG - .61000000001# - .14442945071# * SIGG + 1.0080347435# * SIGG ^ 2 
CCG * CCG - 1.7317560083# * SIGG ^ 3 + 1.1559407939# * SIGG ^ 4 
NARA - NRA * .0001 



■ ■ .jOK 
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NARG - NRG * .0001 

KDA! - 1.1063960104# - .13322445533# ♦ NARA + .11885428625# ^ NARA ^ 2 
KDAI - KDAI - .033170530592# ^ NARA ^ 3 

KDGI - 1.1063960104# - .13322445533# * NARG + .11885428625# * NARG ^ 2 
KDGI - KDGI - .033170530592# * NARG ^ 3 

KCAI - (II - 21 ^ CCA + CCA 2 * (21 * KDAI - II)) / CCA ^ 2 

KCG! - (II - 2! ^ CCG + CCG ^ 2 * (21 * KDGI - II)) / CCG ^ 2 

KEAI - II - 21 * KDAI ^ SIGA + SIGA ^ 2 

KEG I -= II - 21 ♦ KDGI * SIGG + SIGG ^ 2 

/*** A A A A ****'**** PRESSURE DROPS ★***★*★***★★*****♦*★★♦*★♦★★★★★ 

t 

RHOEXA - PEXA / TEXA / RU * 1441 

RHOING - PING / TING / RU ♦ 144! 

RHOEXG - PEXG / TEXG / RU * 1441 

VEXA - .636 ♦ VINA * SQR(RHOINA / RHOEXA) 

DEXA - SQR(4I * WA / RHOEXA / VEXA /PI) 

HINA - RHOINA / 21 / GC * VINA ^ 21 
DELPAH - .595 ★ HINA / PINA / 1441 



ILAI - 


11 - 


SIGA ^ 


21 


+ KCAI 


ILGI - 


11 - 


SIGG ^ 


21 


+ KCG I 


SPVA - 


PINA 


/ PEXA 


★ 


TEXA / TINA 


SPVG - 


PING 


/ PEXG 


■k 


TEXG / TING 


ELA - 


(11 - 


SIGA ^ 


21 


- KEA I ) * SPVA 


ELG - 


(11 - 


SIGG 


21 


- KEG I ) ^ SPVG 



SPVAM - 21 ^ PINA / (PINA + PEXA) ^ TAVA / TINA 

SPVGM - 21 * PING / (PING -f PEXG) * TAVG / TING 

CFA - FA ^ AA / ACA ♦ SPVAM 
CFG - FG * AG / ACG * SPVGM 
FAA - 21 ^ (SPVA -II) 

FAG - 21 * (SPVG - II) 

TLA - IIAI + FAA + CFA - EIA 
TLG -= ILGI + FAG + CFG - ELG 
LHGING - SQR(XNCFL 2 1 - DEXA ^ 21) 

LHGEXG - SQR (XNCFL ^21- DINA 21) 

ANGEXG - ATN(DINA / LHGEXG) 

ANGING - ATN (DEXA / LHGING) 

ANGINA * PI / 21 - ANGEXG 
ANGEXA - PI / 2! - ANGING 
HFXG - II + ((DINA + DEXA) / 21 / FLA) 

HFXA - (LHGING / 21 + FLA + LHGEXG / 21) / FLA 



***** CORE PRESSURE LOSS 



DELPAC - (GA 


/ 


144! 


/ 


PINA) 


^ 21 


/ 


21 


/ 


GC ♦ 


15451 


/ 


MA * 


TINA * 


TLA ^ 


HFXA 


DELPGC - (GG 


/ 


1441 


/ 


PING) 


^ 21 


/ 


21 


/ 


GC * 


15451 


/ 


MG * 


TING ♦ 


TLG ^ 


HFXG 



AHCMINA - SIGA ^ DINA ^ XNCFL 
AHCMING - SIGG * XNCFL ♦ LHGING 
VINAHl - WA / RHOINA / AHCMINA 
VINGHl - WG / RHOING / AHCMING 
VINACl « VINAHl ^ COS (ANGINA) 
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21 ) / 21 ) 

21 ) / 21 ) 

♦ VINAM ^ 2! 

* VINGM ^ 21 



VINGCl - VINGHl * COS (ANGING) 

CALL BENDLOS (ANGINA, KBINA!) 

CALL BENDLOS (ANGING, KBING!) 

VINAM - SQR((VINAH1 ^ 21 + VINACl 
VINGM « SQR( (VINGHl ^ 2! + VINGCl ^ 

DELPABl - RHOINA * KBINA! / 21 / GC 

DELPGBl - RHOING ^ KBING! / 2! / GC 

AHCMEXA - SIGA * DEXA * XNCFL 
AHCMEXG - SIGG ^ LHGEXG ♦ XNCFL 
VEXAH2 - WA / RHOEXA / AHCMEXA 
VEXGH2 - WG / RHOEXG / AHCMEXG 
VEXAC2 - VEXAH2 * COS (ANGEXA) 

VEXGC2 - VEXGH2 ♦ COS (ANGEXG) 

CALL BENDLOS (ANGEXA, KBEXA!) 

CALL BENDLOS (ANGEXG, KBEXG 1 ) 

VEXAM - SQR((VEXAH2 ^ 2 + VEXAC2 
VEXGM - SQR ( (VEXGH2 ^ 2 + VEXGC2 
DELPAB2 - RHOEXA ♦ KBEXA 1 / 21 / 

DELPGB2 - RHOEXG ^ KBEXG! / 2! / 

DELPAB - (DELPABl + DELPAB2) / PINA / 1441 

DELPGB - (DELPGBl + DELPGB2) / PING / 144! 

DELPA - DELPAC + DELPAH + DELPAB 
DELPG - DELPGC + DELPGB 
PEXA - PINA * (1! - DELPA) 

PING - PEXG ♦ (1! + DELPG) 

PCDELPA - 1001 ♦ DELPA 
PCDELPG - 1001 ♦ DELPG 
DELPT - DELPA + DELPG 
PCDELPT - 100! * DELPT 
PCNE « 100! ♦ E 

LOOP WHILE (ABS(DPT / DELPT - II) > ERRLIM) : ' 



^ 2) 


/ 


21) 




A 2) 


/ 


2!) 




GC * 


VEXAM ^ 


2 1 


GC ^ 


VEXGM ^ 


2! 



CHECK FOR CONVERGENCE 



^ BOTTOM OF PRESSURE LOOP*^*** 

' **<***«^****ik**« WEIGHT CALCULATIONS OF THE HEAT EXCHANGER ************ 



WTA -= 


AA ^ 


RHO ^ 


(FRA * DELA * 


WFA + 


Tp * 


(1! - FRA) 


♦ WPA) 


/ 


24! 


WTG - 


AG ♦ 


RHO ♦ 


(FRG * DELG ♦ 


WFG + 


Tp * 


(11 - FRG) 


* WPG) 


/ 


241 


WPLA ' 


- 151 


♦ (41 


♦ XNCFL * FLA 


+ PI 


* XNCFL ♦ (DINA 


+ DEXA) 


+ 


LHGING * DEXA + LHGEXG * 



DINA) 

WIE » (WTA + WTG) / VOL * (LHGING * DEXA * XNCFL + LHGEXG * DINA * XNCFL) / 2! 
WHXT •= WTA + WTG + WPLA + WIE 
OVALHT - DINA + DEXA + FLA 

WALL TEMPERATURE CALCULATIONS AT CORE EXIT 
' ★★★ TEMPERATURES ARE CALCULATED USING BULK AVERAGE ENTRANCE 

' *** AND EXIT TEMPERATURES, NOT PEAK CHANNEL TEMPERATURES *** 

r 

TA(1) - TEXA 
TA(2) - TINA 
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TG<1) -= TING 
TG(2) - TEXG 

FOR J - 1 TO 2 

CALL TRANSP (TA(J) , FARA, CPA, KA I , MUA, MA, FTYPA%) 

CALL TRANSP (TG(J) , FARG, CPG, KG!, MUG, MG, FTYPG%) 

XNRA(J) - 41 * RHA * GA / MUA 

XNRG(J) « 4! * RHG * GG / MUG 

CALL INTERP (STATA 0 , XNRA < J) , COLBFA, FA) 

IF COLBFA - 0! GOTO ERRl 
NPRA - CPA * MUA / KA! 

NSTA « COLBFA / NPRA ^ . 666 

HA - NSTA * GA CPA * 3 6001 

CALL INTERP (STATG 0 , XNRG ( J) , COLBFG, FG) 

IF COLBFG - 0! GOTO ERRl 
NPRG - CPG * MUG / KG! 

NSTG - COLBFG / NPRG ^ .666 

HG * NSTG GG * CPG * 3 6001 

MAL - SQR((2! HA) / (K! * DELA / 12!)) * <BA / 24 ! ) 

MGL - SQR((2! * HG) / (K! * DELG / 12!)) * (BG / 24!) 

ETAFA - TANH! (MAL) / MAL 

ETAFG - TANH! (MGL) / MGL 

ETAOA - 1! - FRA * (1! - ETAFA) 

ETAOG - 1! - FRG * (1! - ETAFG) 

RA - 1! / (ETAOA * HA) + 1! / (AG / AA * ETAOG * HG) 

UA - 1! / RA 

TW(J) - TG(J) - UA / (AG / AA * ETAOG * HG) * (TG ( J) 
TGRW(J) - TG(J) - TW(J) 

NEXT J 



+ 1! / (K! / (Tp / 12!) ) 

- TA(J)) 



★ ★★★★ OUTPUT SECTION*’^*** 



XCDELPA - PINA . 01 * PCDELPA * 1728! / 62.4 
XCDELPG - PING * . 01 * PCDELPG 1728! / 62.4 
OPEN ”SCRN:** FOR OUTPUT AS #1 
CLS 

I - 1 : GOSUB OUTPRINT : 

GOSUB HOLD 

RUN$ - "OUTPUT.” + RIGHT$ (STR$ (IRUN) , LEN (STR$ (IRUN) ) - 1) 

CLS : LOCATE 8, 1 

GOSUB CHECKOUT 
LOCATE 10, 1 

PRINT ’’RESULTS FROM RUN NUMBER IRUN; ” HAVE BEEN PLACED IN FILE RUN$ 
OPEN RUN$ FOR OUTPUT AS #2 
1*2: GOSUB OUTPRINT: 

CLOSE 2 

PRINT : PRINT ”DO YOU WANT A HARDCOPY? ’’; 

IF YN% * YES THEN 

PRINT ’’SET UP PRINTER." 

PRINT "PRESS <RETURN> WHEN READY.” 
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WHILE INKEY$ - WEND 

OPEN ”LPT1:” FOR OUTPUT AS #3 
1-3: GOSUB OUTPRINT: 

CLOSE 3 
END IF 

ERRl: 

PRINT ”DO YOU WANT TO SAVE THIS RUN AS AN INPUT FILE? 

IF YN% - YES THEN 
RUN$ - DOC$ 

GOSUB CHECKOUT 
DOC$ - RUN$ 

OPEN DOC$ FOR OUTPUT AS #2 
PRINT 

PRINT "SAVING RUN # IRUN; ” AS INPUT FILE DOC$ : PRINT 



PRINT 


#2, 


COMMENT$ 








WRITE 


#2, 


UT%; ' 




type of units 




WRITE 


#2, 


TYPA%, NSA%, 


FTYPA% : ' 


cold side surface 




WRITE 


#2, 


TYPG%, NSG%, 


FTYPG%: ' 


hot side surface 




WRITE 


#2, 


RHO / C(UT%, 


7) * C (UT%, 4) , K1 


/ C (UT%, 8) , Tp / C (UT%, 


1) 


WRITE 


#2, 


AFRA / C (UT% 


, 3) , RLENI / C (UT%, 


2) 




WRITE 


#2, 


WA / C(UT%, 


7) , TINA / C (UT%, 5) 


, PINA / C(UT%, 6), VINA 


/ C(UT% 


WRITE 


#2, 


WG / C(UT%, 


7), PEXG / C(UT%, 6) 






WRITE 


#2, 


NH%, NV% 








FOR J 


- 1 


TO NV% 








FOR I 


- 1 


TO NH% 









WRITE #2, VING(J, I) / C(UT%, 2), TIN ( J, I) / C(UT%, 5) 

NEXT I, J 

WRITE #2, FARA, FARG 
CLOSE 2 
END IF 

PRINT : PRINT "DO YOU WISH TO MAKE ANOTHER RUN? 

CLOSE 

IF YN% - YES THEN IRUN - IRUN + 1: GOTO TOP: 

END 

GASVELPRO : 

' GAS VELOCITY PROFILE WILL BE ENTERED 

PRINT "HOT SIDE VELOCITY DISTRIBUTION ” 

DO: INPUT "HOW MANY HORIZONTAL NODES? "; NH%: LOOP WHILE NH% < 1 
DO: INPUT "HOW MANY VERTICAL NODES? NV% : LOOP WHILE NV% < 1 

REDIM VING(NV%, NH%) , TIN(NV%, NH%), MDOT (NV%, NH%) 

GASVELPR02 : 

IF NH% + NV% > 2 THEN 

PRINT " IS THE HOT SIDE INLET TEMPERATURE THE SAME FOR ALL NODES? "; 
II - YN% 

IF II - YES THEN 

PRINT "ENTER THE HOT SIDE INLET TEMPERATURE (”; C$ (UT%, 5); ”) 

INPUT TING 

TING - TING * C(UT%, 5) 
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PRINT “ENTER THE VELOCITY (“; C$(UT%, 2); “/sec) FOR EACH NODE . “ 

PRINT 

ELSE 

PRINT “ENTER THE VELOCITY (“; C$ (UT%, 2); "/SEC) AND TEMPERATURE ("; C$(UT%, 5) ; 

« j « 

PRINT "FOR EACH NODE POINT SEPARATED BY COMMAS . " : PRINT 
END IF 

FOR J - 1 TO NV% 

FOR I - 1 TO NH% 

NODE: PRINT "NODE J; I; ") - 

IF II - YES THEN 

TIN(J, I) * TING 
INPUT VING(J, I) 

VING(J, I) * VING(J, I) ♦ C(UT%, 2) 

ELSE 

INPUT VING(J, I), TIN(J, I) 

IF TIN(J, I) - 0! THEN BEEP: GOTO NODE 
VING(J, I) - VING(J, I) * C(UT%, 2) 

TIN(J, I) - TIN(J, I) ^ C(UT%, 5) 

END IF 

MDOT(J, I) - VING(J, I) / TIN(J, I) 

XMAX - XMAX + MDOT(J, I) 

NEXT I , J 

' NORMALIZE THE NODAL MASS FLOW RATES & FIND AVE INLET TEMP 
TING - 01 
FOR J * 1 TO NV% 

FOR I - 1 TO NH% 

MDOT(J, I) - MDOT(J, I) ♦ WG / XMAX 
TING - TING + TIN(J, I) ♦ MOOT { J, I) 

NEXT I , J 
TING - TING / WG 
ELSE 

PRINT "ENTER THE HOT SIDE INLET TEMPERATURE C$ (UT%, 5); ") "; 

INPUT TING 

TING - TING * C(UT%, 5) 

PRINT "ENTER THE VELOCITY ("; C$ (UT%, 2); "/sec) "; 

INPUT VING (1, 1) 

MDOT (1, 1) - WG 
PRINT 
END IF 
RETURN 

/ ♦★★★★★★★★♦★★★★★★★★★★OUTPUT SUBROUTINE'*^***^***^^^*^^^*’'^****** 



OUTPRINT: 

PRINT #1, "COMPACT HEAT EXCHANGER ANALYSIS": PRINT #1, 
PRINT #1, "RUN #"; IRUN; " INPUT FILE "/ DOC$ 

PRINT #1, NH%; "x"; NV%; "NODES" 

PRINT #1, COMMENT$ 
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PRINT #1, 

PRINT #1, ”CORE HEAT TRANSFER SURFACE"; 

PRINT #1, TAB(35); "COLD-SIDE"; TAB(50); "HOT-SIDE" 

PRINT #1, 

PRINT #1, "FIN TYPE"; TAB(35); SURFACE$ (TYPA%) ; TAB(50); SURFACE$ ( TYPG% ) 

PRINT #1, "FIN AND SURFACE NUMBER"; TAB(35); 

PRINT #1, USING "#&## "; TYPA%; S; NSA%; TYPG%; S; NSG% 

PRINT #1, "Fins/"; C$ (UT%, 1); TAB (35); SFA * C (UT%, 1); 

PRINT #1, TAB(50); SFG * C(UT%, 1> 

IF AXA <> 0 OR AXG <> 0 THEN 

PRINT #1, "FIN LENGTH ("; C$ (UT%, 1); ")"; TAB (35); AXA * C(UT%, 1) ; 

PRINT #1, TAB (50); AXG * C(UT%, 1) 

END IF 

PRINT #1, "PLATE SPACING ("; C$ (UT%, 1); ")"; TAB (35); BA / C (UT%, 1); 

PRINT #1, TAB(50); BG / C (UT%, 1) 

PRINT #1, "HYDRAULIC RADIUS ("; C$ (UT%, 1); ")"; TAB(35); 

PRINT #1, RHA * 12 / C (UT%, 1); TAB (50); RHG 12 / C(UT%, 1) 

PRINT #1, "FIN THICKNESS ("; C$(UT%, 1); ")"; TAB(35); DELA / C(UT%, 1); 

PRINT #1, TAB (50); DELG / C(UT%, 1) 

PRINT #1, "COMPACTNESS ("; C$ (UT%, 3); "/"; C$ (UT%, 4); ")"; TAB(35); 

PRINT #1, BETA * C (UT%, 2); TAB (50); BETG ♦ C (UT%, 2) 

PRINT #1, "FIN/TOTAL AREA ("; C$ (UT%, 3); "/"; C$ (UT%, 3); ")"; TAB (35); 

PRINT #1, FRA; TAB (50); FRG 

PRINT #1, "FLUID TYPE"; TAB(35); FTYPE$ (FTYPA%) ; TAB(50); FTYPE$ (FTYPG% ) 

PRINT #1, 

PRINT #1, "LENGTH - "; TAB (19); RLENI / C(UT%, 2); C$ (UT%, 2); TAB (39); 

PRINT #1, "VOLUME - "; TAB (48); VOL / C(UT%, 4); C$ (UT%, 4) 

PRINT #1, "FRONTAL AREA - "; TAB(19); AFRA / C(UT%, 3); C$ (UT%, 3); TAB (39); 

PRINT #1, "HEIGHT * "; TAB (48); OVALHT / C (UT%, 2); C$(UT%, 2) 

PRINT #1, "METAL DENSITY - "; TAB (19); RHO / C(UT%, 7) * C(UT%, 4); C$ (UT%, 7); "/"; 

C$(UT%, 4); 

PRINT #1, TAB (39); "WEIGHT -"; TAB (48); WHXT / C (UT%, 7); C$(UT%, 7); "s" 

PRINT #1, "PLATE THICKNESS TAB(19); Tp / C(UT%, 1); C$ (UT%, 1); 

PRINT #1, TAB(39); "K -"; TAB(48); K! / C(UT%, 8); C$(UT%, 8) 

IF I - 1 THEN GOSUB HOLD ELSE PRINT #1, STRING$ (70, "-") 

PRINT #1, "COLD HEADER DESIGN DETAILS" 

PRINT #1, TAB (35); "INLET"; TAB (61); "EXIT" 

PRINT #1, "DIAMETER ("; C$(UT%, 2); ")"; TAB (29); DINA / C (UT%, 2) ; 

PRINT #1, TAB (48); DEXA / C(UT%, 2) 

PRINT #1, "VELOCITY ("; C$ (UT%, 2); "/s©c)"; TAB (29); VINA / C (UT%, 2); 

PRINT #1, TAB (48); VEXA / C (UT%, 2) 

PRINT #1, STRING$(70, "-") 

PRINT #1, "HEAT EXCHANGER CONDITIONS" 

PRINT #1, TAB(35); "COLD-SIDE"; TAB(57); "HOT-SIDE" 

PRINT #1, TAB (30); "INLET"; TAB (43); "EXIT"; TAB (53); "INLET"; TAB (65); "EXIT" 

PRINT #1, "PRESSURE ("; C$ (UT%, 6); ")"; TAB (28); PINA / C (UT%, 6); TAB (39); 

PRINT #1, PEXA / C(UT%, 6); TAB (50); PING / C (UT%, 6); TAB (61); PEXG / C (UT%, 6) 

PRINT #1, "TEMPERATURE ("; C$ (UT%, 5); ")"; TAB (28); TINA / C(UT%, 5); TAB (39); 

PRINT #1, TEXA / C(UT%, 5); TAB (50); TING / C(UT%, 5); TAB (61); TEXG / C (UT%, 5) 
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PRINT #1, 
PRINT 
IF FTYPA% 
PRINT 
END IF 
PRINT #1, 
PRINT #1, 
PRINT #1, 
PRINT 
PRINT #1, 
PRINT #1, 
PRINT 
PRINT #1, 
PRINT #1, 
PRINT #1, 
PRINT 
PRINT #1, 
RETURN 

riNTYPE : 
PRINT 



"MASSFLOW (•’; C$(UT%, 7); "/see)”; TAB (33); WA / C(UT%, 7) ; 

#1, TAB(55); WG / C(UT%, 7) 

- 1 OR FTYPG% - 1 THEN 

#1, "FUEL-AIR RATIO"; TAB(33); FARA; TAB(55); FARG 
"PRESSURE DROP (%)"; TAB(33); PCDELPA; TAB(55); PCDELPG 

"PRESSURE DROP (in H20)"; TAB (33); XCDELPA; TAB (55); XCDELPG; PRINT #1, 
"HOT SIDE WALL TEMPERATURE ("; C$(UT%, 5); ")”; 

#1, TAB(50); TW(1) / C (UT%, 5); TAB(61); TW(2) / C(UT%, 5) 

"REYNOLDS NUMBER"; TAB(50); XNRG(l); TAB(6l); XNRG(2) 

”Tw-Tg ("; C$(UT%, 5); ")"; TAB (50); 

#1, TGRW(l) / C(UT%, 5); TAB(61); TGRW(2) / C(UT%, 5) 

STRING$(70, "-") 



"EFFECTIVENESS - "; TAB (19); PCNE; "%"; 
#1, "NTU - "; TAB(48); NTU r PRINT #1, 
"TOTAL PRESSURE DROP « "/ PCDELPT; ”%" 



TAB(39) . 



PRINT 


” TYPE 


1 - PLAIN FIN 


SURFACE 


NUMBER - 


1 


TO 


18 


PRINT 


" TYPE 


2 - LOUVERED FIN 


SURFACE 


NUMBER - 


1 


TO 


14 


PRINT 


" TYPE 


3 - STRIP-OFFSET FIN 


SURFACE 


NUMBER - 


1 


TO 


18 


PRINT 


It type 


4 - WAVY FIN 


SURFACE 


NUMBER * 


1 


TO 


3 



PRINT 

RETURN 

HOLD: 



*★★★★ Suspends program execution until a )cey is pressed 



LOCATE 25, 10 

PRINT "PRESS ANY KEY TO CONTINUE"; 

DEF SEG - 0: POKE 1050, PEEK (1052) : DEF SEG 
DO: LOOP WHILE INKEY$ - "” 

CLS 

RETURN 



CHECKOUT : 

' CHECKS FOR A BAD OUTPUT FILE OR WARNS OF OVERWRITE 

PRTNT "WHAT IS THE NAME OF THE FILE <"; RUN$; ">"; 

INPUT F$ 

IF F$ - "” THEN F$ - RUN$ 

EFLAG - FALSE 
ON ERROR GOTO FILEERR2 
OPEN F$ FOR INPUT AS #2 
CLOSE 2 

ON ERROR GOTO 0 
IF NOT EFLAG THEN 

PRINT "WARNING - FILE "; F$; ” EXISTS. OVERWRITE?"; 

IF YN% - NO GOTO CHECKOUT 
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END IF 
RUN$ = F$ 

RETURN 

rilXERRl: 

IF ERR - 53 THEN 

' get another file name 

PRINT "File UCASE$ (DOC$ ) ; " not found." 

INPUT "FILE NAME"; DOC$ 

RESUME 

ELSE 

' some other error, so print message and abort 
PRINT "Unrecoverable error--"; ERR 
ON ERROR GOTO 0 
END IF 

riLEEKR2 : 

' BAD OUTPUT FILE TRAP 

IF ERR - 64 THEN 

PRINT "BAD FILE NAME - RETYPE" 

ELSEIF ERR - 53 THEN 
EFLAG - TRUE 
RESUME NEXT 
ELSE 

BEEP 
END IF 

RESUME CHECKOUT: 

BADSURTNUM: 

' bad SURFACE NUMBER TRAP 

CLS 

PRINT "BAD SURFACE OR TYPE NUMBER” 

GOSUB FINTYPE: 

EFLAG « TRUE 
RESUME NEXT 

END 

SUB BENDLOS (X, Y) STATIC 

f ★★★★★ subroutine to calculate the pressure loss at a bend 

2 - X * 57.29578 

Y - 2.922713E-02 - 2.639695E-03 ^ 2 + 2.272872E-04 ^ 2 2 - 1.850293E-06 ^2^3+ 

3.655184E-08 ^2^4- 4.49784E-10 ^2^5+ 2.088911E-12 ^2^6 
END SUB 

SUB INTERP (A(), RE, COL, F) STATIC 
J - 1 

WHILE RE < A(J, 1) : J » J + 1: WEND 
I - J - 1 
IF J » 1 THEN 

PRINT " REYNOLDS NUMBER OUT OF RANGE OF PROGRAMMED TABLES - " ; RE 
COL * 01 





|i 
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F «= 0 1 
EXIT SUB 
ELSE 

IF A(J, 1) - 01 THEN 
F - A(2, 4) / RE 

COL - A (1, 4) / RE ^ .7 
ELSE 

Z - (A(I, 1) - RE) / (A(I, 1) - A(J, D) 

COL - Z ♦ (A(J, 2) - A(I, 2) ) 

COL - A (I, 2) + COL 

F - Z * (A(J, 3) - A(I, 3)) 

F - A(I, 3) + F 
END IF 
END IF 
END SUB 

FUNCTION MINI (Xl, YI ) STATIC 

IF XI > YI THEN MINI - Yl ELSE MINI - XI 
END FUNCTION 

SUB STAT (TTPE$, nn%, A()) STATIC 

t ★★★★★★★*★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★ 

' SUBROUTINE STAT RETURNS STANTON NUMBERS AND FRICTION FLOW DATA 

' FOR THE TYPE HEAT EXCHANGER SPECIFIED 

OPEN TYPE$ FOR RANDOM AS #10 LEN - 16 
FIELD #10, 16 AS Z$ 

' read data from data base 

I - (nn% - 1) * 18 
FOR J - 1 TO 18 

GET #10, J + I 
FOR K - 1 TO 4 

A(J, K) - CVS (MID$(Z$, (K “ 1) *4 + 1, 4)) 

NEXT K, J 
CLOSE 10 
END SUB 

SUB SURF (T%, NS%, SIDE$, AA, BB, SF, PS, RH, DEL, BET, FR, WF, WP) STATIC 

/ ***** Subroutine to return heat exchanger surface properties from the 

' data base files. 

SHARED C(), C$0, SURFACE$(), UT%, EFLAG 

DO: CLS 

PRINT "THE CURRENT SIDE$; " SIDE HEAT EXCHANGER SURFACE IS AS FOLLOWS:” 
ON ERROR GOTO BADSURFNUM: 

EFLAG - FALSE 

OPEN SURFACE$(T%) + ”.DAT” FOR INPUT AS #10 
IF T% < 1 OR NS% < 1 THEN EFLAG - TRUE 
FOR J - 1 TO NS% 

INPUT #10, A, AA, BB, SF, PS, RH, DEL, B, BET, FR, WF, WP 
IF EFLAG THEN CLOSE 10: EXIT FOR 
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NEXT J 

ON ERROR GOTO 0 
IF EFIAG THEN 

I - NO: ' CATCHES A BAD SURFACE OR TYPE NUMBER 

ELSE 

RH - RH / 4! 

CLOSE 10 

PRINT : PRINT , SURFACE$ (T%) ; ” FINS” 

PRINT 

PRINT , "FIN TYPE - T% 

PRINT , "SURFACE NUMBER - "/ J - 1 
PRINT 

PRINT "PLATE SPACING”, , PS / C(UT%, 1); C$ (UT%, 1) 

PRINT "FIN THICKNESS", , DEL / C(UT%, 1); C$(UT%, 1) 

PRINT "FINS/"; C$(UT%, 1), , SF * C(UT%, 1) 

IF AA > 0 THEN PRINT "FIN LENGTH”, , AA / C(UT%, 1); C$(UT%, 1) 

PRINT "HYDRAULIC RADIUS”, RH * 12 / C(UT%, 1); C$(UT%, 1) 

PRINT "COMPACTNESS", , BET * C(UT%, 2); C$ (UT%, 3); "/"; C$(UT%, 4) 

PRINT "FIN/TOTAL AREA”, FR; C$(UT%, 3); "/”; C$ (UT%, 3) 

PRINT 

PRINT "IS THE ABOVE CORRECT? : I - YN% 

END IF 

IF I « NO THEN 

PRINT : PRINT "ENTER THE "; SIDE$; " SIDE HEAT EXCHANGER SURFACE" 

PRINT 

INPUT "FIN TYPE"; T% 

INPUT "SURFACE NUMBER"; NS% 

END IF 

LOOP WHILE I - NO 
END SUB 

FUNCTION TANH! (X) STATIC 

TANH! - (EXP (X) - EXP(-X)) / (EXP (X) + EXP(-X)) 

END FUNCTION 

SUB TRANSP (T, FAR, CP, TK, MU, MW, GTYPE%) STATIC 

f ★★★*★ SUBROUTINE TO GET TRANSPORT PROPERTIES AT SPECIFIED TEMP 

SELECT CASE GTYPE% 

CASE 1 : ' AIR-FUEL MIXTURE TRANSPORT PROPERTIES 

IF T < 5001 OR T > 2000! GOTO OUTOFRANGE; 

IF FAR < 01 OR FAR > .034826 GOTO OUTOFRANGE; 

CP - ( ( ( ( ( (A1 * T + A2) * T + A3) * T + A4) * T + A5) * T + A6) * T + A7 ) * T + A8 

IF FAR <> 0! THEN 

CPF - ( ( ( ( ( (B1 * T + B2) * T + B3) * T + B4) * T + B5) * T + B6) ♦ T + B7) T + B8 
CP - (CP + FAR * CPF) / (1! + FAR) 

END IF 

TK - ( ( ( ( (Cl T + C2) * T + C3) * T + C4) * T + C5) * T + C6) * T + C7 
IF FAR <> 0! THEN 
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TKF - ( ( ( ( (D1 * T + D2) * T + D3) * T + D4) * T + D5) * T + D6) * T + D7 
TK - (TK + FAR * TKF) / (1! + FAR) 

END IF 

TK - TK / 3600! 

MU - (((((El * T + E2) * T + E3) * T + E4) * T + E5) ♦ T + E6) * T + E7 
IF FAR <> 0! THEN 

MUF - (((((FI * T + F2) ^ T + F3) ♦ T + F4) ^ T + F5) * T + F6) * T + F7 
MU - (MU + FAR * MUF) / (1! + FAR) 

END IF 

MU - MU / 36001 

MW - 28.97 - .946186 ♦ FAR 

CASE 2 : ' HELIUM TRANSPORT PROPERTIES 



CP - 1.2404 

MU - .0006388 ★ T ^ .687 / 3600! 

TK - .001062 ♦ T .687 / 3600! 

MW - 4.002602 

CASE ELSE 

PRINT "INCORRECT FLUID NUMBER", GTYPE% 

END 

END SELECT 
EXIT SUB 

OUTOFRANGE : 

PRINT SPC(IO); "TRANSP INPUT OUT OF RANGE: TEMP - T; " FAR - FAR 

STOP 
END SUB 

FUNCTION YN% STATIC 

ROW - CSRLIN: COL - POS(O) 

DO 

DEF SEG - 0: POKE 1050, PEEK (1052) ; DEF SEG : ' clear )coyboard buffer 

X$ - ANS$ - 

DO: X$ « INKEY$: LOOP WHILE X$ « "" 



IF 


X$ - "Y" 


OR 


X$ - "y" 


THEN 


ANS$ « "Y" ! 


! YN% - YES 


IF 


X$ - "N" 


OR 


X$ - "n" 


THEN 


ANS$ - "N" : 


: YN% - FALSE 



IF ANS$ - "Y” OR ANS$ - "N" THEN 
LOCATE 25, 10; PRINT SPC(50); 

LOCATE ROW, COL: PRINT ANS$ 

ELSE 

BEEP 

LOCATE 25, 10: PRINT "PLEASE ANSWER YES OR NO"; : ANS$ - 
END IF 

LOOP UNTIL ANS$ « "Y" OR ANS$ • ”N” 



END FUNCTION 
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B.2 Precooler Design Program PRECOOL. BAS. 



PRECOOL.BAS uses the design method, as opposed to the analysis method used in 
COMPHX.BAS, to size the cross-flow precoolers. It uses many of the same steps as outlined 
above in COMPHX.BAS however in this case effectiveness and pressure drop are given and 
size is iterated until the calculated effectiveness and pressure drop match the desired values. 

PRECOOL.BAS sizes the cross-flow heat exchangers used for the precoolers. It was 
taken, with very few modifications, directly from Appendix C in Reference 8. The program 
is line for line the same as Staudt’s program except I substituted an alternate effectiveness 
equation (Equation B.l instead of B.2) which allowed me to model lower water mass flow 
rates than were allowed in Staudt’s program. 1 made water mass flow rate an input w'hich 
allowed me to calculate the capacity ratio (Equation B.IO) and NTU (Equation B.ll) instead 
of inputting them. Finally, I included tabulated values for the fluid properties of water (vis- 
cosity, thermal conductivity, and Prandtl number) so the user does not have to perfonn the 
interpolation manually.* 

The equation used for NTU is: 



For a more detailed description of the program theory see chapter five of Reference 8. 
A sample of the output and the program listing follow. 



C, 



B.13. 




B.14. 



Water properties are from Pitts & Sissom. 
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B.2.1 Sample Output from PRECOOL. BAS 



★ ★★★★★★★★★★★★★★★★★★★★★★★★★★★Hr*********** 

PRECOOLER DESIGN PROGRAM 



HELIUM TEMPS: INLET = 167 *C OUTLET 

HELIUM INLET PRESSURE = 4.01 MPa 
HELIUM MASS FLOW =27.4 KG/SEC 
WATER: INLET TEMP =20 ’C OUTLET TEMP 

MASS FLOW = 250 KG/SEC 3970.525 GPM 
HEAT EXCHANGER SURFACE: S 1.50--125(s) 

RATIO OF HEAT CAPACITIES, (CW/Ch) - 7.360062 



30 ’C 



38.61398 *C 



CROSSFLOW HEAT EXCHANGER 
★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★ 

RESULTS OF CALCULATIONS 

CROSSFLOW HEAT EXCHANGER DIMENSIONS 
VOL = 1 . 077307 m3 
TUBESHEET 

# OF TUBES = 11901 TUBE LENGTH = 1.193399 m 

WIDTH = 1.1811 m DEPTH = .7643056 m 

GAS FLOW 

PRESSURE DROP = 2.601114E-03 FRONTAL AREA = 1.409524 

GEOMETRY RATIO = .6437702 

EFFECTIVENESS = .9319728 NTU - 3.344671 
★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★ 
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B.2.2 Source Code Listing for FRECOOL.BAS 

DECLARE FUNCTION INTERP ! (Tl, Al()) 

CONST PIE - 3.14159, E « 2.71828 
DIM WPR(5) , WK (5) 

FOR I - 0 TO 5 

READ WPR(I) , WK(I) 

WK<I) - WK(I) * 1.729577; ' convert to metric units 

NEXT I 

' TABULATED DATA FOR WATER PRANDEL NUMBER AND THERMAL CONDUCTIVITY 
' Pr k 



DATA 


13.6, 


.319 


DATA 


7.02, 


.345 


DATA 


4.34, 


.363 


DATA 


3.02, 


.376 


DATA 


2.22, 


.386 


DATA 


1.74, 


.393 



PRINT "Preoooler design prograim*' 

***nrw*HEAT EXCHANGER DATA IS PROVIDED IN THE DATA STATEMENT WHICH FOLLOW 
THE DATA NEEDED AN THE UNITS REQUIRED ARE LISTED BELOW 



PARAMETER 



UNITS 



VARIABLE USED 



OUTSIDE DIAMETER 


m 


OD 


EFFECTIVE FIN LENGTH 


m 


DF 


FIN PITCH FiNS/m 


FP 


HYDRAULIC DIAMETER 


m 


RH 


FIN THICKNESS 


m 


DEL 


FREE FLOW/ FRONTAL AREA 


m^2/m^2 


SIGMA 


FIN AREA/TOTAL AREA 


m"‘2/m^2 


AFA 


HEAT TRANSFER AREA/VOLUME 


m''2/m^3 


ALFA 


LONGITUDINAL TUBE SPACING 


m 


DT 


TRANSVERSE TUBE SPACING 


m 


WT 


Jl, J2, FI, F2 






SURFACE NAME 




S$ 


READ OD, LF, FP, RH, DEL, 


SIGMA, AFA, 


ALFA, 



★★★★data for surface CF-8.72 

DATA .01067, .0056,343, .004425, .00048, . 494, . 876, 446, .02032, .024765 
DATA .222, .402, .29, , 246, ”CF-8 . 72 

★★★★data for surface S 1.50-1.00 
DATA .009525,0,0, .006071,0, .333,0,220.1, .009525, .014288 
DATA .299, .399, .384, .221, "S 1.50-1.00” 

★★★★data for surface S 1.50-1.25 (s) 

DATA .00635,0,0, .005029,0, .333,0,263.5, .007938, .009525 
DATA .3460, .4153, .2999, .1842, ”S 1.50-125 (s)" 

★★★★data for surface 8.0 3/8T 

DATA .01021, .007595,314.96, .00363, .00033, .534, .913,587., .0220, .0254 
DATA .1735, . 4069, .13568, .2218, ”8. 0 3/8 T” 

RH « RH / 4 
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' HEAT EXCHANGER PROPERTIES 
ID - . 95 * OD 
KM - 16 

' helium and water properties 

CPW - 4189 
PRH - .668 
CPH - 5193 
R * 2077 

' DATA INPUT 

INPUT "HELIUM MASS FLOW RATE (KG/SEC)"; MH 
INPUT "HELIUM INLET PRESSURE (MPa)"; P6 
P6 - P6 * 10 ^ 6 

INPUT "HELIUM INLET TEMPERATURE (C)"; T6 
T6 - T6 + 273.15 

INPUT "HELIUM OUTLET TEMPERATURE TO COMPRESSOR"; T1 
T1 - T1 + 273.15 
Q - MH ♦ CPH * (T6 - Tl) 

INPUT "HELIUM SPECIFIC PRESSURE DROP"; DP 

TOP: 

INPUT "INLET WATER VELOCITY (m/sec)"; VW 
GW * VW * 980 

INPUT "INLET WATER TEMPERATURE (C)"; TWI 

TWI - TWI + 273.15 

EFF - (T6 - Tl) / (T6 - TWI) 

PRINT "WHAT IS THE DESIRED WATER FLOW RATE (KG/SEC) "; 

INPUT MW 

CWCH » MW * CPW / MH / CPH 

TWO - ((T6 - Tl) / CWCH) + TWI 

PRINT "EFFECTIVENESS - "; EFF 

NTU - -CWCH * LOG (1 + LOG(l - EFF) / CWCH) 

PRINT "NTU NTU 

UA - NTU * MH ♦ CPH 

TWAVG - .5 * (TWO + TWI) - 273.15 

PRINT "AVERAGE WATER TEMPERATURE IS"; (.5 * (TWO + TWI)) - 273.15 

GOSUB VISCWATER: ' FIND VISCOSITY OF WATER 

PRW - INTERP (TWAVG, WPR()): ' FIND PRANDTL NUMBER 

IF PRW - 0 THEN INPUT "PRANDTL NO."; PRW 

KW - INTERP (TWAVG, WK ( ) ) 

IF KW - 0 THEN INPUT "THERMAL CONDUCTIVITY (W/mK) " ; KW 
PRINT VISW, PRW, KW 
TAV - .5 * (Tl + T6) 

DENSH * P6 / R / TAV 

VISH - (6.7 + .044 * TAV) * .000001 

GH * ((.006 * P6 * DP * DENSH) / ( . 02 NTU * .7631)) ^ .5 

' START LOOP 
DO 

REH - (GH ^ 4 ♦ RH) / VISH 
FH « FI * REH ^ (-F2) 
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JH » J1 * REH ^ (-J2) 

HH - (JH * GH * CPH) / (PRH ^ (2/3)) 

PRINT ”H he - HH; "W/m2K”, ”G he GH; "KG/m2S*' 

PRINT *’Re he - REH, ”JH - ” ; JH 

HW - (.023 ^ KW ♦ (PRW ^ .4) * GH ^ .8) / ( (VISW .8) * (ID ^ .2)) 

PRINT ”HW - HW 
IF LF <> 0 THEN 

M - SQR((2 ^ HH) / (KM * DEL)) 

ML « M * LF 

NF - ((EXP (ML) - EXP (-ML)) / (EXP (ML) + EXP (-ML))) / ML 

NOF - 1 - AFA * (1 - NF) 

ELSE 

NF » 0 
NOF - 1 
END IF 

A - (OD - ID) * .5 
AWA - 1 - AFA 

U - 1 / ( (1 / HH / NOF) + (A / KM / AWA) + (1 / HW / AWA) ) 

PRINT ”FIN EFF. NF, "U U 

AH - UA / U 

PRINT *‘Q Q 

VOL - AH / ALFA 

AC - MH / GH 

AFR - AC / SIGMA 

L - VOL / AFR 

DELP - FH * AH * GH * GH / 2 / P6 / AC / DENSH 

GH - GH ^ ((DP / DELP) ^ .5) 

LOOP UNTIL ABS((DELP - DP) / DP) < .05 
END OF LOOP 

NT * INT(4 * MW / (GW * PIE * ID * ID) + .5) 

PRINT "OUTPUT TO SCREEN (1) OR PRINTER (2)”; 

INPUT 1% 

IF 1% - 2 THEN U$ - "LPTl:" ELSE U$ - *'SCRN:" 

OPEN U$ FOR OUTPUT AS #1 

#1, STRING$(40, "*•') 

#1, "PRECOOLER DESIGN PROGRAM" 

#1, 

#1, "HELIUM TEMPS: INLET T6 - 273.15; "*C", "OUTLET -"; Tl - 273.15; "’C" 

#1, "HELIUM INLET PRESSURE P6 * .000001; "MPa" 

#1, "HELIUM MASS FLOW MH; "KG/SEC" 

#1, "WATER: INLET TEMP TWI - 273.15; "*C", "OUTLET TEMP •*"; TWO - 273.15; "’C" 

#1, " MASS FLOW MW; "KG/SEC "; MW * 15.8821; "GPM" 

#1, "HEAT EXCHANGER SURFACE: S$ 

#1, "RATIO OF HEAT CAPACITIES, (CW/Ch) -"; CWCH 
#1, 

#1, "CROSSFLOW HEAT EXCHANGER" 

#1, STRING$(40, "*") 

#1, " RESULTS OF CALCULATIONS " 

ND - INT(L / DT + .5) 



PRINT 

PRINT 

PRINT 

PRINT 

PRTNT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 
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NW - INT(NT / ND + .5) 
W - NW WT 
LT - AFR / W 



PRINT 


#1, 


"CROSSFLOW HEAT EXCHANGER DIMENSIONS" 


PRINT 


#1, 


"VOL VOL; "m3" 




PRINT 


#1, 


tt tubesheet 


II 


PRINT 


#1, 


"# OF TUBES NT, 


"TUBE LENGTH LT; "m" 


PRINT 


#1, 


"WIDTH W; "m". 


"DEPTH -"/ L; "m" 


PRINT 


#1, 


•t FLOW 


II 


PRINT 


#1, 


"PRESSURE DROP 


DELP, "FRONTAL AREA ; AFR, 


PRINT 


#1, 


"GEOMETRY RATIO - " 


; L / SQR(AFR) 


PRINT 


#1, 


"EFFECTIVENESS 


EFF, "NTU NTU 


PRINT 


#1, 


STRING$(40, "*") 





CLOSE 

PRINT "DO YOU WANT TO TRY AGAIN WITH THE SAME EXCHANGER AND HELIUM CONDITIONS 
INPUT N$ 

IF ASC(N$) - 121 OR ASC(N$) - 89 GOTO TOP 

CLOSE 

END 

VTSCWATER: 

f ★★★★★ FINDING THE VISCOSITY OF WATER 

IF TWAVG <- 100 OR TWAVG > 0 THEN 
IF TWAVG <-20 THEN 

VISW - .001 * 10 ^ (1301 / (998.333 + 8.1855 * (TWAVG - 20) + .00585 * 

^ 2) - 1.30233) 

ELSEIF TWAVG <- 100 THEN 

VISW - .001 * 10 ^ ((1.3272 ^ (20 - TWAVG) - .001053 (TWAVG - 20) ^ 
105)) ★ 1.002 
END IF 
ELSE 

PRINT "WATER VISCOSITY (PaS 10^-5)"; 

INPUT VISW 

VISW - VISW * .00001 
END IF 
RETURN 

END 

FUNCTION INTERP (T, A()) STATIC 

IF T < 0 OR T > 100 THEN INTERP - 0 
IF T - 100 THEN INTERP - A (5): EXIT FUNCTION 
IF T - 0 THEN INTERP - A(0) : EXIT FUNCTION 
U - T / 20 
A% - INT(U) 

INTERP - (U - A%) (A(A% + 1) - A(A%)) + A (A%) 

END FUNCTION 



(TWAVG - 20) 



) / (TWAVG + 
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Appendix C Heat Exchanger Surfaces Characteristics. 

The following table gives the characteristics of the heat exchanger surfaces in the sur- 
face data base for COMPHX.BAS. 

Table C-1. Heat Exchanger Surfaces. [1 1] 



FIN 


PLATE 


FIN 


FINS/in 


FIN 


HYD. 


BETA 


AREA 


TYPE 


SPACING 


THICKNESS 


LENGTH 


RADIUS 




RATIO 




in 


in 




in 


in 


ftVcu ft 


ftVft 


PLAIN FINS 














1- 1 


0.500 


.00600 


16.00 


0.000 


.02535 


423 


.897 


1- 2 


0.750 


. 03200 


3.01 


0.000 


.10638 


98 


.706 


1- 3 


0.750 


. 03200 


3.97 


0.000 


. 08460 


119 


.766 


1- 4 


0.470 


.00600 


5.30 


0.000 


.06048 


188 


. 719 


1- 5 


0.405 


. 01000 


6.20 


0.000 


. 05460 


204 


.728 


1- 6 


0.823 


. 00800 


9.03 


0.000 


.04566 


244 


. 888 


1- 7 


0.250 


. 00600 


11.10 


0.000 


.03036 


367 


.756 


1- 8 


0.480 


. 00800 


11.11 


0.000 


.03459 


312 


. 854 


1- 9 


0.330 


.00600 


14.77 


0 . 000 


.02544 


420 


.844 


1-10 


0 . 418 


. 00600 


15.08 


0.000 


. 02628 


414 


.870 


1-11 


0.250 


.00600 


19.86 


0.000 


.01845 


561 


. 849 


1-12 


0.544 


.01000 


10.27 


0.000 


.03777 


290 


.863 


1-13 


0 . 249 


. 00600 


11. 94 


0.000 


. 02820 


393 


.769 


1-14 


0.250 


.00600 


12.00 


0.000 


. 02823 


393 


.773 


1-15 


0.256 


.00600 


16.96 


0.000 


.01695 


608 


. 861 


1-16 


0.204 


. 00600 


25.79 


0.000 


.01131 


856 


. 884 


1-17 


0.345 


. 00600 


30.33 


0.000 


.01203 


813 


. 928 


1-18 


0 . 100 


. 00200 


46.45 


0 . 000 


. 00792 


1333 


. 837 


LOUVERED 


FINS 














2- 1 


0.250 


. 00600 


6.06 


0.375 


. 04380 


256 


. 640 


2- 2 


0.250 


. 00600 


6.06 


0.375 


. 04380 


256 


. 640 


2- 3 


0.250 


. 00600 


6.06 


0.500 


. 04380 


256 


. 640 


2- 4 


0.250 


. 00600 


6.06 


0 . 500 


. 04380 


256 


. 640 


2- 5 


0.250 


. 00600 


8.70 


0.375 


. 03588 


307 


.705 


2- 6 


0.250 


.00600 


8.70 


0.375 


. 03588 


307 


.705 


2- 7 


0.250 


. 00600 


11.10 


0.188 


.03036 


3 67 


.756 


2- 8 


0.250 


.00600 


11.10 


0.250 


.03036 


367 


. 756 


2- 9 


0.250 


.00600 


11.10 


0.250 


.03036 


367 


.756 


2-10 


0.250 


. 00600 


11.10 


0.375 


.03036 


367 


. 756 


2-11 


0.250 


. 00600 


11.10 


0.375 


. 03036 


367 


.756 


2-12 


0.250 


.00600 


11.10 


0.500 


. 03036 


367 


.756 


2-13 


0.250 


. 00600 


11.10 


0.750 


.03036 


367 


.756 


2-14 


0.250 


. 00600 


11.10 


0.750 


. 03036 


367 


.756 
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STRIP FINS 



3 - 1 


0.250 


. 00600 


11 


3 - 2 


0 .485 


. 00400 


12 


3 - 3 


0.414 


. 00600 


15 


3 - 4 


0.375 


.01000 


13 


3 - 5 


0.237 


. 00600 
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Appendix D Heat Exchanger Surface Performance Data. 




Fin pitch = 46.45 per in = 1829 per m 

Plate spacing, ^ = 0.100 in = 2.54 x 10"^ m 

Fin length flow direction = 2.63 in = 66.8 x ICT^m 

Flow passage hydraulic diameter, 4r^ = 0.002643 ft = 0.805 x 10"^ m 

Fin metal thickness = 0.002 in, stainless steel = 0.051 x 10"^ m 

Total heat transfer area/volume between plates. 0 - 1332.45 * 4372 m^/m^ 

Ftn area/total area = 0.837 



Figure D-1. 



Plain-Fin plate-fin surface 46.45T. 
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Fin pitch = 24. 12 per in = 950 per m 
Plate spacing, b - 0.075 in = 1 .91 x 1 0*^ m 
Fin length = 0. 1 1 1 m * 2.8 x 1 0’^ m 

Flow passage hydraulic diameter, 4r^ - 0.003966 ft = 1.209 x 10'^m 
Fin metal thickness = 0.004 in = 0.102 x IC^m 

Total heat transfer area/volume between plates, P = 862.7 ft^/ft^ = 2,830 m^/m^ 
Fin area/total area = 0.857 



Figure D-2. 



Strip-fin plate-fin surface 1/9-24.12 
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Tube cxitside diameter * 0.375 in * 9.525 x 10"^ m 
Hydraulic diameter, * 0.0196 ft = 6.071 x 10"^m 
Free-flow area/frontal area, o = 0.333 

Heat transfer area/total volume, a - 67. 1 ft^/ft^ = 220.144 m^/m^ 
Note: Minimum free-flow area is in spaces transverse to flow. 



Figure D-3. 



Flow normal to a staggered tube bank, surface S 1.50-1.00 
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Tube outside diameter = 0.250 in = 6.35 x 10"^m 
Hydraulic diameter, 4r^ = 0.0166 ft = 5.029 x 10"^m 
Free-flow area/frontal area, o = 0.333 

Heat transfer area/totai volume, a * 80.3 * 263.451 m^/m^ 

Note: Minimum free-fiow area is in spaces transverse to flow. 



Figure D-4. 



Flow normal to a staggered tube bank, surface S 1.50-1.25 (s) 
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Tube outside diametef = 0.42 in * 1 0.67 x 1 0*^ m 
Fin pitch = 8,72 per in * 343 per m 

Flow passa 9 e hydraulic diameter, 4r^ = 0.01452 ft = 4.425 x 10'^ m 
Fin thickness (average) ^ = 0.019 in, copper = 0.48 x 10*^ m 
Free-flow area/fronial area, <7 = 0.494 

Heat transfer area/total volume, a = 136 ft^/ft^ = 446 m^/m^ 

Fin area/ioial area = 0.876 
f Fins slightly tapered. 



Surface CF-8. 72(c), finned circular tubes. Circular fins. 



Figure D-5. 
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Tube outside diameter = 0.402 in = 10.2 x lO'^m 
Fin pitch = 8.0 per in = 315 per m 

Flow passage hydraulic diameter, 4r^ = 0.01 192 ft = 3.632 x IC^m 
Fin thickness = 0.013 in = 0.33 x lO'^m 
Free flow area- frontal area, n - 0.534 

Heat transfer area/total volume, a = 179 fl^/ft^ = 587 m^/m^ 

Fin area/total area = 0.913 

Note: Minimum free flow area in spaces transverse to flow. 



Figure D-6. 



Surface 8.0-3/8T, Firmed circular tubes, continuous fins. 
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